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SUMMARY 

This report describes a FORTRAN coded computer program and method to 
predict the reaction control fuel consumption statistics for a three axis 
stabilized rocket vehicle upper stage. It uses a Monte Carlo approach which 
13 made more efficient by using closed form estimates of impulse useage. It 
includes effects of rocket motor thrust misalignment, static unbalance, 
aerodynamic disturbances, and deviations in trajectory, mass properties and 
control system characteristics. This routine has been used for over a decade 
to accurately predict the control fuel consumption statistics for the Scout 
launch vehicle second and third stage reaction control systems. 

By selection of input data and options this routine can be applied to many 
types of on-off reaction controlled vehicles. 

The psuedo random number generation and statistical analyses subroutines 
including the output histograms can easily be used for other Monte Carlo 
analyses problems. 

A typical run of 200 samples requires 2 seconds of central processor time 
on a CDC CYBER 175 computer. 
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1.0 INTRODUCTION 


Reaction control system fuel consumption predictions for missiles, launch 
vehicles and space vehicles is usually required to be accurate so that weight 
can be minimized without unduly sacrificing mission success probability. 

Almost all reaction controlled missiles, launch vehicles and spacecraft would 
result in loss of control and mission failure if the control fuel useage were 
seriously underpredicted. On many vehicles overprediction of fuel 
requirements results in costly weight and volume penalties. The control 
system designer requires an accurate prediction of the statistical 
distribution of required impulse or control fuel to size the system. For 
missions with different coast periods a prediction of fuel useage is needed 
for propellant loading. 

The method and computer program for predicting reaction control system 
fuel consumption statistics presented in this report has been used over the 
past fifteen years for the Scout launch vehicle second and third stage 
reaction control systems. Predictions using this method correlate well with 
flight data statistics. There is a slight tendency of the method to be 
conservative particularly at the low end of the probability distribution. 

This is due to an assumption of symmetric attitude limit cycle motion when 
undisturbed. In actuality a certain small percentage of vehicles achieve a 
near zero duty cycle in an axis rather than a symmetric limit cycle. As a 
result fuel consumption can be lower than the minimum predicted by the method. 

Prediction of control fuel by this method is only as good as the accuracy 
of the input data; i. e., "Garbage in - garbage out." Scout experience in 
accurate prediction results from accurate knowledge of disturbance statistics, 
particularly rocket booster motor thrust misalignment. These are derived from 
detailed post-flight data analyses. 
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2.0 METHODOLOGY 


The Monte Carlo method, used herein, is a relatively straightforward 
approach for predicting the statistical outcome of a highly non-linear process 
involving many variables. On-off reaction control system fuel consumption 
prediction for an upper stage rocket booster vehicle is a good candidate for 
this method. Large amounts of computer time are normally associated with the 
Monte Carlo approach. The method described herein is made efficient by using 
a series of closed form approximations to the fuel consumption during 
transient maneuvers and steady state operation. As a result several hundred 
samples can be computed with a few seconds of computer time on a high speed 
digital computer. Assumptions and equations used in the computer program are 
described in the following paragraphs. 


2.1 Fuel Consmption Assumptions 


Many assumptions are made to simplify the control fuel consumption 
prediction process. These include; 

. Uncoupled axes (pitch, yaw, roll) each having its own set of control 
motors. 

. On-off control has a simple deadband with constant switching slope for 
each phase of the flight (boost and coast) (see phase plane Figure l). 

. Iiftien undisturbed each axis control achieves- a symmetric limit cycle 
motion. 

. Disturbing torques and impulses are balanced by the control system 
without exceeding the zero rate deadband attitude error. 

. Disturbing torques of sufficient magnitude to cause deadband crossing 
results in fuel consumption based on angular impulse balance plus 
symmetric limit cycle motion (Figure l) . 

. Vlhen sets of reaction control motors are shared by two or more axes, 
each axis has full use of the motors in an uncoupled sense (e.g., 
yaw-roll mixing with four (4) jets is decomposed into a yaw-axis 
control with 2 motors in each sense and roll axes control with 2 motors 
forming a couple in each sense or direction). 

. Sensor, computational, filter, and control motor switching hardware 
responses are treated as a simple equivalent time transport lag. 

. Structural mode coupling with the control system is limited. An 
effective time transport lag is assumed to adequately account for 
structurally coupling. This is the case for the Scout launch vehicle 
(see Appendix D of Reference 1). 

. Reaction control motor pulse shapes can be expressed as an equivalent 
square wave with appropriate thrust level and time delays. 

. Aerodynamic stability derivatives are linear and can be expressed as a 
constant normal force coefficient slope with angle of attack and a 
constant aerodynamic center location. 
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. The vehicle is axisymetric in terms of moments of xnertia and 
aerodynamic configuration. 

. Trajectory parameters are as input and are not significantly perturbed 
by the off-nominal control behavior. 

. Changing control fuel consumption does not alter the vehicle mass 
properties. 

. Aerodynamic angle of attack includes the attitude error equal to the 
deadband halfwidth. 


2.2 Equations 


Equations describing the disturbances and fuel consumption are presented 
below. Trajectory geometry and definition of angles are presented in Figure 2. 


2.2.1 Attitude Control Acceleration 


Control angular accelerations are: 


( 2 - 1 ) 


12I„ 


( 2 - 2 ) 


where , 


<f>c = 


2F„R 


c c 


121 


*^6 ° [(Xg- Xcg) COSX+ Zp SIN M/12 
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2.2.2 Disturbances 


Disturbing accelerations are due to booster thrust misalignment > 
static unbalance and aerodynamic torques. 
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Angles of attack and sideslip are; 


(2-6) g = 0pR - TpR - Tel + “w + ^“db 

(2-7) P = /^PR fe1 '>'PR ^db 


(see Hgure 2) 


where Vq is the perturbation flight path angle from previous stages plus 
^ the integrated change in flight path angle due to pointing error 
(deadband) . 



the equivalent flight path deviation in the yaw plane. 
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Angles of attack and sideslip due to wind are; 


(2— lO) ^ 

V.V„C0S7pRC0S(r„-fpR) 


(2-11) p = 


57.3 


VwSiN(ry,-rpR) 


V>Vv,C0S7prC0S(?„-?Pr) 


- 4 - 


Incremental angles of attack and sideslip due to pointing errors are: 


(2-12) - ±d-(y^~ y^j)/2 

(2-13) = ±d + (f,-fei)/2 


2.2.3 Control System Duty Cycles 


Control system duty cycles during boost include those due to disturbance 
torq.ue balance and undisturbed limit cycle behavior. 


(2-14) 

(2-15) 


DCd = 


; 0n > 0 


TEST 


° (LIMIT CYCLE) ' I^dI <®TEST 


Test accelerations are those disturbing accelerations which define the 
boundaiy between a one-sided motor firing and crossing the deadband as shown 
in Figure 1 and 3* 


2.2.4 Limit Cycle Motion 


Undisturbed limit cycle motion results in a duty cycle dependent on 
control system characteristics. 


(2-16) 'J*^P(LIMIT CYCLE) ° 


where , 
(2-17) 
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for square wave pulse of control motor 
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for thrust with overshoot (Figure 4 ) 
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With control motor thrust overshoot the above equations for duty cycle are 
computed with one iteration starting with = 1. 

When control acceleration is veiy high, long delay times can result in a 
deadband overshoot condition; i. e., the error signal crosses the total 
deadband before the opposite motor actually turns off. The control 
acceleration at which this occurs is; 


( 2 - 21 ) 


MAX 


dt2 Kp/Kp - T, - T 2 - H (Kp/Kp - T2/2)l 
T2(Kp/Kp-T2/2)(Kp/Kp-T,/2) 


If 0- > 0 , the duty cycle becomes, 

c Ci 


(2-22) 

where, 

(2-23) 
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Equations for yaw and roll duty cycles are the same. 


2 . 2.5 Fuel Consumption 


Control fuel consumption is dependent on the total impulse and effective 
specific impulse of the control motors. The specific impulse of the control 
motors used for pitch and yaw is assumed to be the same. A separate roll 
motor specific impulse can be used. In most cases the specific impulse during 
coast IS different than during boost. Boost fuel consumption is: 


(2-24) WbOOST “ ''''B00STp^''''B00STy^''''B00STp 

For pitch and yaw the fuel consumed is: 


(2-25) Wgpo5Tp = ^ DCp dt + lp^pj 

where ^cap additional impulse needed to "capture" at 

separation/ignition from the initial conditions on rate 
and displacement. 


(2-26) 


IcAPp 



2(go-d) 

Kr/Kq 


+ 2 



- 6 - 



Roll fuel consumption Includes an amount for booster roll torque (RTl) and 
pitch and yaw motor misalignment. 


(2-27) 
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During the coast phase a symmetric limit cycle is assumed. Coast fuel 
consumption also includes (l) an incremental aunount for transients due to 
deadband reduction from boost to coast system, ( 2 ) guidance progreun attitude 
rate changes in pitch and yaw, and ( 5 ) an optional retro maneuver with the 
boost pitch and yaw control motors. 
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This assumes that the four boost pitch and yaw motors fire continuously 
and the system operates at a 100 percent duty cycle. Additional fuel for 
deadband reduction and torquing rate changes for each axis is 


(2-30) 
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where A d 
^b 

<ib»^c 
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•• 


13 the change in guidance pitch or yaw program rate 
(deg/sec) 

is the boost system limit cycle rate (rad/sec) 

is the boost and coast system deadband halfwidth 
(degrees) 

is the coast system rate to displacement gain ration 
(sec) 

is the coast system control acceleration (rad/sec^) 
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2.2.6 Trajectory and Vehicle Characteristics 


Time histories of pertinent trajectory and vehicle characteristics are 
input to use in the calculations. These include dynamic pressure, velocity, 
booster thrust and propellant weight remaining, altitude, flight path angle, 
and azimuth. Wind velocity versus altitude is also input. Thrust 
misalignment is assumed to be described as an initial value and a slope which 
IS correlated (see section 2.3)* 

Mass properties (center of mass station and transverse moment of inertia 
is assumed to be a second order polynomial of fraction of booster propellant 
consumed. 

Given a three point curve (x^, y^), (x 2 > 72 ) and (xj, 75 ) the 
second order function is. 


(2-31) V =■ + K^X + K 2 X^ 


where, 

(2-32) 


K2 = 


V3-yi-^(y2-v,) 

22 

(Xj^ - X2X3) 


and, 

(2-33) K, = 


<V2-Yi> 

X, 


-K2X2 


The independent variable (x) is the fraction of booster propellant 
consumed, and the dependent variable is either center of mass station or the 
moment of inertia. The three input points should include the 0, 0.5 and 1.0 
fractional propellant consumed values for best results. The burnout, or 1.0 
fractional propellant consumed, values are used for coast calculations. 


2.2.7 Optional Trajectory Characteristics 


An optional input for certain trajectory calculations is available to 
approximate dynamic pressure, altitude and flight path angle. This option is 
used primarily where many parametric runs are required for which there are no 
trajectories. Given the initial flight path angle (Tq)» velocity (Vq ) and 
dynamic pressure (Qq ) the density at ignition is calculated. 




- 8 - 


(2-34) 



This density is used with a standard atmosphere logarithmic density versus 
altitude data set to obtain initial altitude* Thrust, weight, and drag 
coefficient (Cj)S) are input and used to propagate a gravity turn trajectory. 

The set of first order equations used to propagate the trajectory is; 


(2-35) 

y _ 9o<T-CqSQ) 

W 

(2-36) 

h»-57.3 9 COS 7 

(2-37) 

7 = V SIN 7 

where. 


(2-38) 

Q =■ J4pV^ 

(2-39) 


and. 



= 32.1 74 ft/sec^ 
Rg - 20,919,668 ft 


A Runge-Kutta integration subroutine is used to integrate these 
during boost to obtain the nominal trajectory parameters (V, h, T » 


equations 

Q)- 


2*3 Statistical Analysis 


The analysis method used in this routine is the Monte Carlo Technique. It 
is based on random sampling of a set of random input variables for the 
vehicle, trajectory, disturbances, and control system and computing the 
resulting total fuel consumption, boost fuel, coast fuel, and retro time. By 
computing a large number of such flights the statistical variation of fuel 
consumption is obtained. The method of treatment of input variables and 
analysis of output is listed below. 


2.3.1 Input Random Variables 


There are many random variables required as input. These variables are 
assvimed to be Gaussian (Normal Distribution Function). A mean and standard 
deviation is entered for each of these variables. In all cases except booster 
thrust misalignment each variable is uncorrelated with all other variables. 
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For uncorrelated variables the random sample is: 

(2-40) X = X + RND 

where x is the mean value of x 

<7^ is the standard deviation of x 

RHD IS a random normal deviate (computed by routine) 

Where there is a correlation between two variables as in the case of 
thrust misalignment the random samples are, 

X = X + RND,Ojj 


(2-41) 


V = y + RND,ay + R N D 2 


'xy 


where 

N 


is the correlation coefficient of the x and y variables 
is the number of samples used to derive the correl-ition 
coefficient 


The following is a list of random variables for input; 


4,<Ao.0o 

>'eXe 

"^FRACT 


eg 

SU ,SU 

P y 

RTI 

^JET 

^w 

F 

c 

^sp 

^p/^c 

t 

p 

^Im 

^2m 

^f 

^RG 

Kjj/Kjj 

d 

H 


initial value of booster thrust misalignment 

f 

slope of booster thrust misalignment with time 

initial pitch, yaw and roll rate 

initial attitude error in pitch, yaw, and roll 

initial flight path and azimuth error 

dynamic pressure fractional deviation (Q/Q ) 

transverse moment of inertia 

center of mass station 

pitch and yaw static unbalance 

booster roll torque impulse 

boost control motor misalignment angles 

wind velocity versus altitude 

wind direction 

control motor force levels 

control motor specific impulse, boost and coast 

pitch and yaw control motor overshoot ratio 

pitch and yaw control motor time from on to peak thrust 

control motor effective turn on time delay 

control motor effective turn-off time delay 

autopilot filter effective time delay 

autopilot rate gyro effective time delay 

autopilot rate to displacement gain ratio 

control system deadband halfwidth 

switching hysteresis fraction 
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other input variables are not treated as random. 


2.3.2 Output Statistical Analysis 


The routine calculates boost, coast, and total fuel consumption for N 
samples (flights). The output includes the values of fuel consumption with 
the calculated discrete probability levels for the N samples. In addition to 
these discretes, the sample population mean, standard deviation, skewness, and 
kurtosis are output for the parameters and the logarithm of the parameters. 

The probability distribution function parameters and key probability levels 
are based on three distribution functions; (l) Normal, ( 2 ) Log-Normal, and (3) 
Weibull. An equivalent level of significance for each of these three 
distributions functions is estimated and output baaed on the CHI-Squared 
GJoodness of Fit Test. 

When the fuel consumption samples are arranged in ascending order, the 
probability that the fuel consumption will be under a particular value is 
given by the following equation: 

N 

(2-42) P (fuel consumption< Wi ) = — ! — 

N + 1 


whe re , W^ 
Ni 

N 


is a particular value of fuel consumption. 

is the sample number corresponding to ¥3^ in the 
arranged fuel consumption distribution. 

is the total number of samples calculated. 


It IS desirable to obtain confidence intervals about a given probability 
level. This is also a function of the sample size. The two sided confidence 
interval for the standard deviation ( tr ) of a population of sample size N 
extracted from a Normal Distribution Function is provided by the following 
equation: 


( 2 - 43 ) 



where , 

Sjj 13 the sample standard deviation calculated from N 

samples. 

a is the true standai^ deviation. 

ky is the number of standard deviations on a Normal 

Distribution which provide a probability level equal to 
the confidence coefficient (i.e., 0.95 confidence - 
1.96a , ky = 1.96). If only one side of the confidence 
interval (high side) is desired, the values of ky 
become based on cumulative Normal Distribution 
probabilities. 
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The confidence limits for other probability levels can be established 
assuming normality. That is, once the standard deviation is established with 
an upper confidence limit (say 95 percent) ^then the value at another 
probability level can be expressed as a number of standard deviations from the 
mean. 


(2-44) 



(95^ confidence ky = 1.645 for one 
sided upper limit of confidence 
interval. ) 


where y’ is the value of parameter 'y* at a probability level 

equivalent to the mean plus 'K* standard deviations on a 
Normal Distribution with the desired confidence 
coefficient. 

y is the sample mean (assumed to be the true mean value of 

the distribution function) . 


K is the number of standard deviations required to provide 

the desired probability level. 

The following is a short table for values of 'K' and 'k' for various 
probability levels and confidence coefficients. 


Probability Level 

K, ky (one-sided) 

0.500 

0 

0.900 

1.282 

0.950 

1.645 

0.995 

2.575 

0.999 

5.090 


For example, to find the value of 'y' at 99*5^ probability at 95 percent 
confidence, given sample mean, y, and variance, si , from N samples 


(2-45) 


2.575 s. 


V = V + 




In the case of fuel consumption, the distribution function is not usually 
Gaussian (Normal Distribution). Therefore, several approximate distribution 
functions are tested in the routine based on the Chi -Squared Goodness of Fit 
Test. These include the Normal, Log-Normal and three parameter Weibull 
distribution. The Log-Normal Distribution is simply a transformation of the 
'y' parameter to the logarithm of 'y'* As shown in the sketch below, if the 
probability density fxinction is one sided, the density function of the log 'y* 
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may approximate a Normal Distribution. If this occurs the properties of the 
Normal Distribution function such as probability level and confidence 
coefficient, can be applied to log 'y'* 

Log-Normal Density Function 




The three parameter Weibull distribution is also used. This distribution 
function is described in detail in Reference (2). The cumulative probability 
distribution function and density function can be described in terms of three 


parameters 'a' 

, 'b' , and 'c 

' by the following equations. 

(2-46) 

f(y) » - Ae'* 
b 

density function 

2-47) 

F(v) = 1 - 

distribution function 

where, 

‘ ■ (^)' 


and where 

• « * • V • * • 

SL f Of C 

are values greater than zero. 


The probability density function can take on many shapes by varying the values 
of 'a', 'b' , and *c' and thus has application to a wide variety of problems. 
For this application it will have a shape similar to that shown in the sketch. 

Weibull Density Function 
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For this distribution function the value of 'y' at a given probability 
level is given by: 


(2-48) V = a - h [In II - P)] 

where P is the probability level. 


The routine prints out the parameters 'a', 'b', and 'c' determined from 
the fit based on McClintock' s method of Reference (2). It also prints out the 
values of the parameter at several key probability levels. Calculation of 
confidence intervals for this distribution function is not included in this 
routine . 

In addition to the statistical analyses described above, histograms of the 
output are made on the line printer for ready comparison. 
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3.0 PROGRAM DESCRIPTION 


3.1 General 


This routine is programmed xn FORTRAN IV for the GDC CYBER 175 system. 
The only non- ANSI code used xs the pseudo random number generator (RANP). An 
ANSI option for a replacement is identified in the program Ixsting. The 
routine requires approximately 16 K words. Program flow, and user 
instructions are presented in the following paragraphs. 


3.2 Program Flow 


The program flows straight forward in five basxc parts. 

. Input optxon selection 
. Input of all data 

. Precalculation of boost phase txme histories at 
xntegration steps 

. Monte Carlo calculation of specified number of cases 
of fuel consumptxon or coast txme 
. Statxstical analysis of results 

Flow charts showing the general sequence are presented in Figure 5* A 
list of subroutines, their cross references and common blocks are presented in 
Figure 6. A description of the subroutines is given in the following 
paragraphs. 


3.3 Subroutine Descriptxon 


Most of the repeated specialized calculations are performed in the 
seventeen subroutines. Subroutxne and common reference map is presented xn 
Figure 6. The descriptions of each follow in alphabetical order. 

ALTI 


’ALTI’ xs used by subroutxne TABGEN when the xnput option requires 
generatxon of certain trajectory parameters. Gxven the dynamxc pressure and 
velocity it computes the altitude based on the logarithm of inverse densxty 
versus altxtude table (RHO) supplied in labeled common "DRH" defxned in the 
'TABGEN* subroutine. 

The call to the subroutine is: 

CALL ALTI (ALTO, VO, QO) 

alto - is the altxtude in feet computed by this subroutine. 

VO - is the velocxty xnput (ft/sec). 

QO - is the dynamxc pressure xnput (lbs/ft2). 
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ASCEND 


This subroutine rearranges an array of ntjmbers in monotonically increasing 
or descreasing order. 

The call to this subroutine is: 

CALL ASCEND (L, VAL, M) 

L - is the number of values in array VAL (input). 

VAL - is the array to be rearranged (input and output) . 

M - input option for ascending or descending order, 

M = 0 ascending 
M = 1 descending 

CALCU 

This subroutine is used with "RUNGE" to calculate the trajectory state 
variables and their first derivatives in TABGEN. It is based on equations 
2-35 through 2-39» 

The call to this subroutine is: 

CALL CALCU 

Input and output is passed through labeled common blocks 'DDI', 'DDY', and 
'DHH'. 'DDY* contains the state variables and time information. 

CHISQ 


This subroutine is based on the Chi-Squared Goodness of lit Test to a 
known distribution. It tests for level of significance against a Normal, Log 
Normal, and a Weibull distribution function. 

The call to this subroutine is: 

CALL CHISQ (N, A, CR, NG, AM, AS, ¥A, WB, VC, LD) 


N 

A 

CR 

NG 

AM 

AS 

WA 

WB 

VC 


is the total sample size input. 

IS the input array of values. 

is the output computed level of significance. 

is the number of cells for the CHI -Squared Goodness of 

Pit Test. 

is the computed mean, 
is the computed standard deviation, 
is the Weibull 'a' parameter, 
is the Weibull *b' parameter, 
is the Weibull 'c' parameter. 


LD = 0 for a two parameter distribution (the Normal and Log-Normal) . 
LD = 1 for a three parameter distribution (Weibull). 


Tables of the cumulative Normal probability distribution *C' versus number 
of standard deviations 'Z' are embedded in this subroutine and are part of 
labeled common block 'CZ'. 


- 16 



CYCLE 


'CYCLE' computes the duty cycle for the on-off control during symmetric 
undisturbed limit cycle operation. It is based on equations 2-16 through 
2-20. It Includes the option for deadband overshoot and overshoot type 


control motor 

thrust 

response. 

The call 

to this 

subroutine is: 

CALL CYCLE (THEC, 

GR, DR, H, Tl, T2, DC, A, NDOVER, FCT) 

THEC 


is the control acceleration (rad/sec^). 

GR 


is the rate to displacement gain ratio (Kg/K]) sec). 

DR 

- 

is the deadband halfwidth (rad). 

H 


is the switching hysteresis fraction. 

T1 

- 

is turn-on effective delay time (sec). 

T2 

- 

is turn-off effective delay time (sec). 

DC 

• 

is the computed duty cycle. 

A 

- 

is the computed limit cycle rate (rad/sec). 

NDOVER 


is a deadbeind overshoot condition check 
= 0 no deadband overshoot 
= 1 deadband overshoot 

FCT 


is a factor computed by subroutine OVRF used when second 
order thrust response is indicated. If FCT is less than 
or equal to zero iteration on pulse width is not used. 


HISTO 


This subroutine sorts and plots a histogram on the line printer from an 
array of numbers. It must be given maximum and minimum values, the number of 
values and the increment for each cell of the histogram. The maximum cell 
count is 94 to fit the line printer. Below each cell the numerical value of 
the cell range and the count is printed. Asterisks are used to build the 
histogram. 

The call to this subroutine is: 


CALL HISTO (X, N, XMAX, XMIN, DX) 


X 

N 

XMAX 

XMIN 

DX 


is an array of numbers for the histogram (input), 
is the number of values in the 'X' array (input), 
is maximum value value in *X' (input), 
is minimum value in *X' (input), 
is the histogram cell width (input) . 


OVRF 

This short subroutine computes a factor for control motor thrust when 
characterized by an overshoot response. 


The call to this subroutine is; 


CALL OVRF (T, OR, FAC) 

T - is given time from zero thrust to first overshoot peak 

(sec) . 

OR - IS given overshoot ratio (Fp/Fc). 

FAC - is computed factor to be used for thrust level on short 

pulse widths. 
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PAGEHD 


For output this subroutine ejects a page and prints run number and page 
number at the top of each page. 

The call to this subroutine Is: 

CALL PAGEHD (NRUN, NPAGE, NLINE) 

NRUN - 13 the Input run number to be printed. 

NPAGE - Is the Input page number to be printed. 

NLINE - Is the returned line number count for the main routine. 


PZC 

This function subprogram computes the probability of a parameter falling 
In a specific range for a Normal or Welbull distribution function. It uses a 
table lookup of probability for the Noimal Distribution function. 

The call to this subroutine Is: 


CALL PZC (XL, XU, AH, AS, WA, WB, WC, LD, M) 


XL 

XU 

AM 

AS 

WA, ¥B, WC 
LD 


M 


Is Input lower limit value. 

Is input upper limit value. 

Is Input distribution mean value, 
is input distribution standard deviation, 
is input Welbull distribution parameters if LD > 0. 
is option 

LD < 0 Normal distribution. 

LD > 0 Welbull distribution. 

is current or last table lookup indicator. 


This subroutine uses the Nomal distribution tables in common block 'CZ* 
entered in Subroutine CHISQ. Probabilities for the Welbull distribution are 
based on Equation 2-47. 


RANGE 


'RANGE* is used for histogram preparation. It computes a desireable cell 
definition using the range of the data and a desired number of cells. The 
number of cells will be equal to or greater than that specified depending upon 
the evenness of the cell width. Acceptable engineering divisions are used in 
selecting cell starting and ending values. The array must be in ascending 
order. If it is in descending order it is first rearranged to ascending order. 
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The call to this subroutine is: 

CALL BAN6E (X, N, XHAX, XMIN, LX, NLX, ERRTB) 


is an input array of values in ascending or descending 
order. 

is number of values in 'X' array, 
is calculated upper value of *X' array, 
is calculated lowest value of 'X* array, 
is computed cell width for histrograms. 
is the computed number of cells for the histrogram. 
is an error signal generated. 

- 0 indicates no problems 

>■ 1 Indicates a zero range for the 'X* table 
RNLX 

A pseudo random number generator for a Noxmal Ustribution, N(0, l). It 
uses as computer supplied uniform random number generator (RAITF) to generate a 
number from zero to one. This is transfoimed to give the ^ndom Normal 
Deviate (EXD) from the Normal distribution having zero mean and unity 
variance. The coding presented in this report also shows an alternate pseudo 
random number generator which can be used in place of RANP. 

The call to this subroutine is: 

CALL RNLX (JK) 


X 


N 

XMAX 

XMIN 

LX 

NLX 

ERRTB 


JK 


is the number of random normal deviates to be calculated. 


The input seed *K' and the generated random normal deviates 'RXD' are 
transferred through blank common. 


RUNGE 


This is the basic fourth-order RUNGE-KUTTA integration method wMch is 
used in conjunction with CALCU to integrate the trajectory equations. 

The call to this subroutine is: 

CALL RUNGE (L, l) 

L 

I 

TABGEN 


is the control parameter for the Runge-Kutta Integration 
formulae. 

is the counter for the number of passes made in the 
multi-step process of one integration step. It is set 
to zero after completion of a complete cycle. 


In cases where a ballistic (gravity turn) trajectory is flown following 
booster ignition certain trajectory tables are estimated by this subroutine. 

It uses a fourth order RUNGE-KUTTA integration of the trajectory equations in 
subroutines CALCU and RUNGE. It requires initial conditions of flight path 
angle, dynamic pressure, velocity, predicted attitude error, and drag 
coefficient. Input tables of thrust and weight are used. After filling 
tables of dynamic pressure, velocity, altitude and flight path angle, it 
changes the input pitch attitude error table (TABLE 6) to pitch attitude using 
the newly calculated flight path angle as the reference attitude. 
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The call to this subroutine is: 

CALL TABGEN (QO, VO, GAMO) 

QO - is the input initial dynamic pressure (ibs/ft^). 

VO - is the input initial velocity (ft/sec). 

GAHO - is the input initial flight path angle (degrees) 

Output data is transferred via labeled common block 'LOUT*. Input data is 
transferred by common blocks 'DHH' , 'DDI', and *DDY'. This subroutine 
contains the data table of the logaxrithm of inverse atmospheric density versus 
altitude 'HHO* which is contained in common blocks 


TBLH 

This is a single table lookup subroutine using linear interpolation 
between points. This subroutine requires separate arrays of abscissas and 
ordinates. The abscissas must be in ascending order. 

The call to this subroutine is: 

CALL TBLN (Y, X, T, A, NT, M) 


Y 

X 

T 

A 

NT 

H 


is the ordinate to be found. 

is the given abscissa. 

is the abscissa table. 

is the corresponding ordinate table. 

is the number of values in each table. 

is a current locator for the search of the table. 

'M'must be greater than zero and leas than or equal to 

'NT* . M returns the current location found for the 

abscissa and should be used for the next lookup of the 

same table to reduce the search time. 


TBLU 

This is also a single table lookup. It is based on linear interpolation 
between points for a single array having alternating values of abscissas and 
ordinates. The abscissas must be in ascending order. 


The call to this subroutine is: 


CALL TBLU (NT, Y, X, T, M) 


NT 

Y 

X 

T 

M 


number of values in table 'T' including abscissas and 
ordinates. 

is the ordinate to be found, 
is the given abscissa. 

is the table of alternating abscissas and ordinates, 
is the table locator described under 'TBLN* above. 
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THEMIN 


'THMIN* computes the minimum acceleration for which one sided limit cycle 
motion can occur (ilgure 5). 


The call to this subroutine is: 


CALL OHEMIli (THEC, GR, DR, H, Tl, T2, THMIN) 


THEC 

6R 

DR 

H 

Tl 

T2 

THMIN 


is the given control acceleration (rad/sec^), 

is the rate to displacement gain ratio (Kj^/Kj) sec). 

is the deadband halfwidth (radians). 

is the switching hysteresis fraction. 

is the tum-on time delay (seconds). 

is the turn-off time delay (seconds) . 

is the computed disturbing angular acceleration at which 
the deadband will get crossed (rad/sec^). 


VBL 

'VBL* is the main subroutine which computes the statistical parameters 
from an array of numbers. It computes the mean, standard deviation, skewness 
and kurtosis parameters as well as selected probability levels based on the 
Noimal, Log- Normal and Veibull probability distribution functions. It also 
estimates a Level of Significance based on the Chi-Squared Goodness of Fit 
Test. 


The call to this subroutine is: 

CALL WBL (N, Y, NG, P) 

N - is the number of values in the sample population. 

Y - is the sample population array. 

NG - is the number of cells to be used for the CHI-Squared 

Goodness of Fit Test. 

P - is an array of 20 numbers containing the output 

information. 

P(l) - is sample mean 
f( 2) - standard deviation 

p(3) - skewness 

F(4) - kurtosis 

p(S) - mean of logarithm of values 
P(6) - standard deviation of logarithm of values 

P(?) - skewness of logarithm of values 

P(s) - kurtosis of logarithm of values 
F(9) - 0*995 probability level, 0*95 confidence level, 

for Normal Distribution 

P(lO) - predicted 0.995 probability level, 0*95 
confidence level, based on Log- 
Nozmal distribution 

F(11) - 'a' parameter for Weibull Distribution 
P(l2) - 'b' parameter for Weibull Distribution 
P(l5) - 'c' parameter for Weibull Distribution 
F(14) - 0*5 probability level for Weibull Distribution 
E(15) - 0*990 probability level for Weibull Distribution 

F(16) - 0*995 probability level for Weibull Distribution 

P(l7) - 0.999 probability level for Weibull Distribution 

P(l8) - level of significance for Normal Distribution 

P(l9) - level of significance for Log-Normal Distribution 
F(20) - level of significance for Weibull Distribution 
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3*4 Input Data Description 


Input data is in fixed field format and includes fixed point, hollerith 
and floating point fields. Figure 7 presents a sample run input listing. 

Basic groups of input include the following: 

1) one card of fixed point numbers including all of the run options, run 
number, number of random samples and number of boost phase integration 
steps, 

2) two cards of Hollerith data describing the run for output heading, 

3) initial random sequence integer for the pseudo random number generator, 

4) moment of inertia and center of mass table, 

5) list of single values of vehicle, mission and control system constants, 

6) list of statistical means and standard deviations of random variables 
for the boost phase, 

7) optional separate coast control variables, 

8) tables of vehicle, trajectory time histories and wind versus altitude 
tables, 

9) Optional vehicle and trajectory inputs replacing Group 8 


3.4.1 Group 1 Input (Input Options) 

The first card of input includes eleven fixed point constants in fields of 
five, right justified, starting in columns 1 through 5. These are: 


Card Column 

FORTRAN Variable 

Description 

1-5 

• NRUN 

An arbitrary run number 

6-10 

IT 

Number of random cases to be computed by 
the Monte Carlo procedure (current 
dimensions limit this to 1000 samples) 

15 

10(1) 

Output print option 


1 Gives output information for each sample 
and the statistical distribution 
information 

2 Gives only the statistical distribution 
information 

20 10(2) 2 Assumes the coast control and boost 

controls are the same 

3 Coast control system uses different 
input than boost system (this requires 
additional input under Group 7) 
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Card Coliunn FORTRAN Variable 


Deacription 


' 25 10(3) 1 

2 


30 

10(4) 

1 



2 

34-35 

10(5) 



39-40 

10(6) 


45 

10(7) 

0 



1 

50 

10(8) 

0 



1 

55 

10(9) 

0 


1 


Calculate boost fuel and coast time 
Calculate boost, coast and total fuel 
(If 10(2) * 3; retro time is computed 
based on input total fuel available in 
Group 5) 

Deletes control system pitch and yaw 
filter delay time in coast after time 
specified in Group 5 input 
Control system filter delay used in 
boost and coast 

Number of integration steps (up to 50 ) 
to be used in boost fuel consumption 
during web time (web time is specified 
in Group 5 data) 

Number of integration steps to be used 
during booster motor tail off between 
web time and burnout time specified in 
Group 5 data 

Separate coast control logic if 10(2) = 
3, 

Coast roll motors used for yaw and roll 
Coast pitch and yaw motors are the same 

Input all trajectory tables in Group 8 
Input Group 9 data instead of Group 8 

Print histograms of fuel consumption 
No histogram output 


3*4*2 Group 2 Title 


Group 2 includes two cards of Hollerith arbitrary data for the user to 
describe the run* It includes only the first seventy two (72) columns on each 
card* This information is written on the first page of output of the run* 


3 . 4,3 Group 3 Initial Random Number Seed 


The fourth card of input is an integer number for the start of the 
pseudo-random number generator* It is located in the first eighteen (I 8 ) 
columns on the card and is right justified* On IBM machines having seven 
digit accuracy this should be changed to a seven digit integer (eg* 123456?) 
and a Format change would be required* This also requires the optional coding 
in the subroutine RNDX shown in the FORTRAN listing* 
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3«4*4 Group 4 Mass Properties 


This group provides the transverse moment of inertia and center of mass 
station versus fraction of booster propellant consumed. It uses three points 
(O, 0.5 and 1 . 0 ) of propellant consumed to fit the data to a second order 
polynomial. The nine values are stored in array TABLE 1. Data is entered 
three to a card in three fields of 10 columns as follows. 


Col. 1-10 

11-20 

21-30 

values at 

0.0 

ly (Slug-Ft2) 

Xq 2 (Sta. in.) 

ignition 

0.5 

ly 

V 

^cg 

50^ 

1.0 

ly 

T 

•'^cg 

burnout 


3 . 4.5 Group 5 Single 

Constants 



This group used on all runs consists of single valued constants. There 
are twenty-one cards in this group. Each card has a single value located in 
the first fifteen columns (Format E15»5)» 13ie first sixteen values are stored 
in the "A" array. These start with card eight. 


FORTEAN 


Card No. 

Variable 

Symbol 

Units 

Description 

8 

A (1) 

N 


number of samples for 
correlating thrust 
misalignment initial value 
with slope 

9 

A (2) 



correlation coefficient for 
pitch thrust misalignment 

10 

A (3) 



correlation coefficient for 
yaw thrust misalignment 

11 

A (4) 

Xc 

inches 

body station of control 
motors 

12 

A (5) 

Zc 

inches 

radial location of pitch and 
yaw control motors 

13 

A (6) 

Xj 

inches 

body station of booster 
nozzle throat or action 
point of thrust misalignment 

14 

A (7) 

Rc 

inches 

roll control motor moment 
arm from centerline 

15 

A (8) 

cos \ 

— 

cosine of the pitch and yaw 
control motors forward cant 
angle 
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Card No. 

FORTRAN 

Variable 

Symbol 

Units 

Description 

16 

A (9) 

sin X 

•• 

sine of the pitch and yaw 
control motor forward cant 
angle 

17 

A (10) 

tf 

seconds 

time from ignition during 
which the filter delay is 
operating with the control 
system 

18 

A (11) 

■tcoaat 

seconds 

coast time from burnout 

19 

A (12) 


pounds 

total control fuel available 
(used to compute retro time 
or coast time with 
appropriate options) 

20 

A (15) 

'^'web 

seconds 

booster web time or action 
time from ignition 

21 

A (14) 

^bo 

seconds 

booster burnout time after 
ignition 

22 

A (15) 

Wbo 

pounds 

vehicle weight at booster 
burnout 

25 

A (16) 

Ix 

slug-ft2 

vehicle roll moment of 
inertia at burnout 

24 

CNAS 

Cn«s 

ft2/rad 

aerodynamic normal force 
coefficient slope 

25 

XCP 

Xcp 

inches 

vehicle aerodynamic center 
of pressure station 

26 

GAMEl 

>'ei 

degrees 

standard deviation of 
initial flight path angle 

27 

ZEI 

^ei 

degrees 

standard deviation of 
initial azimuth 

28 

QFRAC 

— 

— 

standard deviation of 
dynamic pressure ratio 
( Qac tual/ ^nominal ) 
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3 •4*6 Group 6 Bandom Variables 


This group includes the mean and standard deviation of related boost system 
random variables. These are read in two "Values per card (format 2E15.5)* 

This group starts with card 29 and continues through card 59. 


Card No. 

FORTRAN 

Variable 

Sjnnbol 

Units 

Description 

29 

0RM,0RS 

Fp/Fc 


boost pitch and yaw, control 
motor overshoot ratio mean 
and standard deviation for 
second order thrust response 
(Pig. 4) (if ORM 1.0, 
constant square wave thrust 
is assumed) 

30 

TIM, TPS 

tp 

seconds 

pitch and yaw control motor 
time to peak overshoot 
thrust mean and standard 
de"viation (Ignored if ORM = 
1) 

31 

T2IM,TRFS 


seconds 

effective turn off delay 
time due to vehicle 
flexibility, mean and 
standard deviation (not used 
when filter is included) 

32 

B(l), B(2) 

^^op 

degrees 

initial value of pitch 
thrust misalignment (mean 
and standard deviation) 

33 

B(3), B(4) 

^'’oy 

degrees 

initial value of yaw thrust 
misalignment (mean and 
standard deviation) 

34 

B(5), B(6) 


deg/sec 

mean and standard deviation 
of time rate of change of 
pitch thrust misalignment 

35 

B(2), B(8) 


deg/sec 

mean and standard deviation 
of time rate of change of 
yaw thrust misalignment 

36 

B(9), B(10) 

Isp 

seconds 

boost control motor specific 
impulse mean and standard 
deviation 

37 

B(ll), B(12) 

00 

deg/sec 

mean and standard deviation 
of initial pitch rate 

38 

B(13), B(14) 


degrees 

mean and standard de-viation 
of initial pitch attitude 
error 


- 26 - 



FORTEAN 


Card No. 

Variable i 

Symbol 

Units 

Description 

39 

B(15), B(16) 

k 

deg/sec 

mean and standard deviation 
of initial yaw rate 

40 

B(17), B(187 

K 

degrees 

mean and standard deviation 
of initial yaw attitude error 

41 

B(19), B(20) 

RTI 

ft-lb-sec 

booster induced roll angular 
impulse mean and standard 
deviation 

42 

B(21), B(22) 

^JET 

degrees 

mean and standard deviation 
of pitch and yaw control 
motor misalignment 
contributing to roll moment 

43 

B(25), B(24) 


degrees 

mean and standard deviation 
of wind direction azimuth 

44 

B(25), B(26) 

Aly 

Axcg 

slug-ft2 

inches 

standard deviation of trans- 
verse moment of inertia and 
standard deviation of center 
of mass station in inches 

45 

C(l), C(2) 

Atf 

seconds 

pitch and yaw filter delay 
time mean and standard 
deviation 

46 

C(3), 0(4) 

■^rg 

seconds 

mean and standard deviation 
of pitch and yaw gyro time 
delay 

47 

0(5), 0(6) 


seconds 

mean and standard deviation 
of pitch and yaw control 
motor tum-on delay time 

48 

0(7), 0(8) 

^2m 

seconds 

mean and standard deviation 
of pitch and yaw control 
motor turn-off delay time 

49 

0(9), 0(10) 

Kr/Kd 

seconds 

boost pitch and yaw control 
system gain ratio -njean and 
standard deviation 

50 

0(11), 0(12) 

db 

degrees 

boost pitch and yaw control 
deadband halfwidth mean and 
standard deviation 

51 

C(13), 0(14) 

H 


mean and standard deviation 
of boost pitch and yaw 
control hysteresis fraction 
(l“doff/don) 

52 

C(15), 0(16) 

Fc 

pounds 

boost pitch and yaw control 


motor thrust mean and 
standard deviation 
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FORTRAN 


Card No. 

Variable 

Symbol 

Units 

Description 

53 

C(17), C(18) 

^RG 

seconds 

roll rate gyro delay time 
mean and standard deviation 

54 

C(19), C(20) 

^IM 

seconds 

boost roll motor tum-on 
delay time mean and standard 
deviation 

55 

C(21), C(22) 

^2M 

seconds 

boost roll motor turn-off 
delay time mean and standard 
deviation 

56 

C(23), C(24) 

Kr/Kj) 

seconds 

boost roll control gain 
ratio mean and standard 
deviation 

57 

C(25), C(26) 


degrees 

boost roll deadband 
halfwidth mean and standard 
deviation 

58 

C(27), C(28) 

H 


boost roll switching 
hysteresis (l - 
<^off/don) 

standard deviation 

59 

C(29), C(30) 

Fc 

pounds 

boost roll control motor 

J.1 - _ j- .j _ J. j j 


thrust mean and standard 
deAn.ation 


3.4*7 Group 7 Separate Coast System Variables (I0(2) = 3 Only) 


This group is an optional input if a separately defined coast control 
system is used. This infoimation must be input if 10 (2) = 3 on the first 
card. This allows separate control motors, gains, time delays and a retro 
function during coast using the canted forward boost control motors. The 
first four cards are input one value per card with format KL5.5. They are; 


Card No. 

FORTRAN 

Variable 

Symbol 

Units 

Description 

60 

PTRC 

Aflc 

deg/sec 

total pitch program rate 
changes during coast 

61 

YTRC 


deg/sec 

total yaw program rate 
change during coast 

62 

TRETRO 

^ retro 

seconds 

retro time required 
following coast with four 
boost pitch and yaw control 
motors 
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FORTRAN 

Card No» Variable Symbol Units Description 

63 RETRSI Igp seconds retro motor specific impulse 

In addition to these four cards, there are nine ( 9 ) additional cards of 
statistical control system values having two values per card (format 2E15»5)» 
These are: 


Card No. 

FORTRAN 

Variable 

Symbol 

Units 

Description 

64 

D(l), 

D(2) 

dc 

degrees 

coast system pitch and yaw 
deadband halfwidth, mean and 
standard deviation 

65 

D(3), 

D(4) 

Fc 

pounds 

coast pitch control motor 
force, mean and standard 
deviation 

66 

D(5), 

D(6) 

^=1M 

seconds 

coast pitch motor tum-on 
delay time, mean and 
standard deviation 

67 

D(7), 

D(8) 

^2M 

seconds 

coast pitch motor turn-off 
delay time, mean and 
standard deviation 

68 

D(9), 

D(10) 

dc 

degrees 

coast roll deadband 
halfwidth, mean and standard 
deviation 

69 

D(ll) 

, D(12) 

Fc 

pounds 

coast roll motor thrust, 
mean and standard deviation 

70 

D(13) 

, D(14) 

^IM 

seconds 

coast roll motor tum-on 
delay time, mean and 
standard deviation 

71 

D(15) 

, D(16) 

■^2M 

seconds 

coast roll motor turn-off 
delay time, mean and 
standard deviation 

72 

D(17) 

CO 

1—1 

^sp 

seconds 

coast control specific 
impulse, mean and standard 
deviation 


3»4»8 Group 8 Trajectory Tables IO(S) = 0 


This group concludes the input data stream for a run. It includes tables 
of time and altitude histories. If 10 (8 ) = 1 certain tables are deleted as 
described in Group 9 m ibe following paragraph. Each table is input in the 
same format. Precediiig each table a card having the integer number of values 
in the table is located in columns 1 through 5» right justified (Format 15). 
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The values are entered alternating abscissas and ordinates In fields of ten 
columns. These are entered six (6) per card until the table is completed. 
Figure 7 is the input for the sample problem which shows this format. The 
format for the whole table is (I5/(6P10.3)) * A description of the tables and 
their order is: 


TABLE NAME 

ORDINATE 

UNITS 

ABSCISSA 

UNITS 

DESCRIPTION 

TABLE 1 

Q 

lbs/ft2 

t 

sec 

dynamic pressure versus 
time after ignition 

TABLE 2 

T 

lbs 

t 

sec 

booster nominal thrust 
versus time after 
ignition 

TABLE 3 

Wp 

lbs 

t 

sec 

booster propellant 
weight remaining versus 
time after ignition 

TABLE 4 

^c 

deg 

t 

sec 

pitch program attitude 
from local horizontal 
versus time after 
ignition 

TABLE 5 

V 

ft/sec 

t 

sec 

nominal velocity versus 
time after ignition 

TABLE 6 

’'pR 

degrees 

t 

sec 

flight path angle from 
local horizontal versus 
time after ignition 

TABLE 7 

^R 

degrees 

t 

sec 

nominal azimuth versus 
time after ignition 

TABLE 8 

^\oM 

degrees 

t 

sec 

nominal angle of 
sideslip versus time 
after ignition on an 
undisturbed trajectory 

TABLE 9 

% 

ft/sec 

h 

ft 

mean value of wind 
velocity versus altitude 
above sea level 

TABLE 10 

% 

ft/sec 

h 

ft 

standard deviation of 
wind velocity versus 
altitude above sea level 

TABLE 11 

h 

feet 

t 

sec 

altitude versus time 
after ignition 


Each of the above tables is limited to ninety (90) values which is a 
maximum of forty-five (45) pairs of abscissas and ordinates. 
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3.4»9 Group 9 Trajectory Tables I0(8) = 1 


For parametric studies the option to generate nominal trajectory 
parameters becomes useful. This option generates a gravity turn trajectory 
from other input data. With this option Group 9 data replaces Group 8 data. 
Figure 8 presents a sample input for the case. The first card includes four 
variables read in four fields of ten columns (format 8F10.3)* These variables 
are: 

FOETHAN 


VAEIABLE 

SYMBOL 

UNITS 

DESCRIPTION 

QO 

Qo 

lbs/ft2 

nominal dynamic pressure at 
Ignition 

VO 

Vo 

ft/sec 

nominal velocity at ignition 

GMO 

^0 

degrees 

nominal flight path angle at 
stage ignition 

CDS 

CpS 

ft2 

aerodynamic drag coefficient 
times reference area 


Following this card seven tables are entered. Each table is preceded by 
a card containing the integer number of values in the table in a field of five 
( 5 ) right justified (format 15)* Following this card the table is entered in 
floating point numbers in fields of ten (lO) columns, six (6) values per 
card. The format is (6F10.3)* The order and description of these tables are 
as follows: 


TABLE NAME 

ORDINATE 

UNITS 

ABSCISSA 

UNITS 

DESCRIPTION 

Table 2 

T 

lbs 

t 

sec 

booster nominal thrust 
versus time after 
ignition 

Table 3 

W 

lbs 

t 

sec 

booster nominal thrust 
versus time after 
Ignition 

Table 4 

0 c~VpR 

degrees 

t 

sec 

difference between 
nominal pitch attitude 
and flight path angle 
versus time after 
Ignition 

Table 7 

*^R 

degrees 

t 

sec 

predicted azimuth time 
history 

Table 8 

^NOM 

degrees 

t 

sec 

nominal angle of 
sideslip versus time 
after ignition 

Table 9 

Vw^i 

ft/sec 

h 

feet 

mean value of wind, 
velocity versus altitude 

Table 10 

VWa 

ft/sec 
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h 

feet 

standard deviation of 
wind velocity versus 
altitude 


3.5 Output Description 


All output is on a line printer. There are several options for the 
output which are selected by the run options presented in paragraph 3*4. 1. 

There are five basic parts to the output, two of which can be deleted by 
ii^ut option. These include: 

1) printout of selected parameters from each random sample (10(1) = 1) , 
can be deleted if (I0(l) = 2), 

2) printout of all fuel consumption values arranged in ascending order 
with associated sample population probability level (always output), 

3) sample distribution statistical parameters, final random number 
generator seed and number of samples having deadband overshoot 
occurrances, 

4) summary statistics and level of significance for Normal, Log-Normal 
and Weibull Distribution functions, 

5) histograms of the fuel consumption (I0(9) = O), can be omitted if 
(10(9) = 1) 

Details of these parts of the output are presented in the following 
subparagra|di s. 


3 . 5.1 Individual Sample Parameters 


If the input option (10(1) = 1), two lines of parameters from each Monte 
Carlo random sample case are printed out after each case. This option results 
in a relatively large amount of printout. A typxcal page of output of these 
parameters is presented in Figure 9» A description of each in the order of 
output is: 

OUTPUT 

NAME 

UNITS 



DESCRIPTION 

SMPL 

— 



an integer sequence of the Monte Carlo 
case nunber 

CAPT IMP 

lb-sec 



staging "capture" transient impulse 

ROLL TOBQIMP 

lb-sec 



roll control motor impulse due to 
rocket booster induced roll torque 
(RTI) 

BOOSTIMP 

lb-sec 



total control impulse required during 
boost phase 

t(web) 

sec 



booster web time 

T(T0) 

sec 



booster tail off time 
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OUTPUT 

HAKE 

UNITS 

DESCRIPTION 

BOOST FUEL 

lbs. 

control fuel consumption during boost 

ISP 

sec 

control fuel specific impulse 

CST FUEL 

lbs 

control fuel consumption during coast 
(not including retro fuel) 

FLOW RATE 

Ibs/sec 

coast fuel consumption flow rate with 
boost system - filter in 

DB RED TORQ 
FUEL 

lbs 

fuel consumption required for deadband 
reduction transient and pitch program 
rate changes during coast 

FLOW RATE 
COAST 

Ibs/sec 

fuel flow rate during coast 

RETRO TIME 

sec 

retro time achieved with remaining 
fuel consumption after coast 

NO. OF DB OVERSHOOTS 
BOOST P-Y ROLL CST P- 

-Y 

number of times deadband overshoot 
occured during boost in pitch or yaw, 
in the roll channel, and during coast 
in pitch or yaw 

TOTAL FUEL 

lbs 

total fuel consumption including 
boost, coast and retro 


When the coast time option is used, (I 0 ( 3 ) = 1 ), (TOTAL FUEL) and (COAST 
fuel) is deleted and (COAST TIKE) is output in seconds. 


3.5.2 Discrete Probability Levels 


Fuel consumption in boost and coast, total fuel consumption and under 
some options retro time or coast time is printed versus probability of 
occurrence after being arranged in ascending order. A sample output is 
presented in Figure 10 . 


3*5.3 Sample Distribution and Final Seed 


This page of output Is presented in Figure 11. It includes the computed 
mean, standard deviation and estimated 99.5 percent probability with 95 
percent confidence level of fuel consumption. 

The last value calculated of the pseudo random number generator seed is 
output. This can be used as a starting value if an additional run is to be 
made to increase the sample size with the least risk of repeating parts of the 
sequence generated in the previous run. 
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Ube total number of cases having exceeded the deadband overshoot criteria 
at least once is also output on this page. Differentiation between boost and 
coast phases is made as well as roll axes during boost. The deadband 
overshoot criteria is an indicator of a very poorly designed on-off control 
system since it approaches a 100 percent duty cycle. 


3.3.4 Summary Statistics and Level of Significance 


A page of statistical information is output for comparison of the sample 
population with known distributions. Figure 12 presents a typical example. 
Sample population mean, standard deviation, skewness and kurtosis parameters 
are output under Normal Distribution. The estimated upper 95 percent 
confidence level at the 99*5 percent probability level is based on a 
prediction if the distribution were truly Gaussian (which is usually not the 
case) . Level of Significance is computed based on the Chi-Squared Goodness of 
Pit Test. The level of significance is output for each of the distributions. 

A higher level of significance generally indicates a better fit. However, 
histogram output should also be considered in the statistical analysis. 

Under the Log-Normal Distribution the sample distribution parameters for 
the natural logarithm of the fuel consumption is shown. The fuel consumption 
at the upper 95 percent confidence level at 99*5 percent probability level 
predicted by the Log-Nomal distribution is also shown in pounds. Note that 
the parameters are expressed as the logarithm but the 99.5 percent level is 
not the logarithm. Level of significance is also presented from the 
Chi-Squared Goodness of Fit Test. 

Weibull Distribution parameters (a, b, c) are presented which are based 
on the sample population. Values of fuel consumption at several selected 
probability levels as predicted by the fit to a Weibull Distribution are also 
output. Level of Significance from a Chi-Squared Goodness of Fit Test to the 
ideal Weibull Distribution having these parameters is also output. 


3.5.5 Histograms 


Histograms of the sample populations of fuel consumption and retro times 
are an optional output. Samples of this output are presented in Figure 13 . 
Each output histogram includes: 

. Identification 
. Total number of samples 
. Minimum and maximum values 
. Number of values in each cell printed 
numerically below each cell 
. Numerical value of starting and end points 
of each cell printed under each cell 
. Graphical display of each cell by number 
of stars or asterisks rising from the base 
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Maximum number of samples in a cell is 94 for display purposes. If a 
cell is completely filled the actual height can be obtained from the number 
printed at the base of the cell. The number of cells is automatically set at 
50 or more over the range of data. In order to select aesthetically pleasing 
cell divisions the number of cells will vary somewhat eind occassionally spill 
over more than a single page of output. 
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Figure 5 (Cont.) 

General Flow Chart of MAIN Routine 
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Figure 7 

Sample Problem Input 10(8) » 0 


1 499 2 3 

2 1 

12 4 1 0 0 


SAMPLE PROBLEM SCOUT 

SECOND STAGE, FUEL OPTION, FILTER OUT IN 


COAST, SEPARATE COAST 

CONTROL, 

NO VAU-ROLL MIXING IO(8)-0 


625699253066457 


INITIAL RANDOM SEQ NO RANF 


0. 49241.45 

282.07 

IVY AND XCG 0IGNITION (SLUG- 


0.5 41574.21 

258.47 

050X PROP FT**2) 


1.0 27991.96 

206.73 

0BURNOUT 


19. 


NO. OF SAMPLES FOR CORR COEFF 


-0.573 


CORR COEFF PITCH 


-0.756 


CORR COEFF VAU 


467.76 


CONTROL MOTOR STATION 


14.93 


CONTROL MOTOR LAT. LOC. (IN.) 


448.51 


BOOSTER NOZZLE THROAT STATION 


16.89 


ROLL CONTROL MOTOR ARM (IN.) 


1. 


COS OF CANT ANGLE (NO CANT) 


0 . 


SIN OF CANT ANGLE (NO CANT) 


42.32 


FILTER IN TIME (SEC) 


34.92 


COAST TIME (SEC) 


167. 


CONTROL FUEL (LB) 


35.86 


BOOSTER UEB TIME (SEC) 


39.32 


BOOSTER BURNOUT TIME (SEC) 


6968.26 


UEIGHT B/0 (LB) 


205.40 


IXX B/0 (SLUG-FT**2) 


30.32 


CNAS FT2/RAD 


134.2 


XCP STATION 


0.30 


FLIGHT PATH ERROR (DEG) 

SIGMA 

0.38 


AZIMUTH ERROR (DEG) 

SIGMA 

0. 12 


DVN. PRESS. RATIO(Q OUER QNOM ) 

SIGMA 

1.25 

0.080 

PITCH/YAU MOTOR OUERSHOOT RATIO 

MEAN, SIGH 

0.019 

0.0047 

PITCH/YAU MOTOR TIME TO PEAK (SEC) 

MEAN, SIGN 

0.012 

0.001 

T2 FLEX FILTER OUT COAST (SEC) 

MEAN,SIGM 

0.0816 

0.0271 

INITIAL ET PITCH (DEG) 

MEAN,SIGM 

0.0098 

0.0239 

INITIAL ET YAU (DEG) 

MEAN, SIGN 

-0.000193 

0.000753 

SLOPE OF ET PITCH (DEG/SEC) 

MEAN,SIGM 



FIGURE 7 (Continued) 
Sample Problem Input I0(8)=0 


000308 

0.000830 

SLOPE OF ET VAU < DEG/SEC) 

MEAN^SIGM 

185. 

0. 

BOOST CONTROL ISP (SEC) 

MEAN^SIGM 

-0.098 

0.386 

INITIAL PITCH RATE (DEG^SEC) 

MEAN, SIGN 

0.898 

0.340 

INITIAL PITCH ERROR (DEG) 

MEAN, SIGN 

0.085 

0.197 

INITIAL VAU RATE (DEG/SEC) 

MEAN, SIGN 

0.0105 

0.843 

INITIAL VAU ERROR (DEG) 

MEAN,SIGM 

0. 

0. 

BOOSTER ROLL IMPULSE (FT-LB-SEC) 

MEAN,SIGM 

0. 

0.3817 

PITCH/YAU JET MISALIGNMENT (DEG) 

MEAN,SIGM 

859.0 

80.0 

UIND AZIMUTH (DEG) 

MEAN,SIGM 

0.0166 

0.0033 

IVY AND XCG STANDARD DEUIATION 

SIGMA, SIGN 

0. 

5. 

PITCH STATIC UNBALANCE (FT-LB) 

MEAN,SIGM 

0. 

5. 

YAU STATIC UNBALANCE (FT-LB) 

MEAN, SIGN 

0.0‘459 

0.0034 

FILTER DELAY TIME (SEC) 

MEAN,SIGM 

0.0093 

0.0017 

PITCH /YAU GYRO DELAY TIME (SEC) 

MEAN, SIGN 

0.0715 

0.0088 

PITCH/YAU MOTOR T1 (SEC) 

MEAN,SIGM 

0.0355 

0.0058 

PITCH/YAU MOTOR T8 (SEC) 

MEAN,SIGM 

0.5 

0.01667 

KR/KD (SEC) 

MEAN, SIGN 

0.808 

0.0867 

BOOST PITCH/YAU DEADBANDS (DEG) 

MEAN,SIGM 

0.035 

0.0167 

BOOST PITCH/YAU HYSTERESIS (FRACTION )MEAN,SIGM 

517.3 

14.5 

BOOST PITCH YAU CONTROL FORCE (LB) 

MEAN,SIGM 

0.0069 

0.0008 

ROLL RATE GYRO DELAY (SEC) 

MEAN, SIGN 

0.0851 

0.0048 

BOOST ROLL MOTOR T1 (SEC) 

MEAN,SIGM 

0.0168 

0.0030 

BOOST ROLL MOTOR T8 (SEC) 

MEAN,SIGM 

0.45 

0.015 

ROLL GAIN RATIO KR/KD (SEC) 

MEAN, SIGN 

1 .438 

0.0478 

BOOST ROLL DEADBAND (DEG) 

MEAN, SIGN 

0.035 

0.0167 

ROLL HYSTERSIS (FRACTION) 

MEAN,SIGM 

45.45 

1.68 

BOOST ROLL MOTOR FORCE (LB) 

MEAN,SIGM 

0.45 


PITCH PROGRAM RATE CHANGE (DEG/SEC) 


0. 


YAU PROGRAM RATE CHANGE (DEG/SEC) 


0. 


RETRO TIME (SEC) 


134.0 


RETRO ISP (SEC)MUST BE NON-ZERO 


0.785 

0.0867 

PITCH YAU COAST DEADBANDS (DEG) 

MEAN,SIGM 

517.3 

14.5 

COAST PITCH MOTOR THRUST (LB) 

MEAN,SIGM 

0.0715 

0.0088 

COAST PITCH MOTOR T1 (SEC) 

MEAN,SIGM 



FIGURE 7 (Continued) 
Sample Problem Input 10(8) =0 



0.0355 

0.0053 COAST 

PITCH MOTOR T3 (SEC) 

MEAN,SIG 


1*433 

0.0478 COAST 

ROLL DEADBAND (DEG) 

MEAN«SIG 


45.45 

1.68 COAST 

ROLL MOTOR THRUST (LB) 

MEAN.SIG 


0.0351 

0.0048 ROLL 

MOTOR T1 

(SEC) 

MEAN^SIG 


0.0163 

0.0030 ROLL 

MOTOR T3 

(SEC) 

MEAN.SIG 


135. 

0 

. COAST 

ISP (SEC) 

MEAN.SIG 

30 

NOMINAL 

DYN. PRESS. 

(PSF) US. 

TIME (SEC) 


0. 

99.93 

3.84 

76.65 

8.84 

56.36 


13.84 

43.73 

18.84 

33.18 

33.84 

35.49 


38.84 

18.30 

33.84 

11.95 

38.84 

5.89 


43.84 

3.38 






34 

BOOSTER 

THRUST (LBS) 

US. TIME 

(SEC) 



0. 

40000. 

0.43 

40643. 

6.75 

51961.4 


13.53 

61736.3 

30.34 

71354.0 

35.54 

73568.9 


30.84 

70483.1 

34.5 

69967.6 

35.86 

63794.0 


36.77 

10389.4 

38.06 

1039.0 

39.33 

0. 


13 

BOOSTER 

PROPELLANT 

UEIGHT REMAINING (LBS) US. TIME (SEC) 


0. 

8374.34 

9.63 

6548.78 

30.34 

4077.35 


38.91 

1844.78 

35.65 

188*53 

39.33 

0. 


18 

NOMINAL 

PITCH ATTITUDE (DEG) 

US. TIME(SEC) 


0. 

30.58 

3.84 

39.73 

8.84 

38.31 


13.84 

36.97 

18.84 

35.73 

33.84 

34.50 


38.84 

33.60 

33.84 

33.86 

43.84 

31.34 


13 

UELOCITV 

' (FT/SEC) US. TIME (SEC) 



0. 

4394. 

8.84 

5183. 

18.84 

6666. 


38.84 

8849. 

38.84 

10897. 

43.84 

10839. 


18 

NOMINAL 

FLIGHT PATH 

ANGLE (DEG) US. TIME (SEC) 


0. 

31.05 

3.84 

39.75 

8.84 

38.18 


13.84 

36.79 

18.84 

35.59 

33.84 

34.53 


38.84 

33.63 

33.84 

33.86 

43.84 

31.58 



6 NOMINAL AZIMUTH(DEG) US. TIME(SEC) 

0. 183.79 18.84 183.843 43.84 

6 NOMINAL ANGLE OF SIDESLIP < DEG ) US. TIME 

0.16 18.84 0.30 43.84 


0 


183.90 

(SEC) 

0.33 



FIGURE 7 CConcluded) 
Sample Problem Input I0C8)^0 


8 

tlEAN UIND 

UELOCITV(FT/SEC ) 

US. 

ALTITUDE(FEET) 

130000. 

130. 

160000. 

360. 


330000. 

310. 

400000. 

310. 






8 

SIGflA UIND UELOCITY(FT/SEC) 

US 

. ALTITUDE(FEET) 

130000. 

75. 

160000. 

100. 


330000. 

100. 

400000. 

100. 






10 

ALTITUDE 

(FEET) US. 

TINE (SEC) 



0. 

134795. 

8.84 

145378. 


18.84 

171644 

38.84 

303371. 

43.84 

363430. 
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Figure 8 

Sample Problem Input 10(8) = 1 


1 799 a 3 

3 1 

13 4 1 1 0 

000001 

SAMPLE PROBLEM SCOUT 

SECOND STAGE, FUEL OPTION, FILTER OUT IN 

000003 

COAST, SEPARATE COAST 

CONTROL, 

NO VAU-ROLL MIXING 10(8 )-l 

000003 

374635699353066457 


INITIAL RANDOM SEQ NO RANF 

000004 

0. 49341.45 

383.07 

IVV AND XCG 0IGNITION (SLUG- 

000005 

0.5 41574.31 

358.47 

050?« PROP FT**3) 

000006 

1.0 37991.96 

306.73 

0BURNOUT 

000007 

19. 


NO. OF SAMPLES FOR CORR COEFF 

000008 

-0.573 


CORR COEFF PITCH 

000009 

-0.756 


CORR COEFF YAU 

000010 

467.76 


CONTROL MOTOR STATION 

000011 

14.93 


CONTROL MOTOR LAT. LOG. (IN.) 

000013 

448.51 


BOOSTER NOZZLE THROAT STATION 

000013 

16.89 


ROLL CONTROL MOTOR ARM (IN.) 

000014 

1. 


COS OF CANT ANGLE (NO CANT) 

000015 

0. 


SIN OF CANT ANGLE (NO CANT) 

000016 

43.33 


FILTER IN TIME (SEC) 

000017 

34.93 


COAST TIME (SEC) 

000018 

167. 


CONTROL FUEL (LB) 

000019 

35.86 


T UEB (SEC) 

000030 

39.33 


T BURNOUT (SEC) 

000031 

6968.36 


UEIGHT B/0 (LB) 

000033 

305.40 


IXX B/O (SLUG-FT:r*a) 

000033 

30.33 


CNAS FT3/RAD 

000034 

134.3 


XCP STATION 

000035 

0.30 


FLIGHT PATH ERROR (DEG) 

000036 

0.38 


AZIMUTH ERROR (DEG) 

000037 

0. 13 


Q FRACTION SIGMA 

000038 

1.35 

0.080 

PITCH/VAU MOTOR OUERSHOOT RATIO 

000039 

0.019 

0.0047 

PITCH/YAU MOTOR TIME TO PEAK (SEC) 

000030 

0.013 

0.001 

T3 FLEX FILTER OUT COAST (SEC) 

000031 

0.0816 

0.0371 

INITIAL ET PITCH (DEG) 

000033 

0.0098 

0.0339 

INITIAL ET YAU (DEG) 

000033 

-0.000193 

0.000753 

SLOPE OF ET PITCH (DEG/SEC) 

000034 



000308 

0.000830 

185» 

0. 

-0.098 

0.386 

0.898 

0.340 

0.085 

0.197 

0.0105 

0.843 

0. 

0. 

0. 

0.3817 

859.0 

80.0 

0.0166 

0.0033 

0. 

10. 

0. 

10. 

0.0459 

0.0034 

0.0093 

0.0017 

0.0715 

0.0088 

0.0355 

0.0058 

0.5 

0.01667 

0.808 

0.0867 

0.035 

0.0167 

517.3 

14.5 

0.0069 

0.0008 

0.0851 

0.0048 

0.0168 

0 . 0030 

0.45 

0.015 

1.438 

0.0478 

0.035 

0.0167 

45.45 

1.68 

0.45 


0. 


1.85 


165. 


0.785 

0.0867 

517.3 

14.5 

0.0715 

0.0088 


FIGURE 8 (Continued) 
Sample Problem Input I0(8)=l 


SLOPE OF ET YAU (DEG/SEC) 

BOOST CONTROL ISP (SEC) 

INITIAL PITCH RATE (DEG/SEC) 

INITIAL PITCH ERROR (DEG) 

INITIAL VAU RATE (DEG/SEC) 

INITIAL VAU ERROR (DEG) 

BOOSTER ROLL ANG» INPULSE (FT-LB-SEC) 
PITCH/YAU JET PlISALIGNtnENT (DEG) 

UIND AZIMUTH (DEG) 

lYY AND XCG STANDARD DEVIATION 

PITCH STATIC UNBALANCE (FT-LB) 

YAU STATIC UNBALANCE (FT-LB) 

FILTER DELAY TIME (SEC) 

PITCH /YAU GYRO DELAY TIME (SEC) 
PITCH/YAU MOTOR T1 (SEC) 

PITCH/YAU MOTOR T3 (SEC) 

KR/KD (SEC) 

BOOST PITCH/YAU DEADBANDS (DEG) 

BOOST PITCH/YAU HYSTERESIS (FRACTION) 
BOOST PITCH YAU CONTROL FORCE (LB) 
ROLL RATE GYRO DELAY (SEC) 

BOOST ROLL MOTOR T1 (SEC) 

BOOST ROLL MOTOR T3 (SEC) 

ROLL GAIN RATIO KR/KD (SEC) 

BOOST ROLL DEADBAND (DEG) 

ROLL HYSTERESIS (FRACTION) 

BOOST ROLL MOTOR FORCE (LB) 

PITCH PROGRAM RATE CHANGE (DEG/SEC) 
YAU PROGRAM RATE CHANGE (DEG/SEC) 
RETRO TIME (SEC) 

RETRO ISP (SEC) 

PITCH VAU COAST DEADBANDS (DEG) 

COAST PITCH MOTOR THRUST (LB) 

COAST PITCH MOTOR T1 (SEC) 


000035 

000036 

000037 

000038 

000039 

000040 

000041 

000048 

000043 

000044 

000045 

000046 

000047 

000048 

000049 

000050 

000051 

000058 

000053 

000054 

000055 

000056 

000057 

000058 

000059 

000060 

000061 

000068 

000063 

000064 

000065 

000066 

000067 

000068 



FIGUF£. 8 CConcluded) 
Sample Problem Input 10 CS):*! 



0.0355 

0.0052 COAST 

PITCH MOTOR T2 (SEC) 

000069 


1.432 

0.0478 COAST 

ROLL DEADBAND (DEG) 

000070 


45.45 

1.68 COAST 

ROLL MOTOR THRUST (LB) 

000071 


0.0251 

0.0048 ROLL 

MOTOR T1 

(SEC) 

000072 


0.0162 

0.0030 ROLL 

MOTOR T2 

(SEC) 

000073 


125. 

COAST 

ISP (SEC) 

000074 

100. 

4394. 

31.05 4.2 

INITIAL 

Q,U, GAMMA, CDS 

000075 

34 


BOOSTER THRUST TIME HISTORY 


000076 

0. 

40000 . 

0.43 40643. 

6.75 

51961.4 

000077 

13.53 

61736.2 

20.24 71254.0 

25.54 

73568.9 

000078 

30.84 

70482.1 

34.5 69967.6 

35.86 

63794.0 

000079 

36.77 

10289.4 

38.06 1029.0 

39.32 

0. 

000080 

12 

BOOSTER 

PROPELLANT WEIGHT REMAINING TIME HISTORY 

000081 

0. 

8274.34 

9.63 6548.78 

20.24 

4077.25 

000082 

28.91 

1844.78 

35.65 188.52 

39.32 

0. 

000083 

18 


THETA - GAMMA US. 

TIME 


000084 

0. 

-0.47 

3.84 -0.03 

8.84 

0.03 

000085 

13.84 

0.18 

18.84 0.14 

23.84 

-0.03 

000086 

28.84 

-0.02 

33.84 0.00 

43.34 

-0.34 

000087 

6 


AZIMUTH US. TIME 



000088 

0. 

182.79 

18.84 182.842 

43.84 

182.90 

000089 

6 


BETA NOMINAL US. TIME 


000090 

0. 

0.16 

18.84 0.20 

43.84 

0.23 

000091 

8 


MEAN UIND US. ALTITUDE 


000092 

120000. 

120. 

160000. 260. 

220000. 

310. 

000093 

400000. 

310. 




000094 

8 


SIGMA UIND US. ALTITUDE 


000095 

120000. 

75. 

160000. 100. 

220000. 

100. 

000096 


*EOR 



Hgure 9 

Sample Output - Individual Cases - 10 (l) 


1 


RUN NO. > PAGE NO. 1 

SANPLE PROtLEN SCOUT SECOND STAGE. FUEL OPTION. FILTER OUT IN 
COAST. SEPARATE COAST CONTROL. NO VAU-ROLL NIXING lOIl)*! 


RUN NO. 1 PAGE NO. 2 


SNPL 

CAPT INP 
(LB-SEC) 
CST FUEL 
(POUNDS 

ROLL TORO 
INP(LB-SC) 
FLOU RATE 
FILTER IN 

BOOST INP 
(LB-SEC) 

DB RED TORO 
FUEL(LBS) 

T(UEB) 
(SEC) 
FLOU RATE 
COAST 

T(TO> BOOST FUEL ISP 

(SEC) (LBS) (SEC) 

RETRO TINE NO. OF DB OUERSHOOTS 

(SEC) BOOST P-V ROLL CST P-V 

TOTAL FUEL 
(LBS) 

1 

5.4E4GET62 

1.8943ET61 

a. 

1.72G4E>ea 

S.179aETa3 

i.s8a4EToa 

3.58saETai 

3.9416E-at 

3.46aaETaa 

5.962aEiaa 

4.9432ETat 

a 

i.2saaE«a2 

a a 

6.9956E4ai 

2 

S.4886E>62 

Z.lSSaEtOl 

a. 

1 .4532E«aa 

4.9719ETa3 

i.ssasETaa 

S.SSSBETat 

5.a4S9E-ai 

3.46aaE4B0 
6 567eETaa 

3.9775ETai 

a 

i.25aaE«a2 

a a 

6.3336E*ai 

3 

S.568SET02 

1.5S82E»61 

a. 

i.4S37E*aa 

8.2154E+e3 

1.3367ET0a 

3.586eET01 

3.2a48E-ai 

3.46aaEToa 

5.2454ETaa 

6.6723ETei 

a 

i.25aaE4a2 

a a 

8.2642E*ai 

4 

S.3516ET62 

1.9809E»61 

a. 

1 .62i9ETae 

7.4289ET03 

1.5676ETaa 

3.s8saE4ai 

4.2794E-01 

3.46aaE«aa 

5.2437E4aa 

5.9431E4ai 

a 

i.25aeE4a2 

a a 

s.aaasETai 

5 

4.6822E*62 

1.4933ET61 

a. 

1.277EET0a 

4.5942E'»a3 

i.424iEtaa 

3.586eE*ai 

3.1786E-ai 

3.46aaETaa 

7.2418E«aa 

3.6753E^ai 

a 

i.aseaEtaa 

a a 

5.3110E*ai 

6 

4.7494ET02 

2.1593ET61 

a. 

1 .saesETaa 

6.7367E4S3 

1.5498ETSa 

3.58aoE4ai 

4.8943E-ai 

3.46aeE4aa 

5.4432E*aa 

S.3894ETei 

a 

i.2saaE«a2 

a a 

7.7a37E«ai 

7 

4.4536ET62 

2.2665ET61 

a. 

1 .54fi4E4ee 

7.9478ET63 

1.6S3SE«aa 

3.58saE4ai 

4.eoaie-ei 

3«4saaE48a 

s.aaiSEToa 

6.3582ETai 

a 

i.2saaETa2 

a a 

8.73aiETai 

8 

5.8624ET62 

1.7617ET61 

a. 

i.3i97ETaa 

7.304SE«a3 
1 .5743ETaa 

3.5868ETai 

3.9111E-ai 

3.4saaE4aa 

5.6S74ETaa 

5.8436ETai 

a 

i.25aaE4a2 

a a 

7.7627ETai 

9 

4.1946ET62 

1.3521E»ai 

a. 

1.4828E*ae 

G.SSSaETOS 
1 .4i3iEtee 

3.S860E^ai 

2.5981E-ai 

3.46aeE«aa 

5.9898E*aa 

5.3264Etei 

a 

i.2seaE4a2 

a a 

E.8198Eiai 

16 

4.7358E«e2 

1.7a26E«ei 

a. 

1.2577E*aa 

5.437aET03 
1 .449ac*aa 

3.58saE«ai 

3.7951E-ai 

3.4saaEtaa 

G.47SiETaa 

4.3496ET61 

a 

i.asaBEtaE 
a a 

6.1071ETai 

11 

S.4aS2E«a2 

2.a3asETai 

a. 

i.56t7ETee 

4.7S88ETa3 

1.5992E*aa 

s.sssaETOi 

4.4735E-ai 

3.4saaE*aa 

6.sa78ETaa 

3.8a7iETai 

a 

i.2saaE*a2 

a a 

5.9976Etai 

12 

3.4429ET62 

2.1519E«ai 

a. 

t.i722ETee 

$.2625E*e3 
1 .6ia2£Taa 

3.S8saETai 

5.1S53E-ai 

3.4saaETaa 

6.4662ETaa 

4.2iaaETai 

a 

i.EsaaETaa 

a a 

6.5229Ef01 

13 

6.a959E«a2 
1 .9S98E*ei 

a. 

1 .89eeE«ee 

s.ie23E>a3 
1 .5984ETaa 

3.586aETei 

3.988SE-at 

3.46aeE4^aa 

4.9586Etaa 

6.4819E«at 

a 

i.BseBETaB 

a a 

B.saiSEtai 

14 

6.2931E«a2 

2.26asETat 

a. 

1 .S378E»ea 

5.a2S8ETa3 
1 .6S42E>ea 

3.68seETai 

S.1534E-ai 

3.4saaETaa 

6.2493ETaa 

4.a2i4ETai 

a 

1 .25aaE*B2 

a a 

6.4478E4ai 

15 

4.5577E«a2 

2.1728Etai 

a. 

i.sssaEToa 

8.6983ETe3 

i.627eE«ae 

3.586aEtei 

5.a236E-at 

3.46eBEtaa 

4.5541E*a0 

6.9586ETai 

a 

i.asaaETBa 

a a 

9.2941E*ai 



Figure 10 

Sample Output - Discrete Probability 


RUN NO. 1 

SANPLC PROILEN SCOUT 

SECOND STAGE 

, FUEL OPTION, 

FILTER OUT IN 

PAGE N0< 

COAST, SEPARATE COAST 

CONTROL, NO 

VAU-ROLL NIXING 

I0(8)*e 


RUN NO. 1 

PROtACIUTV 

BOOST FUEL 

COAST FUEL 

RETRO TINE 

PAGE NO 
TOTAL FUEL 

(POUNDS) 

(POUNDS) 

(SECONDS) 

(POUNDS ) 

.aozoo 

22.729 

11.891 

7.8529 

39.521 

.••4«e 

22.936 

12.501 

7.7257 

43.684 

.eoeeo 

24.861 

12.974 

7.7145 

43.999 

.eose* 

25.540 

13.063 

7.6578 

44.486 


26.821 

13.068 

7.6064 

46.934 

.01200 

27.694 

13.135 

7.5769 

47.270 

.01400 

28.001 

13.521 

7.4796 

47.389 

.01600 

28.113 

13.585 

7.4234 

48.100 

.01800 

29.158 

13.699 

7.3787 

49.333 

.02000 

29.747 

13.790 

7.3498 

49.605 

.02200 

29.908 

13.816 

7.3384 

49.613 

.02400 

29.977 

13.831 

7.3149 

49.775 

.02600 

30.128 

13.908 

7.3149 

50.301 

.02800 

30.509 

14.081 

7.2684 

50.362 

.03000 

30.728 

14.214 

7.2443 

50.501 

.03200 

30.789 

14.269 

7.2418 

50.506 

.03400 

31.044 

14.380 

7.2233 

50.588 

.03600 

31.075 

14.438 

7.2182 

50.635 

.03800 

31.220 

14.688 

7.1969 

50.878 

.04000 

31.401 

14.770 

7.1951 

50.890 

.04200 

31.560 

14.787 

7.1945 

51.082 

.04400 

31.592 

14.788 

7.1917 

52.219 

.04600 

32.336 

14.814 

7.1916 

52.263 

.04800 

32.589 

14.889 

7.1857 

52.338 

.05000 

32.621 

14.894 

7.1787 

52.451 

.05200 

32.631 

14.933 

7.1600 

52.632 

.05400 

32.892 

14.943 

7.1310 

52.684 

.05600 

32.958 

14.985 

7.1260 

52.758 

.05800 

32.985 

15.018 

7.1168 

52.878 

.06000 

33.041 

15.029 

7.1086 

52.889 

.06200 

33.046 

15.042 

7.1039 

53.110 

.06400 

33.318 

15.184 

7.0836 

53.213 

.06600 

33.516 

15.25S 

7.0688 

53.277 

.06000 

33.S75 

15.288 

7.0657 

53.459 

.07000 

33.615 

15.332 

7.0648 

53.602 

.07200 

33.705 

15.369 

7.0565 

53.786 

.07400 

33.994 

15.397 

7.0413 

53.943 

.07600 

34.049 

15.416 

7.0124 

54.014 

.07000 

34 125 

15.429 

7.0119 

54.033 

.08000 

34.126 

15.445 

6.9678 

54.088 

.08200 

34.196 

15.469 

6.9607 

54.268 

.08400 

34.445 

15.558 

6.9453 

54.348 

.08600 

34.635 

15.582 

6.9426 

54.569 

.08800 

34.791 

15.583 

6.9342 

54.574 

.09000 

34.815 

15.602 

6.9309 

S4.S96 

.09200 

34.907 

15.631 

6.9213 

54.607 

.09400 

34.927 

15.668 

6.9161 

54.615 

.09600 

35.373 

15.S74 

6.8917 

54.845 


52 



FIGURE 10 CContinued) 

Sample Output - Discrete Probability 


RUN NO. t 


PR09A8ILITV 

BOOST FUEL 
(POUNDS) 

COAST FUEL 
(POUNDS) 

RETRO TINE 
(SECONDS) 

.81488 

39.987 

17.418 

6.5874 

.81688 

39.964 

17.432 

6.5818 

.81888 

39.994 

17. 464 

6.5740 

.88888 

48.863 

17.468 

6. 5678 

.88888 

48.185 

17.485 

6.5666 

.88488 

48.814 

17.494 

6.5665 

.88688 

48.814 

17.588 

6.5686 

.88888 

48.258 

17.514 

6.5589 

.83888 

48.869 

17.553 

6.5488 

.83888 

48.351 

17.554 

6.5394 

.83488 

48.476 

17.567 

6.5388 

.83688 

48.578 

17.573 

6.5342 

.83888 

48.629 

17.598 

6.5338 

.84888 

48.637 

17.688 

6.5862 

.84888 

48.654 

17.688 

6.5288 

.84488 

48.838 

17.615 

6.5168 

.84688 

48.845 

17.617 

6.5158 

.84888 

48.853 

17.618 

6.5145 

.85888 

48.948 

17.681 

6.5135 

.85888 

41.879 

17.658 

6.5183 

.85488 

41.175 

17.686 

6.5887 

.85688 

41.440 

17.788 

6.5888 

.85888 

41.518 

17.783 

6.5876 

.86888 

41.697 

17.786 

6.5868 

.86888 

41.798 

17.785 

6.5851 

.86488 

41.883 

17.787 

6.5881 

.86688 

41 .885 

17.788 

6.4869 

.86888 

41.838 

17.789 

6.4861 

.87888 

41.871 

17.746 

6.4786 

.87888 

41.934 

17.748 

6.4778 

.87488 

41.988 

17.754 

6.4751 

.87688 

48.888 

17.754 

6.4789 

.87888 

42.188 

17.755 

6.4674 

.88888 

42.114 

17.785 

6.4658 

.88888 

42.226 

17.836 

6.4498 

.88488 

48.832 

17.855 

6.4458 

.88688 

42.385 

17.987 

6.4433 

.88888 

42.312 

17.912 

6.4404 

.89888 

42.342 

17.966 

6.4383 

.89288 

48.488 

17.988 

6.4272 

.29488 

48.429 

18.884 

6.4268 

.29688 

42.476 

18.828 

6.4244 

.29888 

48.477 

18.839 

6.4138 

.38888 

48.528 

18.847 

6.4188 

.38888 

42.578 

18.862 

6.4017 

.38488 

42.681 

18.873 

6.3978 

.38688 

42.687 

18.189 

6.3946 

.38888 

48.775 

18.138 

6.3988 

.31888 

48.841 

18.144 

6.3819 

.31888 

42.873 

18.154 

6.3777 

.31488 

43.836 

18.154 

6.3767 

.31680 

43.858 

18.167 

6.3755 

.31888 

43.126 

18. 169 

6.3753 


- 53 - 


PAGE NO. 4 

TOTAL FUEL 
(POUNDS) 

si.asA 
61.314 
61 .388 
61.446 
61.449 
61.453 
61.479 
61.537 
61.663 
61.676 
61.886 

61.838 
61.887 
61.896 
61.964 
61.985 
61.954 
61.971 
68.816 

68.838 

68.313 

68.314 
68.489 
68.781 
68.744 

68.846 
68.989 
68.969 
68.976 

63.846 
63.867 
63.138 
63.151 
63.816 

63.838 
63.336 
63.386 
63.485 
63.598 
63.684 
63.613 
63.645 
63.787 

63.854 

63.855 
63.948 
64.885 
64.117 
64.136 
64.188 
64.188 

64.844 

64.858 



FIGURE. 10 (Continued) 

Sample Output - Discrete Probability 


RUN NO. 


PACE NO. 


PROtAltLlTV 


.38eeo 

.3SE0* 

.38400 

.33600 

.33800 

.33000 

.33300 

.33400 

.33600 

.33800 

.34000 

.34300 

.34400 

.34600 

.34800 

.35000 

.35300 

.35400 

.35600 

.35800 

.36000 

.36300 

.36400 

.36600 

.36800 

.37000 

.37300 

.37400 

.37600 

.37800 

.38000 

.38300 

.38400 

.38600 

.38800 

.39000 

.39300 

.39400 

.39600 

.39800 

.40000 

.40300 

.40400 

.40600 

.40800 

.41000 

.41300 

.41400 

.41600 

.41800 

.43000 

.43300 

.43400 


lOOST rUEL 
(POUNDS) 

43.171 

43.388 

43.431 

43.496 

43.895 

43.917 

43.984 

44.305 

44.311 

44.331 

44.394 

44.330 

44.343 

44.371 

44.439 

44.473 

44.638 

44.679 

44.715 

44.869 

44.938 
44.963 
44.969 
45.046 
45.319 
45.355 
45.386 
45.307 
45.447 
45.485 
45.568 
45.645 
45.776 
45.811 

45.870 

45.937 

45.938 
46.014 
46.110 

46.314 

46.315 

46.383 
46.537 
46.559 
46.561 
46.657 
46 669 
46.733 
46.806 

46.833 
47.033 
47.040 
47.049 


COAST FUEL 
(POUNDS) 

18.184 
18.187 
18. 194 
18.335 
18.331 

18.364 
18.308 
18.337 
18.339 
18.335 
18.355 

18.364 

18.383 

18.384 
18.388 
18.390 
18.405 

18.407 

18.408 

18.408 
18.433 
18.437 
18.453 
18.466 
18.489 
18.493 
18.504 
18.511 
18.535 
18.549 
18.578 
18.585 
18.608 
18.615 
18.636 
18.665 
18.681 
18.688 
18.690 
18.698 
18.733 
18.739 
18.753 
18.766 
18.773 
18.775 
18.785 
18.790 
18.796 
18.803 
18.837 
18.848 
18.856 


RETRO TINE 
(SEC0ND5) 

6.3741 

6.3741 

6.3654 

6.3603 

6.3567 

6.3550 

6.3451 

6.3367 

6.3316 

6.3310 

6.3398 

6.3361 

6.3198 

6.3173 

6.3138 

6.3113 

6.3039 

6.3011 

6.3877 

6.3817 

6.3786 

6.3757 

6.3754 

6.3754 

6.3746 

6.3733 

6.3583 

6.3549 

6.3493 

6.3483 

6.3460 

6.3459 

6.3446 

6.3398 

6.3376 

6.3378 

6.3331 

6.3110 

6.3088 

6.3056 

6.3047 

6.3030 

6.3018 

6.1997 

6.1935 

6.1848 

6.1736 

6.1670 

6.1568 

6.1558 

6.1537 

6.1519 

6.1459 


TOTAL FUEL 
(POUNDS) 

64.396 

64.431 

64.478 

64.633 

64.706 

64.779 

64.779 

64.867 

64.868 
64.879 
64.908 
64.944 
64.973 
65.013 
65.095 
65.116 
65.118 
65.339 
65.359 
65.383 
65.537 
65.583 
65.588 
65.656 
65.744 
65.778 
65.781 
65.805 

65.833 
65.841 

65.939 
65.967 
66.379 
66.590 
66.637 
66.737 
66.773 

66.833 
66.875 
66.887 
66.934 
67.114 
67.149 

67.365 
67.549 
67.553 
67.576 
67.616 
67.851 
67.856 
67.914 
67.916 
68.108 



FIGURE 10 CContlnued) 

Sample Output - Discrete Probability 


RUN NO. 


PACE NO. 6 


PROlAtILlTV 

BOOST ruEL 
(POUNDS > 

COAST FUEL 
(POUNDS) 

RETRO TIRE 
(SECONDS) 

TOTAL FUEL 
(POUNDS) 

.4E6RR 

47.695 

18.896 

6. 1443 

68.198 

.43BM 

47.185 

18.966 

6.1467 

68.368 


47.266 

18.943 

6.1390 

68.313 

.oaaa 

47.272 

18.945 

6.1386 

£8.357 

.434e* 

47.331 

18.946 

6.1323 

68.449 

.43£0» 

47.522 

18.956 

6.1366 

£8.491 

.4380* 

47.525 

18.984 

6.1287 

68.580 

.44eo« 

47.534 

18.986 

£.1262 

£8.636 

.4420* 

47.545 

18.987 

£.1247 

68.646 

.44468 

47.654 

19.618 

6.1187 

£8.849 

.44668 

47.774 

19.619 

6.1134 

68.993 

.44868 

47.831 

19.049 

6.1110 

£9.606 

.45668 

47.916 

19.052 

6.1051 

£9.195 

.45266 

47.958 

19.075 

6.1617 

69.242 

.45466 

48.606 

19.084 

6.6946 

69.356 

.45666 

48.034 

19.691 

6.0963 

69.384 

.45866 

48.235 

19.699 

6.0836 

69.442 

.46666 

48.266 

19.162 

6.6754 

69.566 

.46266 

48.321 

19.269 

6.6613 

£9.531 

.46466 

48.329 

19.234 

6. 6591 

69.569 

.46666 

48.337 

19.248 

6.6567 

69.721 

.46866 

48.435 

19.274 

6.6472 

69.746 

.47666 

48.446 

19.287 

6.6455 

69.863 

.47266 

48.541 

19.292 

6.6444 

69.867 

.47466 

48.596 

19.296 

6.6438 

69.819 

.47666 

48.716 

19.366 

6.6465 

69.956 

.47866 

48.765 

19.316 

6.6283 

69.969 

.48666 

48.769 

19.321 

6.6285 

69.978 

.48266 

48.836 

19.328 

6.6268 

76.662 

.48466 

48.835 

19.355 

6.6256 

76.618 

.48666 

48.835 

19.355 

6.6243 

76.683 

.48866 

48.876 

19.433 

6.6231 

76.144 

.49666 

43.872 

19.454 

6.6267 

76.266 

.49266 

48.984 

19.458 

6.6134 

76.263 

.49466 

49.116 

19.463 

6.6167 

76.363 

.49666 

49.168 

19.476 

6.6677 

76.369 

.49866 

49.172 

19.482 

6.6655 

76.331 

.56666 

49.188 

19.560 

5.9998 

76.342 

.56266 

49.261 

19.521 

5.9969 

76.379 

.56466 

49.432 

19.555 

5.9936 

76.422 

.50606 

49.435 

19.558 

5.9898 

76.453 

.50866 

49.492 

19.564 

5.9888 

76.547 

.51066 

49.566 

19.583 

5.9788 

76.599 

.51266 

49.534 

19.589 

5.9738 

76.643 

.51466 

49.585 

19.598 

5.9626 

76.656 

.51666 

49.601 

19.660 

5.961? 

76.836 

.51866 

49.676 

19.602 

5.9662 

76.856 

.52666 

49.753 

19.667 

5.9569 

76.926 

.52266 

49.796 

19.621 

5.9568 

76.835 

.52460 

49.963 

19.631 

S.9S56 

71.642 

.52606 

50.660 

19.659 

5.9538 

71.127 

.52860 

56.610 

19.668 

5.9483 

71.148 

.53066 

56.094 

19.688 

5.9446 

71.161 
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FIGURE. IQ CGotitlnued) 

Sample Output - Discrete Probability 


RUH NO. 1 


PROtAIILITV 

800ST FUEL 

COAST FUEL 

RETRO TINE 

(POUNDS) 

(POUNDS) 

(SECONDS ) 

.53200 

50.112 

19.779 

5.9369 

.53400 

50.178 

19.787 

5.9324 

.53600 

50.178 

19.789 

5.9304 

.53800 

50.215 

19.809 

5.9277 

.54000 

50.221 

19.819 

5.9274 

.54200 

50.234 

19.827 

5.9217 

.54400 

50.349 

19.832 

5.9212 

.54600 

50.424 

19.835 

5.9181 

.54800 

50.494 

19.839 

5.9160 

.55000 

50.494 

19.870 

5.9067 

.55200 

50.601 

19.880 

5.9060 

.55400 

50.791 

19.916 

5.8960 

.55600 

50.965 

19.948 

5.8929 

.55800 

51.005 

19.956 

5.8915 

.56000 

51.117 

19.974 

5.8840 

.56200 

51.182 

19.987 

5.8829 

.56400 

51.263 

20.004 

5.8825 

.56600 

51.297 

20.020 

5.8770 

.56800 

51.316 

20.032 

5.8722 

.57000 

51.392 

20.057 

5.8647 

.57200 

51.569 

20.093 

5.8637 

.57400 

51.574 

20.122 

5.8626 

.57600 

51.593 

20.129 

5.8583 

.57800 

51.622 

20.130 

5.8541 

.58000 

51 .671 

20. 167 

5.8528 

.58200 

51.749 

20.173 

5.8508 

.58400 

51.858 

20.182 

5.8500 

.58600 

52.034 

20.202 

5.8487 

.58800 

52.036 

20.210 

5.8464 

.59000 

52.052 

20.227 

5.8451 

.59200 

52.146 

20.269 

5.8444 

.59400 

52.219 

20.295 

5.8345 

.59600 

52.260 

20.296 

5.8296 

.59800 

52.270 

20.296 

5.8291 

.60000 

52.299 

20.306 

5.8202 

.60200 

52.358 

20.323 

5.8199 

.60400 

52.364 

20.341 

5.8024 

.60600 

52.449 

20.367 

5.7955 

.60800 

52.580 

20.376 

5.7951 

.61000 

52.606 

20.421 

5.7934 

.61200 

52.609 

20.435 

5.7917 

.61400 

52.652 

20.442 

5.7913 

.61600 

52.705 

20.470 

5.7909 

.61800 

52.711 

20.470 

5.7898 

.62000 

52.748 

20.514 

5.7882 

.62200 

52.754 

20.518 

5.7789 

.62400 

52.760 

20.567 

5.7740 

.62600 

52.764 

20.576 

5.7736 

.62800 

53.038 

20.577 

5.7727 

.63000 

53.048 

20.584 

5. 7558 

.63200 

53.077 

20.591 

5.7548 

.63400 

53.090 

20.606 

5.7527 

.63600 

53.090 

20.617 

5.7440 


PACE NO. 

TOTAL FUEL 
(POUNDS) 

71.219 

71.220 
71.257 
71 .355 
71.357 
71.387 
71.451 
71.463 
71.559 
71.641 
71.672 
71.770 
71.945 
71.982 
71.994 
72.079 
72.273 

72.298 

72.299 
72.312 
72.428 
72.453 
72.518 
72.549 
72.613 
72.757 
72.759 
72.783 
72.793 

72.796 

72.797 
72.890 
72.956 
73.001 
73.027 

73.221 
73.359 
73.362 
73.367 
73.437 
73.603 
73.658 
73.693 
73.713 
73.723 

73.797 
73.808 
73.911 
73.916 
73.916 
73.966 
73.980 
74.136 



figure 10 CContinued) 

Sample Output — Blsctete Probability 


RUN NO. 


PACE NO. 8 


PROIAKIITV 

BOOST FUEL 

COAST FUEL 

RETRO TINE 

TOTAL FUEL 


(POUNDS ) 

(POUNDS) 

(SECONDS) 

(POUNDS) 

.C3888 

53.897 

28.626 

5.7361 

74.184 

.64«6t 

53.189 

28.653 

5.7294 

74.389 

.64800 

53.193 

28.656 

5.7271 

74.446 

.6440R 

53.241 

28.692 

5.7243 

74.514 

.S4Se« 

53.264 

28.738 

5.7239 

74.521 

.648«« 

53.284 

28.749 

5.7237 

74.548 

.65«e« 

53.363 

28.745 

5.7235 

74.576 

.6SS«« 

53.367 

28.777 

5.7175 

74.593 

.6S40O 

53.433 

28.733 

5.7178 

74.855 

.65680 

53.479 

28.786 

5.7167 

74.978 

.65888 

53.498 

28.798 

5.7130 

74.985 

.66888 

53.654 

28.798 

5.7855 

74.986 

.66888 

53.783 

28.826 

5.7824 

75.333 

.66488 

53.712 

28.831 

5.6995 

75.345 

.66688 

53.737 

28.831 

5.6959 

75.583 

.66888 

53.847 

28.835 

5.6846 

75.667 

.67888 

53.894 

28.856 

5.6665 

75.667 

.67388 

53.952 

28.893 

5.6574 

75.683 

.67488 

53.957 

28.982 

5.6522 

75.782 

.67688 

54.287 

28.923 

5.6488 

75.786 

.67888 

54.249 

28.925 

5.6453 

75.798 

.68888 

54.258 

28.926 

5.6299 

75.944 

.68288 

54.268 

28.958 

5.6232 

76.181 

.68488 

54.385 

28.952 

5.6283 

76.199 

.68688 

54.678 

28.958 

5.6198 

76.259 

.68888 

54.688 

28.962 

5.6888 

76.262 

.68888 

54.712 

28.996 

5.6875 

76.296 

.69288 

54.724 

28.996 

5.6868 

76.587 

.69488 

54.778 

21.884 

5.6834 

76.513 

.69688 

55.158 

21.818 

5.6012 

76.528 

.69888 

55.368 

21.858 

5.6811 

76.569 

.78888 

55.581 

21.863 

5.6888 

76.659 

.78288 

55.543 

21.866 

5.5975 

76.856 

.78488 

55.575 

21.899 

5.5908 

77.837 

.78688 

55.589 

21.181 

5.5865 

77.858 

.78888 

55.729 

21.183 

5.5662 

77.139 

.71888 

55.782 

21.113 

5.5648 

77.147 

.71288 

55.791 

21.191 

5.5612 

77.187 

.71488 

55.838 

21.287 

5.5575 

77.277 

.71688 

56.878 

21.216 

5.5558 

77.397 

.71888 

56.189 

21.229 

5.5491 

77.428 

.72888 

56.274 

21.231 

5.5472 

77.576 

.72288 

56.355 

21.254 

5.5467 

77.627 

.72488 

56.449 

21.387 

5.5398 

77.786 

.72688 

56.484 

21.326 

5.5162 

77.797 

.72888 

56.585 

21.334 

5.5137 

77.799 

.73088 

56.655 

21.365 

5.5117 

77.982 

.73288 

56.987 

21.366 

5.4968 

77.947 

.73480 

56.834 

21.371 

5.4824 

78.182 

.73688 

57.841 

21.371 

5.4886 

78.447 

.73888 

57.149 

21.448 

5.4665 

78.452 

.74888 

67.181 

21.468 

5.4653 

78.578 

.74288 

57.382 

21.467 

5.4591 

78.511 
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FIGURE. 10 (Continued) 
Sample Output - Discrete Probability 


RUN NO. I 


PROBAtlLlTY 

BOOST FUEL 
(POUNDS ) 

COAST FUEL 
(POUNDS ) 

RETRO TiriE 
(SECONDS) 

.85999 

62.665 

22.699 

5.1219 

.85299 

63.013 

22.659 

5,1071 

.85499 

63. 188 

22.654 

5.0959 

.85699 

63.270 

22.663 

5.0945 

.85899 

63.337 

22.692 

5.0834 

.86999 

63.506 

22.891 

5.0735 

.86299 

63.582 

22.892 

5.0699 

.86499 

63.624 

22.951 

5.0655 

.86699 

63.799 

22.982 

5.0541 

.86899 

63.813 

23.914 

5.9486 

.87099 

63.826 

23.970 

5.0468 

.37299 

63.930 

23.977 

5.9453 

.87499 

63.937 

23.993 

5.9415 

.87690 

64.992 

23.176 

5. 0198 

.87809 

64.117 

23.291 

5.0158 

.88909 

64.159 

23.322 

5,9147 

.88899 

64.245 

23.331 

5.0141 

.88409 

64.344 

23.413 

5.0107 

.88699 

64.380 

23.428 

5.0969 

.88899 

64.668 

23.474 

5.0052 

.89099 

64.722 

23.482 

4.9907 

.89899 

64.784 

23.521 

4.9894 

.89499 

64.819 

23.577 

4.9651 

.89699 

64.869 

23.641 

4.9637 

.89899 

64.872 

23.791 

4.9586 

.99999 

64.875 

23.891 

4.9474 

.99299 

65.905 

23.815 

4.9413 

.99499 

65.969 

23.831 

4.9369 

.99699 

65.277 

23.862 

4.9348 

.99899 

65.661 

23.874 

4.9287 

.91999 

65.723 

23.877 

4.9129 

.91899 

65.825 

23.894 

4.9122 

.91499 

65.828 

23.916 

4.8993 

.91609 

66.388 

23.936 

4.8838 

.91899 

66.459 

24.961 

4.3807 

.98090 

66.622 

24.101 

4.8519 

.92209 

66.632 

24.128 

4.3263 

.92490 

67.939 

24.165 

4.8256 

.92699 

67.970 

24.168 

4.3201 

.92890 

67.523 

24.258 

4.8171 

.93000 

67.777 

24.644 

4.7769 

.93290 

67.874 

24.648 

4.7689 

.93409 

68.122 

24.660 

4.7554 

.93609 

68.917 

24.708 

4.7S33 

.93800 

69.023 

24.759 

4,7531 

.94000 

69.147 

24.775 

4.7311 

.94200 

69.174 

24.898 

4.7290 

.94400 

69.290 

24,829 

4,7185 

.94600 

69.449 

24.890 

4.7921 

.94809 

69.573 

24.891 

4.6766 

.95009 

69.536 

24.943 

4.6613 

.95200 

79.021 

24.943 

4.6002 

.95409 

70.136 

24.959 

4.5862 


PAGE NO. 

TOTAL PUEL 
(POUNDS) 

84.Z3E 

34.413 

84.691 

84.678 

84.767 

84.984 

85.199 

85.199 

35.198 

85.871 

85.417 

85.481 

85.594 
85.568 
85.699 
85.751 
85.799 
86.915 
86.866 
86.336 

86.846 
86.867 

87.188 

87. 189 
87.849 

87.391 
87.318 
87.683 
87.966 
87.997 
38.938 
88.115 
88.296 
88.879 
38.493 
83.424 

88.595 
88.575 
88.689 
88.695 
88.942 
89.188 
89.894 
99.211 
99.698 
99.716 

99.846 
99.969 
91.955 

91.392 
91.549 
91.644 
91.754 



FIGURE 10 CConcluded) 

Sample Output - Discrete Probability 


RUN NO. 1 PACE NO. It 


PROINOIUTV 

BOOST FUEL 
(POUNDS ) 

COAST FUEL 
(POUNDS) 

RETSO TINE 
(SECONDS ) 

TOTAL FUEL 
(POUNDS) 

.95690 

70.169 

25.029 

4.5736 

92.110 

.95890 

70.372 

25.107 

4.5541 

92.238 

.96990 

70.547 

25.268 

4.5511 

92.941 

.96299 

7t .4t2 

25.312 

4.S382 

93.462 

.96490 

71.468 

25.530 

4.4798 

93.953 

.96690 

71.545 

25.579 

4.4349 

94,409 

.96890 

7t.700 

25.681 

4.3756 

95.466 

.97990 

72.882 

25.694 

4.3737 

95.693 

.97290 

73.257 

25.695 

4.3549 

96.082 

.97460 

73.411 

25.698 

4.3376 

96.688 

.97690 

73.4X7 

25.798 

4.3165 

96.732 

.97800 

73.968 

26.164 

4.2579 

97,255 

.98090 

74.278 

26.615 

4.2546 

97.645 

.98290 

74.386 

26.666 

4.2486 

99.132 

.98400 

76.124 

26.726 

4.1131 

99.869 

.98600 

77.220 

26.766 

4.1107 

199.913 

.98800 

77.593 

26.874 

3.9945 

191.108 

.99000 

77.897 

28.161 

3.9288 

192.675 

.99200 

78.030 

28.779 

3.9255 

102.800 

.99400 

78.164 

28.882 

3.7219 

106.619 

.99600 

85.715 

30.642 

3.5038 

109.539 

.99800 

88.123 

30.970 

3.3252 

113.315 
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RUM MO. 1 


Figure 11 

Output Page of Sample Statistical lata 


PAGE MO. 


SAMPLE PROBLEM SCOUT SECOMD STAGE, FUEL OPTIOM, FILTER OUT IM 
COAST. SEPARATE COAST COMTROL. MO YAU-ROLL MIXING 10(8 )»0 


BOOST FUEL 

MEAN 49.783 

• STD. DEUIATION 11.490 

cr> 

° 0.995 PROB 88.179 

0.95 CONFID 


COAST FUEL 

RETRO TIME 

TOTAL FUEL 

19.660 

5.9573 

70.988 

3.039 

.7685 

13.133 

31.395 

3.3356 

113.183 


FINAL UALUE OF RANDOM SEQUENCE INTEGER - 374469945168434565 

NO. OF SAMPLES UITH DEADBAND OUERSHOOT 
BOOST(PITCH AND YAU ) - 0 

(ROLL) - 0 


COAST 


0 



Figure 12 

Summary Statistics and Level of Significance 


DISTRIBUTION FUNCTION FIT 

BOOST COAST TOTAL 

NORMAL DISTRIBUTION 


MEAN 

49.783 

19.660 

70.988 

STANDARD DEU 

11.478 

3.026 

12.120 

SKEUNESS 

.302 

.363 

.342 

KURTOSIS 

2.744 

3.408 

2.987 

0.995 PROB 0.95 CONF 

81.378 

27.990 

104.350 

LEUEL OF SIGNIFICNCE 

.035 

.044 

.014 

LOG NORMAL DISTRIBUTION 
(UALUES OF LOG) 

MEAN 

3.880 

2.967 

4.248 

STANDARD DEU 

.236 

.154 

.172 

SKEUNESS 

-.275 

-.134 

-.125 

KURTOSIS 

2.792 

3.142 

2.858 

0.995 PROB 0.95 CONF 

92.777 

29.722 

112.184 

LEUEL OF SIGNIFICNCE 

.107 

.756 

.516 

UEIBULL DISTRIBUTION 

PARAMETERS A 

21.6359 

12.5937 

42.2023 

B 

31.6198 

7.9491 

32.3661 

C 

2.6618 

2.5267 

2.5726 

PROBABILITY 0.5000 

49.201 

19.473 

70.284 

0 . 9900 

77.734 

27.136 

100.777 

0.9950 

80.801 

27.974 

104.092 

0.9990 

86.965 

29.667 

110.776 

LEUEL OF SIGNIFICANCE 

.637 

.003 

.293 



RETRO Tine HISTOGRAH 


Figure 13 
Histogram Output 
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BOOST FUEL CONSunPTION HISTOCRAPI 


Figure 13 (Cont.) 
Histogram Output 
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COAST FUEL CONSUMPTION HISTOGRAM 


Figure 13 (Cont.) 
Histogram Output 
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TOTAL FUEL CONSUnPTION HISTOCRAN 


Fi^re 13 (Concluded) 
Histrogram Output ^ 
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APPENDIX A 


FORTHAN Program Listing 


A complete listing of the FORTRAN Source Program is presented in the 
following pages. It starts with the MAIN routine and is followed by the 
seventeen (17) subroutines arranged in alphabetical order. There are a total 
of 947 cards in the MAIN routine. The total program including the subroutines 
contains 1,626 cards. 

If the routine is used on a non-CDC computer the pseudo random number 
generation will change. In this case the calls to RANSET (line 29) and RANGET 
(lines 668, 743, and 829) should be deleted. Subroutine RNDX has an alternate 
pseudo random number generator coding which is listed as comment cards. 


A-1 
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*DECK MAIN 1 

PROGRAd HAIN< INPUT, OUTPUT, TAPE5-INPUT,TAPE6-0UTPUT ) Z 

3 

SCOUT UPPER STAGE STATISTICAL FUEL CONSUMPTION ROUTINE 4 

COMPUTES THE STATISTICAL DISTRIBUTION OF CONTROL FUEL CONSUMPTION 5 
FOR AN ON-OFF REACTION JET ATTITUDE CONTROL SYSTEM BY THE MONTE 6 
CARLO TECHNIQUE. THE EFFECT OF BOOSTER INDUCED DISTURBANCES AND 7 
AERODYNAMIC DISTURBANCES ARE CONSIDERED IN ADDITION TO THE SELF 8 
EXCITED LIMIT CYCLE MOTION. 9 

DIMENSION 10(9), A(16), TABLK9), 8(30), C(30), D(18), BT(IS), 10 
1 PC(8), YC(8), RC(7), PDC(10), UB00ST( 1000 ) , TCOAST( 1000 ) , 11 

a UCOAS(1000),TRET(1000),UTOT(1000),FUB(20),FUC(a0),FUT(a0), 12 

3 TBL1(90),TBL2(90),TBL3(90),TBL4(90),TBL5(90),TBL6(90), 13 

> 4 TBL7(90),TBL8(90),TBL9(90),TBL10(90),TBL11(90), 14 

5 TA1(50),TCX(50),UXX(50),GAP(50),ZPX(50),PSX(50),HX(50), 15 

6 UELX(50),X1KA(50),X1KB(50),UUMX(50),UUSX(50),THX(50),TOU(50), 16 

7UGHT(50),XINER<50),CNASQ(50) 17 

COMMON K, RXD(30) 18 

COMMON/DOUT/NQT, TBLl , NUT, TBL5, NGT, TBL6, NHT, TBLl 1 , TBO, NTC, TBL4 19 

C0MM0N/DDI/TBL2, TBL3, NTHT, NUT, CDS , M2, M3, UBO 20 

DATA NUMCHI/11/ 21 

C READ RUN NO., NUMBER OF SAMPLES, AND RUN OPTIONS 22 

10 READ( 5,1330) NRUN, IT, ( I0( I) , 1=1,9 ) 23 

IF (EOF(5).NE.0) STOP 24 

C READ TUO CARDS OF TITLE INFORMATION FOR OUTPUT ONLY 25 

READ( 5,1350) 26 

C READ INITIAL SEQUENCE INTEGER FOR RANF PSEUDO RANDOM NUMBER 27 

READ( 5,1320) K 28 

CALL RANSET (K) 29 

NPAGE=1 30 

NLINE=0 31 

CALL PAGEHD (NRUN,NPAGE,NLINE ) 32 

URITE( 6,1350) 33 

NLINE-NLINE+2 34 
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C READ IN XCG AND MOMENT OF INERTIA TABLE 35 

READ( 5,1440) ( TABLi ( I) , I - 1 , 9 ) 36 

C READ IN SINGLE UALUED CONSTANTS 37 

READ( 5,1450) < A< I ) , I = 1 , 16 ) , CNAS, XCP, GAMEI , ZEI , QFRAC 38 

TBO»A(14) 39 

UBO-AC15) 40 

GAMEI-GAMEI/57.3 41 

ZEI-ZEI/57.3 42 

READ IN MEAN AND STANDARD DEVIATIONS OF CONTROL MOTOR OVERSHOOT 43 
RATIO, TIME TO PEAK THRUST, ADDED TURN-OFF TIME DELAY DUE TO 44 

STRUCTURAL FLEXIBILITY 45 

READ( 5,1460) ORM, ORS, TPM, TPS, T2FM, T2FS 46 

READ IN MEAN AND STANDARD DEVIATIONS OF BOOST AND CONTROL SYSTEM 47 
READ( 5,1460) (B( I ), I»l,30 ) 48 

READ( 5,1460) ( C ( I ) , I - 1 , 30 ) 49 

IF (IO(2)-3 .LT. 0 ) GO TO 20 50 

SEPARATE COAST CONTROL SYSTEM OPTION 10(2 )»3 51 

READ IN TORQUEING RATE CHANGES AND RETRO PARAMETERS 52 

READ( 5,1450) PTRC, YTRC, TRETRO, RETRSI 53 

C READ IN COAST CONTROL SYSTEM VARIABLES, MEAN AND SIGMA 54 

READ( 5,1460) (D( I ), I-l, 18 ) 55 

C READ IN TABLES 56 

20 M2=l 57 

M3=l 58 

C TEST FOR INPUT OPTION TO CALCULATE OR READ IN Q, VEL, GAMMA, ALTITUDE59 

IF (10(8). EG. 0) GO TO 30 60 

C READ IN INITIAL Q,V, GAMMA AND DRAG COEFFICIENT CDS FOR CALCULATIONS! 

READ( 5,990) QO, VO, GAMO, CDS 62 

READ( 5,1470) NTHT, (TBL2< I ), I-1,NTHT ) 63 

READ( 5,1470) NUT, < TBL3< I ) , I - 1 , NUT ) 64 

READ( 5,1470) NTC, ( TBL4 ( I ) , I - 1 ,NTC ) 65 

READ( 5,1470) NAZT, ( TBL7( I ) , I - 1 , NAZT ) 66 

READ( 5,1470) NSYT, ( TBL8 ( I ) , I - 1 ,NSYT ) 67 

READ( 5,1470) NVMH, ( TBL9( I ) , I - 1 , NVMH ) 68 
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READ( 5,1470) NUSH, ( TBL10 ( I ) , I - 1 , NUSH ) 69 

C CALCULATE Q,UELOCITY,GAni>lA, AND ALTITUDE TABLES 70 

CALL TABGEN (GO, UO, GANG) 71 

GO TO 40 78 

C 10(8 )-0 READ IN ALL TRAJECTORY TABLES 73 

30 READ< 5,1470) NOT, ( TBLl ( I ) , I -1 ,NQT ) 74 

READ( 5,1470) NTHT, ( TBL8 < I ) , I- 1 , NTHT ) 75 

READ( 5,1470) NUT, ( TBL3 ( I ) , I -1 ,NUT ) 76 

READ( 5,1470) NTC, ( TBL4 ( I ) , I -1 ,NTC ) 77 

READ( 5,1470) NUT, ( TBL5< I ) , I -1 , NUT ) 78 

READ( 5,1470) NGT, (TBL6( I ), 1 = 1 ,NGT ) 79 

READ( 5,1470) NAZT, ( TBL7 ( I ) , I- 1 , NAZT ) 80 

READ( 5,1470) NSYT, (TBL8( I ), I- 1,NSYT ) 81 

READ( 5,1470) NUflH, ( TBL9( I ) , I- 1 , NUPIH ) 82 

READ( 5,1470) NUSH, ( TBL10( I ) , I -1 ,NUSH ) 83 

READ( 5,1470) NHT, ( TBLl 1 ( I ) , I- 1 , NHT ) 84 

C COMPUTE CONSTANTS FOR SECOND ORDER CURUE FITTING MASS PROPERTIES 85 

40 QR1=TABL1(7 )/TABLl(4) 86 

QR2=TABL1(7)*(TABL1(7)-TABL1(4) ) 87 

SL2IY= (TABLl (8 )-TABLl (2 )+QRl*(TABLl (2 )-TABLl (5 ) ) )/QR2 88 

SLIIY- (-TABL1 (2 )+TABLl (5 )-SL2IY)t:TABLl (4 )*TABL1 (4 ) )/TABLl (4 ) 89 

SL2CG=(TABL1(9)-TABL1(3)+QR1*(TABL1(3)-TABL1(6) ) )/QR2 90 

SL1CG= ( -TABLl ( 3 )+TABLl ( 6 )-SL2CGXTABLl (4 ):icTABLl ( 4 ) )/TABLl ( 4 ) 91 

NTB=I0(5) 92 

NTT=I0(6) 93 

MER1=0 94 

MRER1=0 95 

MERC-0 96 

NPKl-0 97 

NPK2-0 98 

NPK3-0 99 

NPK4-0 1®® 

SUB-0. 1®1 

STC-0. 1®2 
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SUC-0. 

STR-0. 

SUT-0. 

ANTPU-IO(S) 1®6 

ANTPT-I0<6) 1®7 

TUEB-A<13) 108 

UTC-TBL3(E) 109 

DELT1=TUEB/ANTPU 110 

DELT2=(TB0-TUEB)/ANTPT HI 

PBFCT-0.0 112 

PCFCT=0.0 113 

YCFCT=0«0 114 

YBFCT-0.0 115 

C SET TABLE LOOKUP INDICES FOR INITIAL UALUE 116 

f11 = l 117 

na=i 118 

n 3 =i 119 

t14 = l 120 

f15 = l 121 

f16 = l 128 

n7=l 123 

f18»l 124 

ng=i 185 

M10=l 126 

mi = l 127 

NTOT=NTB+NTT 128 

DELTA=DELT1 129 

130 

CALCULATE AND STORE THRUST PUSALIGNMENT AND AERO COEFFICIENTS 131 

AT TIME POINTS DURING BOOST PHASE 138 

TA-0. 133 

DO 50 J-l,NTOT 134 

IF(J.GT.NTB) DELTA-DELT2 135 

TA-TA+0.5*DELTA 136 
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TA1(J)-TA 137 

CALL TBLU ( NOT, Q, TA. TBLl , fll ) 138 

CALL TBLU ( NTHT, T, TA, TBLa, na ) 139 

CALL TBLU ( NUT, UR, TA, TBL3, f13 ) 14© 

CALL TBLU (NTC, TC, TA, TBL4, f14 ) 141 

CALL TBLU (NUT,U,TA,TBL5,f15 ) 143 

CALL TBLU ( NOT, GAflP, TA, TBL6, MS ) 143 

CALL TBLU ( NAZT, ZP, TA, TBL7, M7 ) 144 

CALL TBLU ( NSYT, BETPR, TA, TBL8, (18 ) 145 

CALL TBLU (NHT, H,TA, TBLU, Nil ) 146 

CALL TBLU ( NUflH, UUN, H, TBL9, (19 ) 147 

CALL TBLU (NUSH,UUS,H,TBL10,(110 ) 148 

TCX(J)=TC/57.3 149 

UXX(J)-U 150 

GAP(J)-GA(1P/57.3 151 

ZPX(J)-ZP/57,3 15 a 

PSX(J)-BETPR/57.3 153 

UUnX(J)-UU(1 154 

UUSX(J)»UUS 155 

HX(J)«H 156 

PERCA=(UTC-UR)/UTC 157 

UTa-UBO+UR 158 

ACC“T*3a. 173/UTa 159 

UELX(J) = ACC*DELTA)CC(11 )/57.3 160 

C CALCULATE MOMENT OF INERTIA AND XCG AUERAGE DURING UEB BURN 161 

ERTIA=(TABLl(a)+SLlIV*PERCA+SLaiY*PERCA*PERCA) 163 

CGA- (TABLl ( 3 )+SLlCG^:PERCA+SLaCG:»:PERCA*PERCA ) 163 

XERT-ia.*ERTIA 164 

THX( J )-< (A(4 )-CGA ):»:A(8 )+A(5 )1CA(9 ) )/XERT 165 

X1KA( J)-CNAS*Q*(CGA-XCP)/XERT 166 

X1KB(J)-(A(6)-CGA):»:T/XERT 167 

XINER( J)-ERTIA 168 

T0U(J)-T/UT3 169 

UGHT(J)-UT3 170 
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CNASQ(J)-CNAS*Q 171 

50 CONTINUE 172 

173 

BEGIN THE MONTE CARLO RUNS 174 

175 

DO 750 ITER-1, IT 176 

C CALCULATE AND TEST CONTROL MOTOR THRUST OUERSHOOT FACTORS 177 

IF (ORM-1.0 ,LE. 0 ) GO TO 140 178 

CALL RNDX (8) 179 

ORPB-ORM+RXD ( 1 )*ORS 180 

ORYB-ORM+RXD(E)*ORS 181 

ORPC-ORM+RXD ( 3 ))|cORS 182 

ORVC-ORM+RXD ( 4 ):*:ORS 183 

TPPB-TPM+RXD ( 5 )*TPS 184 

TPYB-TPM+RXD ( 6 )*TPS 185 

5 TPPC-TPM+RXD(7):icTPS 186 

TPVC=TPM+RXD(8)*TPS 187 

IF (ORPB-1. .LE. 0 ) GO TO 70 188 

CALL OURF (TPPB,ORPB,PBFCT) 189 

GO TO 80 190 

70 PBFCT=0,0 191 

80 IF (ORPC-1.0 .LE. 0 ) GO TO 90 192 

CALL OURF (TPPC,ORPC,PCFCT) 193 

GO TO 100 194 

90 PCFCT-0.0 195 

100 IF (ORYB-1.0 .LE. 0 ) GO TO 110 196 

CALL OURF (TPVB,ORYB,YBFCT) 197 

GO TO 1E0 198 

110 YBFCT-0.0 199 

130 IF (ORYC-1.0 .LE. 0 ) GO TO 130 200 

CALL OURF (TPYC,ORYC,YCFCT) 201 

GO TO 140 203 

130 YCFCT-0.0 203 

140 CALL RNDX (1) 204 
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CALCULATE ADDITIONAL CONTROL MOTOR TURN-OFF DELAY DUE TO EOS 

STRUCTURAL FLEXIBILITY 206 

TaFLX-T2FM+RXD( 1 )*T2FS 207 

COMPUTE BOOST PHASE INITIAL CONDITIONS, ET, DISTURBANCES 208 

CALL RNDX (16) 209 

DO 150 J-1,12 210 

150 BT(J)-B(2*J-1)+RXD< J ):tcB(2*J) 211 

BT(13)-B(25)*RXD(13)+1. 212 

BT(14)«B(26))»cRXD(14) + l. 213 

SUP=B(27)+RXD(15)*B(28) 214 

SUY»B(29)+RXD(16)*B(30) 215 

CALCULATE UIND DIRECTION 216 

DU=BT(12)/57.3 217 

BT(3)“B(5)+B(6)*(RXD(1 )*A(2 )+RXD(3 )»SQRT( ( A( 1 )-l. )*( 1 .-A(2 )*A(2 ) )/218 
l(A(l)-2.))) 219 

BT(4)*B(7)+B(8):»c(RXD(2)*A(3)+RXD(4 )»SQRT( (Ad )-l. )«( 1 .-A(3 )*A(3 ) )/220 
l(A(l)-2,))) 221 

CALL RNDX ( 1 ) 222 

QRD-1 .+RXD( 1 )*QFRAC 223 

SELECT CONTROL SYSTEM RANDOM UARIABLES 224 

CALL RNDX (16) 225 

DO 160 J-1,8 226 

PC(J)-C(2*J-1)+RXD( J )*C(2*J) 227 

160 YC(J)-C(2*J-1 )+RXD(J+8)*C(2*J) 228 

CALL RNDX (7) 229 

DO 170 J=l,7 230 

170 RC(J) = C(2*J+15)+RXD(J)*C(2:»cJ + 16) 231 

TEST FOR STAGE OPTION 2-NO COAST CONTROL 3-COAST CONTROL 232 

IF (I0(2)-3 .LT. 0 ) GO TO 190 233 

CALCULATE RANDOM CONTROL SYSTEM CHARACTERISTICS FOR COAST UITH 234 

SEPARATE COAST SYSTEM USING YAU-ROLL MIXING 235 

CALL RNDX (10) 236 

DO 180 J-1,4 237 

PDC(J)-D(2)t:J-l )+RXD(J))HD(2)t:J) 238 
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180 PDC( J+4)-DCa*J+7)+RXD( J+4):*cD(a*J+8 ) 

PDC(9)-D(1 )+RXD(9)3tcD(a) 

PDC(10)-D(17)+RXD(10)*D(18) 

IF (I0<7)-1 .LT. 0 ) GO TO 190 

COMPUTE COAST CONTROL UARIABLES WITHOUT VAU-ROLL MIXING 
CALL RNDX (3) 

YCF1-D(3)+RXD( 1 )»D<4 ) 

VCTl-D(5)+RXD(a)*D<6) 

YCTa-D(7)+RXD(3)*D(8) 

COMPUTE ROLL TORQUES IMPULSE, CONTROL MOTOR TIME DELAYS, CONTROL 
ACCELERATIONS AT IGNITION 
190 RTI-ABS(BT(10) )*ia./A<7) 

TlP»PC(l)+PC(a)+PC(3) 

Tap»pc(i )+pc(a)+pc(4) 

T1Y=YC(1 )+YC(a)+YC(3) 

TaY«YC(l )+YC(a)+YC(4) 


839 

840 

841 
843 

343 

344 

345 

846 

847 

348 

349 

850 

851 
853 

853 

854 

THEDDA-PC(8)^( <A(4)-TABL1(3))ICBT(14) )*A(8)+A(9)«A(5) )/(18.*TABLl(3)355 
1«BT(13)) 356 

PSDDA-THEDDA*YC(8)/PC(8) 357 

COMPUTE CAPTURE MANUEUER IMPULSE 858 

CAPIMP = PC ( 8 ) Jtc ( 3 . CT3P+ABS ( ( BT ( 6 )+3 . * < BT ( 7 )-PC C 6 ) )/PC ( 5 ) )/ ( 57 . 3*THEDa59 
IDA) ) )+YC(8)*(3.)|:TaY+ABS( <BT(8)+3.!|:(BT(9)-YC(6) )/YC(5) )/(57.3:»:PSDDAa60 
3))) 361 

BIMP-0. 363 

THETP-0.0 863 

THETY=0.0 364 

CALL RNDX (3) 365 

UEL1-GAMEI1cTBL5(3)<:RXD(1) 366 

UEL3-ZEI*TBL5(3)*RXD(3) 367 

UINDEX»RXD(3) 368 

NERl-O 369 

NRRl-0 370 

NERC-0 371 

378 
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CALCULATE IMPULSE IN PITCH AND YAU DURING BOOST PHASE E73 

274 

DELTA-DELTl 275 

DO 300 J-l,NTOT 27S 

IF(J.GT.NTB) DELTA-DELT2 277 

TA-TAKJ) 278 

H-HX(J) 279 

TC»TCX(J) 280 

U=UXX(J) 281 

GAMP=GAP(J) 282 

ZP-ZPX(J) 283 

BETPR-PSX(J) 284 

CTHX=THX( J)/BT(13) 285 

UU=UUnX(J )+UINDEX*UUSX( J) 286 

SGAMP-SINCGAMP) 287 

o CGAMP-COS(GAMP) 288 

UUCSDZ=UU*COS(DU-ZP) 289 

XNUM1=UUCSDZ*SGAMP 290 

XNUM2“UU*SIN(DU-ZP) 291 

XDEN1“U+UUCSDZ*CGAMP 292 

QUDFAC = QRD*(XDEN1*2./U+(UU/U)*:t:2-1. ) 293 

XKA=QUDFAC*XlKA(J)/BT(i3) 294 

xkb=xik:b(J)/bt(13) 295 

C COMPUTE THRUST MISALIGNMENT 296 

ETPA*(BT( 1 )+BT(3)*TA)/57.3 297 

ETYA=(BT(2)+BT(4))»JTA)/57.3 298 

THEDDA»PC(8)*CTHX 299 

PSDDA = THEDDA)t:YC(8)/PC(8) 300 

C COMPUTE ANGLES OF ATTACK IN PITCH AND VAU 301 

GAMEl-UELl/U 302 

ZE1«=UEL2/U 303 

C CALCULATE ALPHA DUE TO WINDS 304 

ALPHU-CXNUMl/XDENl ) 305 

C CALCULATE BETA DUE TO WINDS 306 
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BETAU-(XNUria/XDENl) 307 

ALPHP-TC-GAMP-GAHEl+ALPHU 308 

BETAP-BETPR+BETAU+ZEl 309 

C CALCULATE PITCH AND VAU ACCEL DUE TO THRUST MISALIGNUENT , AERO 310 

C AND CG OFF-SET 311 

TOUI=TOU(J)/(XINER( J)»BT(13) ) 313 

THETP=XKA*ALPHP+XKB*ETPA+SUP3ICT0UI 313 

THETY--XKA*BETAP+XKB*ETYA+SUV*TOUI 314 

C PICK DEADBAND SIDE AND ADD DEADBAND TO ANGLE OF ATTACK 315 

IF (THETP) 310,330,300 316 

C COMPUTE CROSS TRACK UELOCITY ERROR UPDATE 317 

300 UEL1-UEL1+UELX(J) 318 

GAME-UELl/U 319 

DALP = C(11 )/57.3-(GAME-GAPIEl )/3. 330 

DADBP-DALP^:XKA 331 

THETP-THETP+DADBP 333 

GO TO 330 333 

310 UEL1»UEL1-UELX(J) 334 

GAME-UELl/U 335 

DALP--C(11 )/57.3-(GAnE-GAf1El )/3. 336 

DADBP-DALP^tcXKA 337 

THETP=THETP+DADBP 338 

330 IF (THETV) 340,350,330 339 

330 UEL3«UEL3+UELX( J ) 330 

ZE=UEL3/U 331 

DBET--C(11 )/57.3+(ZE-ZEl)/3. 333 

DADBY--DBET*XKA 333 

THETY-THETY+DADBY 334 

GO TO 350 335 

340 UEL3-UELS-UELX( J ) 336 

ZE-UEL3/U 337 

DBET-+C(ll)/57.3+(ZE-ZEl)/3. 338 

DADBY--DBET*XKA 339 

THETY-THETY+DADBY 340 
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CALL SUBROUTINE FOR COMPUTING MINinUM DISTURBING ACCELERATION FOR 341 


ONE SIDED LIMIT CVCLE 342 

250 CALL THEMIN ( THEDDA, PC ( 5 ) , PC ( 6 ) , PC ( 7 ) , TIP, T2P, THMIN ) 343 

IF (ABS(THETP )-THMIN .GE. 0 ) GO TO 260 344 

CALL CYCLE ( THEDDA, PC ( 5 ) ,PC (6 ), PC ( 7 ) ,T1P, T2P, DCP, AP, NERRl , PBFCT ) 345 

GO TO 270 346 

COMPUTE DUTY CYCLE AND INCREMENTAL IMPULSE 347 

260 DCP-ABS(THETP/THEDDA) 348 

270 DPI-DCP*PC(8)*DELTA 349 

CALL THEMIN ( PSDDA, YC ( 5 ) , YC (6 ), YC (7 ) ,T1Y,T2Y, THMIN ) 350 

IF (ABS(THETY)-THMIN .GE. 0 ) GO TO 280 351 

CALL CYCLE ( PSDDA, YC( 5 ) , YC ( 6 ) , YC( 7 ) ,T1Y,T2Y, DCY, AY,NERR1 , YBFCT ) 352 

GO TO 290 353 

280 DCY=ABS<THETY/PSDDA) 354 

290 DYI-DCY*YC(8 )*DELTA 355 

COMPUTE ROLL IMPULSE DUE TO CG OFFSET 356 

THXU = THX( J )»UGHT(J ) 357 

TULTLP»X1KB(J)/THXU 358 

TUULP-TOU(J)/(THXU*XINER(J) ) 359 

CNAUL-QUDFAC*X1KA(J)/THXU 360 

ALP»DALP+ALPHP 361 

BET=DBET+BETAP 362 

RM = TOU( J )C(SUP:*cETYA-SUY^cETPA )+QUDFAC»CNASQ( J )*(SUY*ALP+ 363 

1 SUP*BET)/UGHT(J)+SUY*(TULTLP)»:ETPA+TUULP*SUP+ 364 

2 CNAUL* ALP ) -SUP* ( TULTLP*ETYA+TUULP*SUY-CNAUL*BET ) 365 

COMPUTE INCREMENTAL ROLL IMPULSE 366 

RIMP»12.*ABS(RM)*DELTA/A(7) 367 

300 BIMP-BIMP+ABS(DPI)+ABS(DYI)+RIMP 368 

COMPUTE TOTAL BOOST IMPULSE AND FUEL CONSUMPTION 369 

BIMP»(BIMP+CAPIMP )*( 1 .+(A(5)*A(8)*ABS(BT(11 ) ) )/(57.3*A(7) ) )+RTI 370 
UB00ST(ITER)»*BIMP/'BT<5) 371 

SUB-SUB+UBOOST< ITER ) 372 

COMPUTE LIMIT CYCLE FUEL CONSUMPTION DURING COAST UITH FILTER 373 

IN USING BOOST CONTROL SYSTEM 374 
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THEDDA-PC(8)*( (A(4 )-BT< 14 )*TABL1 (9 ) )*A(8 )+A(5 )*A(9 ) )/( ia.*TABLl (8 )375 


1*BT<13)) 376 

PSDDA-THEDDA*VC<8 )/PC(8 ) 377 

PHIDD-RC<7)*A<7)/(6.*A(16)) 378 

CALL CYCLE ( THEDDA, PC < 5 ) , PC< 6 ) , PC (7 ),T1P, T2P, DCP, AP, NERRl , PBFCT ) 379 

NERl-NERl+NERRl 38© 

CALL CYCLE (PSDDA, YC ( 5 ), YC ( 6 ), YC( 7 ), T1 Y,TEY, DCY, AY, NERRl , YBFCT ) 381 

NERl-HERl+NERRl 382 

T1R»RC(1 )+RC(2) 383 

T2R»RC(1 )+RC(3) , 384 

CALL CYCLE ( PHIDD, RC ( 4 ) , RC< 5 ), RC( 6 ), T1R,T2R, DCR, AR, NERRl , 0 . ) 385 

NRRl-NRRl+NERRl 386 

C COMPUTE CONTROL FUEL FLOU RATE LIMIT CYCLE UITH BOOST CONTROLS 387 

UDFILT=(DCP*PC(8)+DCY*YC(8)+DCR*RC(7)X2. )/BT(5) 388 

> TTO-TBO-TUEB 389 

M C TEST FOR COAST FILTER OPTION 1-COAST FILTER OUT 2-COAST FILTER IN 390 

IF (I0(4)-2 ,GE. 0 ) GO TO 420 391 

COMPUTE FUEL UITH FILTER IN COAST FOR OPTION UITH FILTER SUITCHED 392 
OUT DURING COAST 393 

DELT-A(10)-TBO 394 

UFILT«UDFILT*DELT 395 

GO TO 530 396 

TEST FOR STAGE OPTION 2-NO COAST CONTROL 3-COAST CONTROL 397 

420 IF (I0(2)-3 .GE. 0 ) GO TO 490 398 

BOOST AND COAST CONTROL SYSTEM IS THE SAME 399 

TEST FOR COAST OPTION 1 -COAST TIME 2-COAST FUEL 400 

IF (I0(3)-2 .GE. 0 ) GO TO 460 401 

C STAGE 2 COAST TIME - FILTER IN 402 

TCOAST( ITER ) » < A ( 12 )-UBOOST < ITER ) )/UDFILT 403 

STC-STC+TCOASTCITER) 404 

C TEST FOR PRINT OPTION 1-INDIUIDUAL SAMPLE DATA, 2-ONLY STATISTICS405 

IF (I0(l)-2 .GE. 0 ) GO TO 600 406 

430 IF (NPKl .NE. 0 ) GO TO 450 407 

440 NPAGE-NPAGE+1 408 
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CALL PAGEHD (NRUN,NPAGE,NLINE ) 409 

URITE( 6,980) 410 

URITE( 6,1010) 411 

NLINE-NLINE+6 418 

NPKl-1 413 

450 IF (NLINE+3-68 .GT. 0 ) GO TO 440 414 

URITE( 6,1000) ITER,CAPIPlP,RTI,BIt1P,TUEB,TT0,UB00ST(ITER),BT(5) 415 

URITE( 6,1080) TC0AST(ITER),UDFILT,NER1,NRR1 416 

NLINE-NLINE+3 417 

GO TO 600 418 

C STAGE a COAST FUEL - FILTER IN 419 

460 UCOAS(ITER)-A( 11 )*UDFILT 430 

SUC-SUC+UCOAS( ITER) 431 

UTOT ( ITER ) =UBOOST ( ITER )+UCOAS ( ITER ) 433 

SUT-SUT+UTOTCITER) 433 

C TEST PRINT OPTION 1-INDIUIDUAL SAMPLE DATA, 3-ONLV STATISTICAL 434 

IF (lO(l)-a .GE. 0 ) GO TO 600 435 

IF (NPKl .NE. 0 ) GO TO 480 436 

470 NPAGE-NPAGE+1 437 

CALL PAGEHD (NRUN,NPAGE,NLINE ) 438 

URITE( 6,980) 439 

URITE( 6,1030) 430 

NLINE=NLINE+6 431 

NPK1=1 433 

480 IF (NLINE+3-63 .GT. 0 ) GO TO 470 433 

URITEC 6,1000) ITER,CAPIMP,RTI,BinP,TUEB,TT0,UB00ST(ITER),BT(5) 434 

URITE( 6,1040) UCOAS(ITER),UDFILT,NERl,NRRl,UTOT(ITER) 435 

NLINE-NLINE+3 436 

GO TO 600 437 

SEPARATE COAST CONTROL 10(3 )-3 438 

TEST FOR COAST OPTION 1-COAST TIME 3-COAST FUEL AND RETRO TIME 439 
490 IF (I0(3)-3 .GE. 0 ) GO TO 500 440 

C COMPUTE STAGE 3 COAST TIME - FILTER IN 441 


TCOAST( ITER )“ ( A( 13 )-UBOOST( ITER )-3.»(PC(8 )+YC(8 )+RC(7 ) )*TRETR0/'RET443 
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IRSD/UDFILT 443 

STC-STC+TCOAST(ITER) 444 

C TEST PRINT OPTION 1-INDIOIDUAL SAMPLE DATA, E-ONLV STATISTICAL 445 

IF(IOd)) 430,600,600 446 

C COMPUTE STAGE 3 COAST FUEL, RETRO TIME AND TOTAL FUEL FILTER IN 447 

500 UCOAS(ITER)-A(ll)*UDFILT 448 

URETRO=A( 12 )-UBOOST( ITER )-UCOAS( ITER) 449 

TRET ( ITER ) - URETRO^RETRSI / ( 2 . * ( PC ( 8 )+YC ( 8 ) +RC ( 7 ) ) ) 450 

SUC»SUC+UCOAS(ITER) 451 

STR=STR+TRET ( ITER ) 452 

UTOT( ITER ) = A ( 12 )-URETRO:»c ( 1 . -TRETRO/TRET ( ITER ) ) 453 

SUT-SUT+UTOT ( ITER ) 454 

C TEST PRINT OPTION 1-INDIUIDUAL SAMPLE DATA 2-ONLY STATISTICAL 455 

IF (I0(l)-2 .GE. 0 ) GO TO 600 456 

> IF (NPKl ,NE. 0 ) GO TO 520 457 

' 510 NPAGE-NPAGE+1 458 

CALL PAGEHD (NRUN,NPAGE,NLINE ) 459 

URITEC 6,980) 460 

URITE( 6,1030) 461 

NLINE«NLINE+6 462 

NPKl-1 463 

520 IF (NLINE+3-62 .GT. 0 ) GO TO 510 464 

URITEC 6,1000) ITER,CAPIMP,RTI,BIMP,TUEB,TT0,UB00ST(ITER),BT(5) 465 

URITEC 6,1040) UCOASCITER),UDFILT,NERl,NRRl,UTOTCITER) 466 

NLINE=NLINE+3 467 

GO TO 600 468 

FILTER OUT COAST OPTION 469 

TEST FOR STAGE OPTION 2-NO COAST CONTROL 3-COAST CONTROL 470 

530 IF CIOC2)-3 .GE. 0 ) GO TO 550 471 

C UITH BOOST CONTROL, FILTER OUT IN COAST 472 

T1P-PCC2)+PCC3) 473 

T2P-PCC2)+PCC4)+T2FLX 474 

T1V-YCC2)+VCC3) 475 

T2V-YCC2)+YCC4)+T2FLX 476 
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CALL CYCLE ( THEDDA, PC ( 5 ) , PC ( 6 ) , PC ( 7 ) , TIP, TEP, DCP, AT, NERRl , PCFCT ) 477 


NERC-NERC+NERRl 478 

CALL CYCLE ( PSDDA, YC (5 ) , YC ( 6 ), YC( 7 ) , T1 Y, TEY, DCY, AT, NERRl , YCFCT ) 479 

NERC-NERC+NERRl 480 

UDCST-(DCR*RC(7 ):t:E.+DCP*PC(8 )+DCY*YC(8) )/BT(5) 481 

C TEST FOR COAST OPTION 1-COAST TIME E-COAST FUEL AND RETRO TINE 48E 

IF (I0(3)-S .GE. 0 ) GO TO 540 483 

C CONFUTE COAST TINE 484 

TCOASTC ITER )-DELT+(A( IE )-UBOOST( ITER )-UFILT )/UDCST 485 

STC*STC+TCOAST(ITER) 486 

GO TO 600 487 

CONFUTE COAST FUEL AND BOOST+COAST FUEL 488 

COAST TINE IS ONLY PORTION UITH FILTER SUITCHED OUT; COAST FUEL 489 
INCLUDES FILTER IN AND FILTER OUT COAST TINES 490 

540 UCOAS(ITER)-A(ll)*UDCST+UFILT 491 

SUC-SUC+UCOAS(ITER) 49S 

UTOT( ITER ) = UBOOST< ITER )+UCOAS( ITER ) 493 

SUT»SUT+UTOT(ITER ) 494 

GO TO 600 495 

STAGE 3 OPTION 496 

SEPARATE COAST CONTROL FILTER OUT IN COAST AND YAU-ROLL NIXING 497 
E ROLL NOTORS FOR YAU; S ROLL NOTORS FOR ROLL 498 

550 THEDDA = THEDDA)(CPDC(E)/PC(8) 499 

PSDDA»PSDDA*E.*PDC(6)/YC(8) 500 

PHIDD-PHIDD!|:PDC(6)/RC(7) 501 

T1P-PC(E)+PDC(3) 50B 

TEP-PC(E)+PDC(4)+TEFLX 503 

T1Y-YC(E)+PDC(7) 504 

TSY-YC(E)+PDC(8)+TEFLX 505 

T1R-RC( 1)+PDCC7) 506 

TER-RC( 1)+PDC(8 ) 507 

C CONFUTE PITCH, YAU, ROLL TORQUEING AND DEADBAND REDUCTION FUEL 508 

PRTRU-(PDC<E )*(E.*TEP+ABS( (PTRC+57. 3*AP+E. » (PC (6 )-PDC ( 1) )/PC(5 ) )/(509 
157.3*THEDDA) ) ) ) 510 
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VRTRU -S . *PDC ( 6 )* < a . *T2V+ABS ( ( YTRC+57 . 3»AV+2 . * ( YC ( 6 )-PDC ( 9 ) )/VC ( 5 ) )51 1 
1/(57. 3*PSDDA )) ) 512 

RRTRU -2 . *PDC < 6 ) * ( 2 . *T2R+ ABS ( < 57 . 3*AR+2 . * ( RC ( 5 ) -PDC ( 5 ) ) /RC ( 4 ) )/ ( 57 . 513 
13»PHIDD))) 514 

C TEST COAST SYSTEM 10(7)- 0 YAU-ROLL MIXING, -1 NO MIXING 515 

IF (I0(7)-l ,LT. 0 ) GO TO 560 516 

C NO YAU-ROLL MIXING, SEPARATE PITCH AND VAU MOTORS 517 

PSDDA-PSDDA*YCF1/(2.*PDC(6) ) 518 

T1Y-YC(2)+YCT1 519 

T2V-VC(2)+YCT2+T2FLX 520 

YRTRU - YCF 1* ( 2 . *T2Y+ ABS ( ( YTRC+57 . 3CAY+2 . * ( YC ( 6 ) -PDC ( 9 ) )/YC ( 5 ) )/ ( 57 . 521 
13CPSDDA ) ) ) 522 

C COMPUTE FUEL REQUIRED FOR DEADBAND REDUCTION AND TORQUEING 523 

560 TRDRU-(PRTRU+YRTRU+RRTRU)/PDC(10) 524 

CALL CYCLE ( THEDDA, PC( 5 ) , PDC ( 1 ), PC(7 ),T1P,T2P, DCP, AT,NERR1 ,PCFCT ) 525 
NERC-NERC+NERRl 526 

CALL CYCLE ( PSDDA, YC (5 ) , PDC(9 ), YC(7 ),T1 Y,T2Y, DCY, AT, NERRl, YCFCT ) 527 

NERC-NERC+NERRl 528 

CALL CYCLE ( PHIDD,RC( 4 ), PDC( 5 ) ,RC(6 ),T1R,T2R, DCR, AT,NERR1, 0.0 ) 529 

NERC-NERC+NERRl 530 

C TEST FOR COAST CONTROL SYSTEM 0-YAU-ROLL MIXING, 1- NO MIXING 531 

IF (I0(7)-l .GE. 0 ) GO TO 570 532 

C YAU-ROLL MIXING 533 

UDCST- ( PDC ( 2 )»DCP+2 . *PDC ( 6 )* ( DCY+DCR ) )/PDC( 10 ) 534 

GO TO 580 535 

570 UDCST-(PDC(2 )*DCP+2.*PDC(6)»DCR+YCF1»DCY)/PDC( 10) 536 

C TEST FOR COAST OPTION 1-COAST TIME 2-COAST FUEL AND RETRO TIME 537 

580 IF (I0(3)-2 .GE. 0 ) GO TO 590 538 

C COAST TIME OPTION 539 

TCOAST( ITER )-DELT+(A( 12 )-UBOOST( ITER )-UFILT-TRDRU-(PC(8)+YC(8)+RC(540 
17) )*2.»TRETR0/RETRSI )/UDCST 541 

STC-STC+TCOAST(ITER) 542 

GO TO 600 543 

C STAGE 3 COAST FUEL, RETRO TIME, TOTAL FUEL UITH INPUT RETRO TIME 544 
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590 UC0AS(ITER)-A(11)*UDCST+UFILT 545 

SUC-SUC-t-UCOAS(ITER) 546 

URETRO-AC la )-UBOOST( ITER )-TRDRLI-UCOAS( ITER ) 547 

TRET ( ITER ) -URETRO*RETRSI / ( a . * ( PC < 8 )+YC ( 8 )+RC ( 7 ) ) ) 548 

STR-STR+TRET(ITER ) 549 

UTOT ( ITER ) ■ A ( la )-URETRO* ( 1 , -TRETRO/TRET< ITER ) ) 550 

SUT-SUT+UTOT(ITER) 551 

C UPDATE COUNTER ON NUMBER OF DEADBAND OUERSHOOT CASES 553 

600 IF (NERl-1 .LT. 0 ) GO TO 610 553 

MERl-MERl+l 554 

610 IF (NRRl-1 .LT. 0 ) GO TO 630 555 

MRERl»nRERl+l 556 

630 IF (NERC-1 .LT. 0 ) GO TO 630 557 

nERC»MERC+l 558 

C TEST FOR COAST FILTER OPTION 1-COAST FILTER OUT a-COAST FILTER IN 559 

630 IF (I0(4)-a .GE. 0 ) GO TO 750 560 

C TEST PRINT OPTION 1-INDIUIDUAL SAMPLE RESULTS, a-ONLV STATISTICAL 561 

IF (lO(l)-a .GE. 0 ) GO TO 750 56a 

C TEST FOR STAGE OPTION a-BOOST CONTROL COAST 3-COAST CONTROL 563 

IF (I0(a)-3 .GE. 0 ) GO TO 690 564 

C TEST FOR COAST OPTION 1-COAST TIME a-COAST FUEL AND RETRO TINE 565 

IF (I0(3)-a .GE. 0 ) GO TO 660 566 

IF (NPKl .NE. 0 ) GO TO 650 567 

C PRINT COAST TIME 568 

640 NPAGE>=NPAGE-<-l 569 

CALL PAGEHD (NRUN,NPAGE, NLINE ) 570 

URITE( 6,980) 571 

URITEC 6,1050) 573 

NLINE-NLINE+6 573 

NPKl-1 574 

650 IF (NLINE+3-63 .GT. 0 ) GO TO 640 575 

URITEC 6,1000) ITER,CAPIMP,RTI,BIMP,TUEB,TT0,UB00ST(ITER),BT(5) 576 
URITE( 6,1060) TC0AST(ITER),UDFILT,UDCST,NER1,NRR1 577 

NLINE-NLINE+3 578 
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GO TO 750 579 

C PRINT COAST FUEL AND RETRO TIME 580 

660 IF (NPKl .NE. 0 ) GO TO 680 581 

670 NPAGE-NPAGE+1 588 

CALL PAGEHD (NRUN,NPAGE,NLINE ) 583 

URITEC 6,980) 584 

URITE< 6,1070) 585 

NLINE-NLINE+6 586 

NPKl-1 587 

680 IF (NLINE+3-6B .GT. 0 ) GO TO 670 588 

URITE< 6,1000) ITER,CAPinP,RTI,BIMP,TUEB,TT0,UB00ST(ITER),BT(5) 589 
URITE( 6,1080) UCOAS(ITER),UDFILT,UDCST,NERl,NRRl,NERC,UTOT(ITER) 590 
NLINE-NLINE+3 591 

GO TO 750 593 

C TEST FOR COAST OPTION 1-COAST TINE S-COAST FUEL AND RETRO TINE 593 

690 IF (I0(3)-S .GE. 0 ) GO TO 730 594 

C PRINT COAST TINE 595 

IF (NPKl .NE. 0 ) GO TO 710 596 

700 NPAGE«NPAGE+1 597 

CALL PAGEHD (NRUN,NPAGE,NLINE ) 598 

URITEC 6,980) 599 

URITEC 6,1090) 600 

NLINE-NLINE+6 601 

NPKl-1 603 

710 IF CNLINE+3-63 .GT. 0 ) GO TO 700 603 

URITEC 6,1000) ITER, CAPINP,RTI,BIMP,TUEB,TTO,UBOOST( ITER), BTC5) 604 

URITEC 6,1100) TC0ASTCITER),UDFILT,TRDRU,UDCST,NER1,NRR1,NERC 605 

NLINE-NLINE+3 606 

GO TO 750 607 

C PRINT COAST FUEL AND RETRO TIME 608 

730 IF CNPKl .NE. 0 ) GO TO 740 609 

730 NPAGE-NPAGE+1 610 

CALL PAGEHD (NRUN,NPAGE,NLINE ) 611 

URITEC 6,980) 613 
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URITE( 6,1110) 613 

NLINE-NLINE+6 614 

NPKl-1 615 

740 IF (NLINE+3-ea .GT. 0 ) GO TO 730 616 

URITEC 6,1000) ITER,CAPIPlP,RTI,BinP,TUEB,TT0,UB00ST(ITER),BT(5) 617 
URITE( 6,11H0) UCOASC ITER ),UDFILT,TRDRU,UDCST, TRET ( ITER ),NER1,NRR1618 
l,NERC,UTOT(ITER) 619 

NLINE-NLINE+3 620 

750 CONTINUE 621 

622 

END OF INDIUIDUAL SAMPLE CALCULATION LOOP 623 

624 

BEGIN PART 2 STATISTICAL ANALYSIS 625 

AUI«IT 626 

UI-AUI+1. 627 

UIB-AUI-1. 628 

UBMEAN-SUB/AUI 629 

C TEST FOR COAST OPTION 1-COAST TIME 2-COAST FUEL AND RETRO TIME 630 

IF (I0(3)-2 .GE. 0 ) GO TO 790 631 

C COAST TIME OPTION 632 

TCMEAN-STC/AUI 633 

C CALL SUBROUTINE FOR REARRANGING ARRAYS 634 

CALL ASCEND ( IT, UBOOST, 0 ) 635 

CALL ASCEND ( IT, TCOAST, 1 ) 636 

SDXUB-0. 637 

SDXTC=0, 638 

C COMPUTE DISCRETE SAMPLE PROBABILITY LEUEL,SUM UALUES AND SQUARES 639 

DO 780 KAB-1,IT 640 

XKAB-KAB 641 

PROB-XKAB/UI 642 

DXUB-(UBMEAN-UBOOST(KAB) ):r*2 643 

DXTC-(TCMEAN-TCOAST(KAB ) )*)|:2 644 

SDXUB-SDXUB+DXUB 645 

SDXTC-SDXTC+DXTC 646 
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C PRINT ELEMENTS OF REARRANGED ARRAY AS COMPUTED 647 

IF (NPKa .NE. 0 ) GO TO 770 648 

760 NPAGE-NPAGE+1 649 

CALL PAGEHD ( NRUN, NPAGE> NLINE ) 650 

URITEC 6,1130) 651 

NLINE-NLlNE+4 65E 

NPK3-1 653 

770 IF (NLINE+1-6E .GT. 0 ) GO TO 760 654 

URITEC 6,1140) PROB,UBOOST(KAB),TCOAST(KAB) 655 

NLINE-NLINE+1 656 

780 CONTINUE 657 

PRINT MEAN, SIGH A, LAST RANDOM SEQUENCE INTEGER, NO. OF DEADBAND 0UER658 

SHOOT CASES 659 

SIGB-SQRTCSDXUB/UIB) 660 

SIGT-SQRTCSDXTC/UIB) 661 

NPAGE-NPAGE+1 663 

CALL PAGEHD ( NRUN, NPAGE, NLINE ) 663 

URITEC 6,1150) 664 

URITEC 6,1160) UBMEAN,TCMEAN 665 

URITEC 6,1170) SIGB,SIGT 666 

C FIND LAST RANDOM SEQUENCE INTEGER, CDC RANF RANDOM NUMBER GENERAT0R667 

CALL RANGET CIO 668 

URITEC 6,1180) K 669 

URITEC 6,1190) 670 

URITEC 6,1300) MERl 671 

URITEC 6,1310) MRERl 673 

URITEC 6,1330) MERC 673 

GO TO 10 674 

C OPTION FOR FUEL CONSUMPTION 675 

790 UCMEAN-SUC/AUI 676 

UTOTM-SUT/AUI 677 

C TEST FOR STAGE OPTION 3-COAST UITH CONTROL 3-COAST CONTROL 678 

IF CI0<3)-3 .LT. 0 ) GO TO 870 679 

C SEPARATE COAST CONTROL OPTION 680 
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TRETM-STR/AUI 681 

CALL ASCEND ( IT, UBOOST^C ) 688 

CALL ASCEND (IT,UCOAS,0) 683 

CALL ASCEND (IT,UTOT,0) 684 

CALL ASCEND (IT, TRET, 1) 685 

SDXUB-0. 686 

SDXFC-0. 687 

SDXUT-0. 688 

SDXRT-0. 689 

C CALCULATE PROBABILITIES OF FUEL CONSUMPTION AND RETRO TIME 690 

DO 880 KAB=1,IT 691 

XKAB-KAB 698 

PROB-XKAB/UI 693 

DXUB»(UBMEAN-UBOOST(KAB ) )**8 694 

DXFC“(UCMEAN-UCOAS(KAB) 695 

DXUT“(UTOTM-UTOT(KAB) )**8 696 

DXRT-(TRETM-TRET(KAB) )**8 697 

SDXUB=«SDXUB+DXUB 698 

SDXFC=SDXFC+DXFC 699 

SDXUT-SDXUT+DXUT 70Q 

SDXRT-SDXRT+DXRT 701 

C PRINT ELEMENTS OF REARRANGED ARRAY AS COMPUTED 708 

IF (NPK3 .NE. 0 ) GO TO 810 703 

800 NPAGE=NPAGE+1 704 

CALL PAGEHD ( NRUN, NPAGE, NLINE ) 70S 

URITE( 6,1830) 706 

NLINE-NLINE+4 707 

NPK3-1 708 

810 IF (NLINE+1-68 .GT. 0 ) GO TO 800 709 

URITE( 6,1840) PROB,UBOOST(KAB),UCOAS(KAB),TRET(KAB),UTOT(KAB) 710 
NLINE-NLINE+1 711 

880 CONTINUE 718 

SIGB-SQRT(SDXUB/UIB ) 713 

SIGUC-SQRT(SDXFC/UIB) 714 
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SIGUT-SQRT(SDXUT/UIB ) 715 

SIGRT-SQRT(SDXRT/UIB ) 716 

IF (IT-100 .LT. 0 ) GO TO 850 717 

ANI-0,995*UI 718 

NI-ANI 719 

ANIP-NI 730 

IF ( .5-(ANI-ANIP) .GT. 0 ) GO TO 840 731 

NI-NI+1 733 

840 CONTINUE 733 

C COMPUTE UPPER 9SH CONFIDENCE LEUEL ON 99, SX PROBABILITY LEUEL 734 

DEN-SQRT(1»-1.96*SQRT(3./(A0I-1. ))) 735 

XK-l./DEN 736 

XUB95“UBMEAN+XK*(UB00ST(NI )-UBMEAN) 737 

XUC95“UCMEAN+XK*(UC0AS(NI )-UCMEAN) 738 

XUT95»UTOTM+XK*(UTOT(NI )-UTOTM) 739 

XRET95«TRETM+XK:(c(TRET(NI )-TRETM) 730 

850 NPAGE-NPAGE+1 731 

CALL PAGEHD (NRUN,NPAGE,NLINE ) 733 

URITE( 6,1340) 733 

URITE( 6,1350) 734 

URITE( 6,1340) 735 

URITEC 6,1350) 736 

URITE( 6,1360) UBMEAN,UCMEAN,TRETM,UTOTM 737 

URITE( 6,1370) SIGB,SIGUC,SIGRT,SIGUT 738 

IF (IT-100 .LT, 0 ) GO TO 860 739 

URITE( 6,1310) XUB95,XUC95,XRET95,XUT95 740 

860 CONTINUE 741 

C RETRIEUE LAST RANDOM SEQUENCE INTEGER 743 

CALL RANGET (K) 743 

URITE( 6,1180) K 744 

URITE( 6,1190) 745 

URITE( 6,1300) flERl 746 

URITE( 6,1310) MRERl 747 

URITE( 6,1330) MERC 748 
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COMPUTE STATISTICS AND COMPARE UITH NORMAL, LOG-NORMAL, AND UEIBULL 749 


STATISTICAL DISTRIBUTION FUNCTIONS 750 

CALL UBL ( IT,UBOOST,NUMCHI,FUB) 751 

CALL UBL (IT,UCOAS,NUnCHI,FUC) 753 

CALL UBL (IT,UTOT,NUMCHI,FUT) 753 

URITE( 6,1400) 754 

URITE( 6,1410) ( (FUB< J),FUC( J),FUT(J) ),J-1,4),FUB(9),FUC(9),FUT(9)755 
1,FUB(18),FUC(18),FUT(18) 756 

URITEC 6,1430) 757 

URITEC 6,1410) ( (FUB(J),FUC(J),FUT(J) ),J-5,8),FUB(10),FUC(10),FUT(758 
110),FUB(19),FUC(19),FUT(19) 759 

URITE( 6, 1430) < (FUB( J ),FUC( J ),FUT< J ) ), J-11, 17 ),FUB(30 ),FUC(30 ),FU760 
1T(30) 761 

C TEST FOR HISTOGRAM OUTPUT OPTION 10(9 )-0 HISTOGRAMS, -1 NONE 763 

IF (10(9). EQ.l) GO TO 10 763 

C CALCULATE AND PLOT HISTOGRAMS 764 

NDX-30 765 

C COMPUTE RANGE OF UARIABLES AND PLOT HISTOGRAM FOR RETRO TIME 766 

CALL RANGE (TRET, IT, GREAT, SMALL, DX,NDX,ERRTB ) 767 

IF (ERRTB .GT. 0 ) GO TO 950 768 

URITE( 6,1390) 769 

CALL HISTO (TRET, IT, GREAT, SMALL, DX) 770 

GO TO 950 771 

C STAGE 3- SEPARATE COAST CONTROL OPTION 773 

870 CALL ASCEND ( IT,UB00ST,0 ) 773 

CALL ASCEND (IT,UCOAS,0) 774 

CALL ASCEND (IT,UTOT,0) 775 

SDXUB-0. 776 

SDXFC=0. 777 

SDXUT-0. 778 

C COMPUTE PROBABILITY DISTRIBUTION AND SUM OF SQUARES 779 

DO 900 KAB-1,IT 780 

XKAB-KAB 781 

PROB-XKAB/UI 783 
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DXUB-(UBHEAN-UBOOST(KAB) 783 
DXFC- (UCMEAN-UCOASCKAB ) 784 

DXUT-(UTOTM-UTOT(KAB) 785 

SDXUB-SDXUB+DXUB 786 

SDXFC-SDXFC+DXFC 787 

SDXUT-SDXUT+DXUT 788 

C PRINT ELEMENTS OF REARRANGED ARRAY AS COMPUTED 789 

IF (NPK4 ,NE« 0 ) GO TO 890 790 

880 NPAGE-NPAGE+1 791 

CALL PAGEHD (NRUN,NPAGE,NLINE ) 792 

URITE( 6,1280) 793 

NLINE=NLINE+4 794 

NPK4-1 795 

890 IF (NLINE+1-62 .GT. 0 ) GO TO 880 796 

URITE( 6,1290) PROB,UBOOST(KAB),UCOAS(ICAB),UTOT<KAB) 797 

NLINE=NLINE+1 798 

900 CONTINUE 799 

C COMPUTE STANDARD DEUIATION 800 

SIGB=SQRT(SDXUB/UIB) 801 

SIGUC=SQRT(SDXFC/UIB ) 802 

SIGUT»SQRT(SDXUT/UIB) 803 

IF (IT-100 ,LT. 0 ) GO TO 930 804 

ANI=>0«9g5:|;UI 805 

NI°-ANI 806 

ANIP=NI 807 

IF ( .5-(ANI-ANIP ) .GT. 0 ) GO TO 920 808 

NI=NI+1 809 

920 CONTINUE 810 

C COMPUTE UPPER 9S% CONFIDENCE LEUEL, 99.5J( PROBABILITY LEUEL 811 

DEN-SQRT(1 .-1 .DSJlcSQRTta./'CAUI-l . ) ) ) 812 

XK-l./DEN 813 

XUB95>=UBnEAN+XK*(UB00ST(NI )-UBMEAN) 814 

XUC95-UCMEAN+XK*(UC0AS(NI )-UCMEAN ) 815 

XUT95-UT0TM+XK*(UT0T(NI )-UTOTM) 816 
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930 


> 

I 

iv> 

a\ 


NPAGE-NPAGE+1 

CALL PAGEHD ( NRUN, NPAGE,NLINE ) 
URITE( 6^1340) 


URITE( 

URITE( 

URITE( 

LIRITE( 

URITE( 


6>1350) 
6,1340) 
6, 1300) 
6,1360) 
6, 1370) 


940 


UBnEAN,UCnEAN,UTOTM 
SIGB,SIGUC,SIGUT 
IF (IT-100 .LT, 0 ) GO TO 940 
URITEC 6,1310) XUB95,XUC95,XUT95 
CONTINUE 

FIND LAST RANDOM SEQUENCE INTEGER 
CALL RANGET (K) 

URITEC 6,1180) K 


WRITE ( 
URITE( 
URITE( 
URITE( 
COMPUTE 
UEIBULL 


6,1190) 

6,1300) 

6,1310) 

6,1330) 


MERl 
MRERl 
MERC 

STATISTICAL PARAMETERS AND 
DISTRIBUTION FUNCTIONS 


CALL UBL (IT,UBOOST,NUMCHI,FUB) 

CALL UBL (IT,UCOAS,NUMCHI,FUC) 

CALL UBL (IT,UTOT,NUMCHI,FUT ) 

URITEC 6,1400) 

URITEC 6,1410) C CFUBCJ),FUCC J),FUT(J) ), J 
1,FUB(18),FUCC18),FUTC18) 

URITEC 6,1430) 


817 

818 
819 

830 

831 
833 

833 

834 

835 

836 

837 

838 

839 

830 

831 
833 

833 

834 

FIT TO NORMAL, LOG-NORMAL, AND 835 

836 

837 

838 

839 

840 

1,4),FUB(9),FUCC9),FUT(9)841 

843 
843 


URITEC 6,1410) C (FUB(J),FUC(J),FUT(J) ), J-5,8),FUB(10),FUC(10),FUT(844 
110),FUB(19),FUC(19),FUT(19) 845 

URITEC 6,1430) ( C FUB ( J ) , FUC ( J ) , FUT ( J ) ), J« 1 1 , 17 ) ,FUB C 30 ) , FUC ( 30 ) ,FU846 
1TC30) 847 

TEST FOR HISTOGRAM OUTPUT 10(9 )-0 GIUES HISTOGRAMS, -1 NONE 848 

IF CI0C9).EQ.l) GO TO 10 849 

CALCULATE AND PLOT HISTOGRAMS OF BOOST, COAST, AND TOTAL FUEL 850 
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NDX-30 851 

950 CONTINUE 858 

CALL RANGE (UBOOST, IT, GREAT, SMALL, DX,NDX,ERRTB ) 853 

IF (ERRTB .GT. 0 ) GO TO 960 854 

URITE( 6,1360) 855 

CALL HISTO (UBOOST, IT, GREAT, SMALL, DX) 856 

960 CALL RANGE ( UCOAS, IT, GREAT, SMALL, DX, NDX, ERRTB ) 857 

IF (ERRTB .GT. 0 ) GO TO 970 858 

URITE( 6,1370) 859 

CALL HISTO (UCOAS, IT, GREAT, SMALL, DX) 860 

970 CALL RANGE (UTOT, IT, GREAT, SMALL, DX, NDX, ERRTB ) 861 

IF (ERRTB .GT. 0 ) GO TO 10 862 

URITE( 6,1380) 863 

CALL HISTO (UTOT, IT, GREAT, SMALL, DX) 864 

GO TO 10 865 

866 

END OF MAIN ROUTINE ONLY FORMAT STATEMENT FOLLOU 867 

868 

980 FORMAT (/85H SMPL CAPT IMP ROLL TORQ BOOST IMP T(UEB) 869 

1 T(TO) BOOST FUEL ISP/8X,77H(LB-SEC ) IMP(LB-SC) (LB-S870 

2EC) (SEC) (SEC) (LBS) (SEC)) 871 

990 FORMAT (8F10.3) 872 

1000 FORMAT (/I5, 1P7E12.4 ) 873 

1010 FORMAT (8X,44HCST TIME FLOU RATE NO. OF DB OUERSHOOTS, /10X, 39874 
IH(SEC) FILTER IN BOOST P-Y ROLL) 875 

1020 FORMAT (5X,1P2E12.4,I8,I12) 876 

1030 FORMAT (8X,43HCST FUEL FLOU RATE NO. OF DB OUERSHOOTS, 3X, 10HTO877 
ITAL FUEL,/,10X,40H(LBS) FILTER IN BOOST P-Y R0LL,7X,5H(L878 

2BS),/) 879 

1040 FORMAT (5X,1P2E12.4,I8,I12,4X,1PE12.4) 880 

1050 FORMAT (8X,55HCST TIME FLOU RATE FLOU RATE NO. OF DB 0UERSH0881 

1OTS,/,10X,53H(SEC ) FILTER IN COAST BOOST P-Y ROLL 882 

2 / ) 883 

1060 FORMAT (5X,1P3E12. 4, 18,112) 884 
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1070 FORMAT (8X,66HCST FUEL FLOU RATE FLOU RATE NUMBER 

lAND OUERSHOOTS,4X,10HTOTAL FUEL^/,8X,78H(P0UNDS ) FILTER 


OF 

IN 


aOAST 

1080 FORMAT 
1090 FORMAT 
IBER OF 
EAST 

1100 FORMAT 
1110 FORMAT 

ITIME NO. OF DB 
2 IN FUEL(LBS) 


BOOST P-Y ROLL COAST P-V (LBS) ,/) 

(5X, 1P3E12.4, 18. 21 12, 4X, 1PE12.4) 

(8X.69HCST TIME FLOU RATE DB RED TORO 
OUERSHOOTS,/,7X.81H(SECONDS) FILTER IN 
BOOST P-Y ROLL COAST P-Y,/) 

(5X,1P4E12.4, 18,2112) 

(8X,93HCST FUEL FLOU RATE DB RED TORO 


FLOU RATE 
FUEL(LBS) 


FLOU RATE 


OUERSHOOTS 

COAST 


TOTAL FUEL, /,8X,89H( POUNDS 
(SEC) BOOST P-Y ROLL CST 


FUEL,7X, 10HCOAST 


3 (LBS) ,/) 

1120 FORMAT (5X,1P5E12.4,I8,2I6,4X,1PE12.4) 

1130 FORMAT ( //5X, 1 IMPROBABILITY, 6X, 10HBOOST 
1, 8H(P0UNDS ),9X,9H( SECONDS ),/ ) 

1140 FORMAT ( 5X, F10 . 5, F16 . 3, F19 . 4 ) 

1150 FORMAT (//22X, 10HBOOST FUEL, 7X, 10HCOAST TIME) 

1160 FORMAT (/5X,4HMEAN,11X,F11.3,F19.4) 

1170 FORMAT (/5X,14HSTD. DEUI ATION, F12 .3, F19. 4 ) 

1180 FORMAT (//5X,42HFINAL UALUE OF RANDOM SEQUENCE INTEGER 
1190 FORMAT (/5X,38HN0. OF SAMPLES UITH DEADBAND OUERSHOOT ) 
1200 FORMAT (/10X,25HBOOST(PITCH AND YAU) »,I8) 

1210 FORMAT (/15X,20H(ROLL ) -,I8) 

1220 FORMAT (/10X,25HCOAST =,I8) 

1230 FORMAT ( /5X, 1 IMPROBABILITY, 6X, 10HBOOST FUEL,5X, 10MCOAST 
IHRETRO 


DEADB885 
C886 

887 

888 
NUM889 

CO890 

891 

892 

RETRO 893 
FILTER894 
P-Y 895 

896 

897 

TIME,/23X898 

899 

900 

901 

902 

903 

> ,120) 904 

905 

906 

907 

908 

FUEL.5X, 10909 


TIME,6X, 10HTOTAL FUEL, /,23X, 8H(P0UNDS ), 7X, 8H( POUNDS ), 7X, 9M ( 910 
2SEC0NDS),7X,8H( POUNDS),/) 911 

1240 FORMAT ( 5X, F10 . 5, F16 . 3, F15 . 3 , F16 . 4, F15. 3 ) 912 

1250 FORMAT ( 22X, 10HBOOST FUEL,5X, 10HCOAST FUEL, 5X, 10HRETRO TIME, 6X, 10H913 
ITOTAL FUEL) 914 

1260 FORMAT ( /5X, 4HMEAN, 1 IX, Fl 1 . 3, F15 . 3, F16. 4,F15 . 3 ) 915 

1270 FORMAT (/5X,14HSTD. DEUIATION, F12 . 3, F15 . 3, F16 . 4, F15 . 3 ) 916 

1280 FORMAT ( /5X, 1 IMPROBABILITY, 6X, 10HBOOST FUEL,5X, 10HCOAST FUEL, 6X, 10917 
IMTOTAL FUEL,/,23X,8M(P0UNDS),7X,8H(P0UNDS),8X,8H(P0UNDS) ) 918 
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ia90 FORMAT (5X, F10 . B, F16 . 3, FIS . 4>F15 . 3 ) 919 

1300 FORMAT (/aaX.EBHBOOST FUEL COAST FUEL, 6X, 10HTOTAL FUEL) 9a0 

1310 FORMAT (/5X, )*:0 . 995 PROB*, 5X, FI 1 . 3,F15* 3, F16» 4,F15 . 3, /, 5X, *0 . 95 CONBEl 
iFiD^t;) 9aa 

1330 FORMAT (118) 933 

1330 FORMAT (1115) 934 

1340 FORMAT (//) 935 

1350 FORMAT (73H 936 

1 /73H 937 

3 ) 938 

1360 FORMAT ( IHl , 1 IX, *BOOST FUEL CONSUMPTION HISTOGRAM*,/) 939 

1370 FORMAT ( IHl , 1 IX, *COAST FUEL CONSUMPTION HISTOGRAM*,/) 930 

1380 FORMAT ( IHl , 1 IX, *TOTAL FUEL CONSUMPTION HISTOGRAM*,/) 931 

1390 FORMAT ( IHl , 1 IX, *RETRO TIME HISTOGRAM*,/) 933 

1400 FORMAT ( IHl ,5X, *DISTRIBUTION FUNCTION FIT*, /34X, *BOOST COAST 933 

1 TOTAL*, //5X,*N0RMAL DISTRIBUTION*,/) 934 

1410 FORMAT ( 10X, *MEAN*, 16X, 3F10 . 3, /10X, *STANDARD DEU*,8X,3F10.3,/10X,*935 
ISKEUNESS*, 13X,3F10.3,/10X,*KURTOSIS*,iaX,3F10.3,/10X, *0.995 PROB 0936 
3.95 CONF*,3F10.3,/10X,*LEUEL OF SIGNIFICNCE*, 3F10 . 3 ) 937 

1430 FORMAT (/,5X,*L0G NORMAL DISTRIBUTION*, /6X, *(UALUES OF LOG)*) 938 

1430 FORMAT (/,5X, *UEIBULL DISTRIBUTION*, /10X, *PARAMETERS A*,8X, 3F10.4,939 
l/aiX,*B*,8X,3F10.4,/31X,*C*,8X,3F10.4,/10X,*PROBABILITY 0 . 5000*, 3X940 
a,3F10.3,/33X,*0.9900 *, 3F10 . 3, /EEX, *0. 9950 *, 3F10. 3, /33X, *0.9990941 

3 *,3F10.3,/10X,*LE<JEL OF SIGNIFICANCE*, 3F10 . 3 ) 943 

1440 FORMAT (3(E10.3)) 943 

1450 FORMAT (E15.5) 944 

1460 FORMAT (3(E15.5)) 945 

1470 FORMAT (I5/(6F10.3 ) ) 946 

END 947 



o o 


*DECK ALTI 948 

SUBROUTINE ALTI ( ALTO,UO,QO ) 949 

THIS SUBROUTINE COMPUTES THE ALTITUDE (ALTO ) FROM THE 950 

UELOCITY AND DYNAMIC PRESSURE USING THE ALTITUDE-DENSITY TABLE 951 

DIMENSION DUM(ia) 952 

C0MM0N/DRH/RH0(18),MR 953 

MR-1 954 

RHI-ALOG(UOlcUO/(2,3»:QO) ) 955 

C EXCHANGE THE ORDINATES AND ABSCISSAS IN THE ALTITUDE-DENSITY TABLE956 

DO 10 J-1,18,2 957 

DUM(J)»RHO( J + 1 ) 958 

10 DUM( J+1 ) = RHOC J) 959 

CALL TBLU ( 18, ALTO,RHI , DUM,MR ) 960 

RETURN 961 

END 962 
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^tcDECK ASCEND 963 

SUBROUTINE ASCEND (L,UAL,M) 964 

SUBROUTINE FOR REARRANGING ARRAY IN ASCENDING ORDER 965 

OR DESCENDING ORDER IF N-1 966 

DIMENSION UAL(l) 967 

K=L~1 968 

DO 40 J-1,K 969 

KB-J+1 970 

DO 40 JL=KB,L 971 

IF(M) 10,10,30 973 

10 IF (UAL(J)-UAL( JL) .LE. 0 ) GO TO 40 973 

GO TO 30 974 

30 IF(UAL( JL)-UAL( J ).LE.0) GO TO 40 975 

30 TEMP-UAL(J) 976 

UAL( J )-OAL(JL) 977 

UAL (JL) “TEMP 978 

40 CONTINUE 979 

RETURN 980 

END 981 
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*DECK CALCU 983 

SUBROUTINE CALCU 983 

THIS SUBROUTINE COMPUTES THE STATE UARIABLES AND THE FIRST 984 

DERIUATIUES FOR USE UITH THE RUNGE SUBROUTINE 985 

COnnON/DDI/TBLa(90 ),TBL3(90),NTHT,NUT,CDS,na,M3,UBO 986 

C0P1M0N/'DDY/YD0T(4),Y(4),T,DT 987 

C0tiri0N/DRH/RH0(18),nR 988 

DATA G0,R0,C0N/3a, 174,a0919668.,57.a9577951/' 989 

CALL TBLU (NTHT,TH,T,TBLa,fia) 990 

CALL TBLU ( NUT, UT, T, TBL3, M3 ) 991 

ALT=Y(a) 99a 

CALL TBLU ( 18, RHI , ALT, RHO, MR ) 993 

Y(4)-0.5*Y(1 ):|CY(1)/EXP(RHI ) 994 

G-G0/(1.+Y(a)/R0)«*3 995 

GAMA=Y(3)/C0N 996 

YDOT ( 1 ) -GO* ( TH-CDS*Y < 4 ) )/ ( UT+UBO )-G*SIN ( GAMA ) 997 

YDOT(3)=-CON*(G*COS(GAMA) )/Y(l ) 998 

YDOT(a)-Y(l )*SIN(GAMA) 999 

RETURN 1000 

END 1001 
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:i:DECK CHISQ 1008 

SUBROUTINE CHISQ ( N> A> CR« NQ, An> AS, UC, LD ) 1003 

THIS SUBROUTINE TEST THE SAMPLE DISTRIBUTION FOR GOODNESS 1004 

OF FIT TO A NORMAL, LOG-NORMAL, AND UEIBULL DISTRIBUTION FUNCTION 1005 

DIMENSION A(1000),JCEL(100),Z(47),C(47) 1006 

C0MM0N/C2/2,C 1007 

DATA <2(1 ), 1=1, 47)/ 1008 

1— 1000 , , -5» ,— 4. ,— 3*4, —3. 2, —3. , — S ♦ 9, — S. 8, — 3 » 7, —3 • 6 1009 

3, —3 •5,-3 ♦ 3, — 3 ♦ 1,— 3.0,— 1 .8,-1 .6,-1 .4,-1 .3,-1 .0,-0 .8,-0 «6,— 0.4, —0.3, 1010 

30 . . 0.3. 0.4. 0.6. 0.8. 1 .0. 1 .3. 1 .4. 1 .6. 1 .8. 3. 0.3. 1.3.3.3.5.3.6.3.7.3.81011 

4.3.9.3.0. 3.3.3.4.4.0.5.0.1000.0/ 1013 

DATA (C(I ), 1-1,47)/ 1013 

1 0., 0.0000006, 0.00006, 0.0003, 0.0007, 0.0013, 0.0019 1014 

3.0. 0036.0.0035.0.0047.0.0063.0.0107.0.0179.0.0338.0.0359.0.0548, 1015 

3 0.0808,0.1151,0.1587,0.3119,0.3743,0.3446,0.4307,0.5000,0.5793, 1016 

4 0.6554,0.7357,0.7881,0.8413,0.8849,0.9193,0.9453,0.9641,0.9773, 1017 

5 0.9831,0.9893,0.9938,0.9953,0.9965,0.9974,0.9981,0.9987,0.9993, 1018 

6 0.9997,0.99994,0.9999994,1.0/ 1019 

PN-N 1030 

NPARM-4 1031 

IF(LD.LE.0) NPARM-3 1033 

XN-NG 1033 

AB=A(10) 1034 

AT=A(N-9) 1035 

DU- (AT-AB)/(XN-3. ) 1036 

DO 10 J-1,NG 1037 

10 JCEL(J)=0 1038 

JCEL(1)-10 1039 

JCEL(NG)-10 1030 

NR-N-10 1031 

C COUNT THE NUMBER OF UALUES IN EACH GROUP 1033 

DO 40 J-11,NR 1033 

X-A(J)-AB 1034 

IF(X) 30,30,30 1035 
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c 


c 


B0 

30 

40 


50 


IC-1 

1036 

GO TO 40 

1037 

IC-X/DO+a 

1038 

IF(IC.GT.NG) IC-NG 

1039 

JCEL( IC)“JCEL( lO + l 

1040 

CALCULATE PROBABILITIES UITH EACH GROUP AND COMPARE 

ACTUAL UITH EX1041 

SPECTED FREQUENCV AND COMPUTE CHI-SQUARED UALUE CS 

1048 

COMPUTE LOUER CELL PROBABILITV 

1043 

XL“-1 .E10 

1044 

XU“AB 

1045 

n»i 

1046 

FE=PN*PZC ( XL, XU, AM, AS, UA, UB, UC, LD, M ) 

1047 

FA=JCEL(1 ) 

1048 

CS- ( (FA-FE )*(FA-FE ) )/FE 

1049 

COMPUTE UPPER CELL PROBABILITY 

1050 

XL-AT 

1051 

XU-1 .E10 

105a 

M-46 

1053 

FE-PN*PZC ( XL, XU, AM, AS, UA, UB, UC, LD, M ) 

1054 

FA-JCEL(NG) 

1055 

CS=CS+( (FA-FE )»(FA-FE ) )/FE 

1056 

NL-NG-1 

1057 

XL-AB 

1058 

M»1 

1059 

DO 50 J-a,NL 

1060 

XU-XL+DU 

1061 

FE=PN*PZC(XL,XU,AM,AS,UA,UB,UC,LD,M) 

1068 

FA-JCEL(J) 

1063 

CS=CS+( (FA-FE )*(FA-FE ) )/FE 

1064 

XL-XU 

1065 

TEST CHI-SQUARED UALUE FOR LEUEL OF SIGNIFICANCE 

1066 

F-NG-NPARM 

1067 

B = a./(9.:tcF ) 

1068 

ZT- ( 1 .-B-(CS/F )*)t:0. 3333333333 )/SQRT(B ) 

1069 



n»i 1070 

CALL TBLN(CR,ZT,Z,C,47,n) 1071 
RETURN 1072 
END 1073 
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4cDECK CYCLE ivfn 

SUBROUTINE CYCLE ( THEC, GR, DR,H, T1 , T2, DC, A, NDOUER, FCT ) 107S 

SUBROUTINE FOR COMPUTING THE LIMIT CYCLE 1076 

THEC-CONTROL ACCELERATION 1077 

GR-RATE TO DISPLACEMENT GAIN RATIO 1078 

DR-DEADBAND HALFUIDTH 1®79 

H-HYSTERSIS 

Tl-EFFECTIUE TURN ON DELAY TIME 1081 

Ta=EFFECTIUE TURN OFF DELAY TIME 1083 

DC-DUTY CYCLE 1®83 

A-LIMIT CYCLE RATE 1®84 

NDOUER-0 FOR NO DEADBAND OUERSHOOT 1085 

-1 FOR DEADBAND OUERSHOOT CONDITION 1086 

IF FCT IS GREATER THAN ZERO, PULSE UIDTH IS ESTIMATED AND USED 1087 
TO ESTIMATE EFFECTIUE CONTROL ACCELERATION 1088 

NDOUER-0 
FK-1.0 
D-DR/57.3 
TDD-ABS(THEC) 

K-0 1093 

IF (FCT .LE. 0 ) GO TO 10 

1095 

10 A=<TDD*T2*(GR-Ta/3. )+D*H)/(a.*GR-Tl-Ta) 1096 

ALPHA=D*(a.*GR-Tl-Ta-H*(GR-Tl/a. ) )/(Ta*(GR-Ta/a. )*(GR-Tl/a. ) ) 1097 

IF (TDD-ALPHA .LE. 0 ) GO TO 30 1098 

NDOUER-1 

C DUTY CYCLE APPROXIMATION FOR DEADBAND OUERSHOOT 1100 

H-H<^3. 1101 

AD-a.*D*(l .0-H)/(TDD:»:GR) H02 

DN0M-(Tl-T3+AD)/<a.*Ta-AD+( ( AD/2.+Ta-AD )*(T1-T2+AD )/GR ) ) 1103 

DC-1 ./( 1 .+DNOM) 1104 

GO TO 30 1105 

C COMPUTE DUTY CYCLE UITHOUT DEADBAND OUERSHOOT 1106 

E0 DC-FK/( 1 .+TDD5KD/A+T1-GR )/A ) 1107 



IF (K .LE, 0 ) GO TO 30 1108 
PU-2.*A/TDD 1109 
FK-1.0+FCT/PU 1110 
TDD=ABS(THEC*FK ) 1111 
K-0 Ilia 


P) V Ui 
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*DECK HISTO 1116 

SUBROUTINE HISTO ( X, N, XI1AX, XMIN, DX ) 1117 

THIS SUBROUTINE SORTS AND PLOTS A HISTOGRAM ON THE LINE PRINTER 1118 
X-ARRAV OF SAMPLES 1119 

N-NUMBER OF SAMPLES IN ARRAY 1130 

XMAX-MAXIMUM UALUE 1121 

XMIN-=MINIMUM UALUE 1122 

DX- INCREMENT FOR EACH CELL 1123 

K-NUMBER OF SAMPLES IN EACH CELL 1124 

DIMENSION X(1),KAXIS(94) 1125 

DATA ISTAR,IBLANK/1H*,1H / 1126 

C PRINT NUMBER OF SAMPLES, MIN AND MAX UALUES 1127 

URITE( 6,110) N,X(1),X(N) 1128 

L-1 1129 

IC»0 1130 

AXl-XMIN 1131 

> 10 AX2-AX1+DX 1132 

^ DO 20 J»L,N 1133 

IF (X(J)-AX2 ,GE, 0 ) GO TO 30 1134 

K-K+1 1135 

20 CONTINUE 1136 

30 IF (K-94 .GE, 0 ) GO TO 50 1137 

Kl-K+1 1138 

DO 40 J“K1,94 ' 1139 

40 KAXIS(J ) = IBLANK 1140 

K2»K 1141 

IF (K .LE. 0 ) GO TO 80 1142 

GO TO 60 1143 

50 K2-94 1144 

60 DO 70 J-1,K2 1145 

70 KAXIS(J)-ISTAR 1146 

80 URITE( 6,100) AX1,AX2,K, (KAXIS(J),J-1,94) 1147 

IF (AX2-XMAX .GE. 0 ) GO TO 90 1148 

L-L+K 1149 
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K-C 1150 

AXl-AXa 1151 

GO TO 10 115a 

90 RETURN 1153 

100 FORNAT <lX,F6,a,:|c TO*, F7. a, 15, « I*,94A1,*I») 1154 

110 FORMAT (1X,*N0 OF SAMPLES-*, 15, * MIN-*,F6.a,» MAX-*,F7.a,/4X,*RA1155 

INGE NUMBER 0 5 10 15 aO a5 30 35 40 45 1156 

a 50 55 60 65 70 75 80 85 90 95*, /a3X, *I . , . . I . , . , 1 157 

31. . ..1. ..1*. ..I. ...I I. ...I. ...I. ...I. ...I. ...I.. . .1. . ..11158 

4. . ..1. ...1. ...1. ...1*) 1159 

END 1160 
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:|:DECK OURF 

THIS SUBROUTINE COMPUTES A FACTOR(FAC) ON CONTROL MOTOR 116S 

THRUST FOR SHORT PULSES OF AN OUERSHOOT RESPONSE 1163 

T » TIME FROM ZERO THRUST TO FIRST OUERSHOOT PEAK 1164 

OR- THRUST OUERSHOOT RATIO FMAX/FPEAK 1165 

SUBROUTINE OURF (T, OR, FAC) 1166 

A-0.31831*ALOG(1.0/(OR-1.0) ) 1167 

P-ATAN(1.0/A) 1168 

U-6.a832/T 1169 

FAC- (1. 570796-0. 5*P-SIN( a. :t:P ) )/U 1170 

RETURN ll'^l 

END H'^2 
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*DECK PAGEHD 1173 

SUBROUTINE PAGEHD ( NRUN, NPAGE, NLINE ) 1174 

THIS SUBROUTINE EJECTS AN OUTPUT PAGE AND PRINTS RUN NO* 1175 

AND PAGE NO. AT THE TOP OF EACH NEU PAGE 1176 

URITE( 6,10) 1177 

URITE( 6,E0) NRUN,NPAGE 1178 

NLINE=5 1179 

RETURN 1180 

10 FORMAT (IHl ) 1181 

80 FORMAT (5X,7HRUN NO . , 15, 54X, 8HPAGE NO*, 15) 1188 

END 1183 
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*DECK PZC 

FUNCTION PZC(XL,XU, AM, AS,UA,UB,UC,LD,n) 

CONNON/CZ/Z < 47 ) , C ( 47 ) 

THIS FUNCTION SUBPROGRAM COMPUTES THE PROBABILITV OF A UALUE 
BETUEEN A LOUER AND UPPER LIMIT FOR A NORMAL OF UEIBULL 
DISTRIBUTION 
IF(LD.GT,0) GO TO 10 
C NORMAL DISTRIBUTION PROBABILITY 

SF=(XL-AM)/'AS 
CALL TBLN(P1,SF,Z,C,47,M) 

SF-(XU-AM)/AS 

CALL TBLN(P2,SF,Z,C,47,M) 

PZC-P2-P1 

RETURN 

C UEIBULL DISTRIBUTION PROBABILITY 

10 Pl-0. 

IF(XL-UA) 20,20,15 
15 AG=((XL-UA)/UB)**UC 
Pl=l .-EXP(-AG) 

20 P2*0. 

IF(XU-UA) 30,30,25 
25 AG=( (XU-UA)/UB )**UC 
P2 = l . 

IF(AG.GT. 100 ) GO TO 30 
P2=1.-EXP(-AG) 

30 PZC-P2-P1 
RETURN 
END 
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*DECK RANGE lEia 

SUBROUTINE RANGE (X,N,XNAX,XriIN,DX,NDX,ERRTB ) 1313 

THIS SUBROUTINE DIUIDES THE RANGE OF A TABLE OF NUMBERS INTO 1314 

N OR GREATER EQUAL PARTS UITH ESTHETICALLV PLEASING SCALES 1315 

TEST FOR INCREASING OR DECREASING FUNCTION. IF DECREASING 1316 

THE TABLE IS CHANGED TO AN INCREASING ARRANGEMENT 1317 

X - INPUT ARRAY 1318 

N - NUMBER OF UALUES IN ARRAY 1319 

XMAX= UPPER END OF LAST CELL 1330 

XMIN- LOUEST UALUE 1331 

DIMENSION X(l) 1333 

ERRTB=0 1333 

C TEST FOR ASCENDING OR DESCENDING ARRAY 1334 

IF (X(l)-X(N)) 40.10,30 1335 

10 ERRTB=1 1336 

RETURN 1337 

4, C DECREASING ARRAY, CHANGE TO ASCENDING ORDER 1338 

^ 30 LN-N/3 1339 

DO 30 J=1,LN 1330 

Z=X(J) 1331 

NF-N+l-J 1333 

X(J)=X(NF) 1333 

30 X(NF)»Z 1334 

C COMPUTE AESTHETICALLY PLEASING INCREMENTS OUER RANGE 1335 

40 XMIN=X(1) 1336 

XMAX=X(N) 1337 

DELT=XMAX-XMIN 1338 

ZDX»NDX 1339 

XINC-DELT/ZDX 1340 

XLIN«ALOG10(XINC) 1341 

LIN-XLIN 1843 

IF (XLIN .GE. 0 ) GO TO 50 1343 

LE--LIN+1 1844 

GO TO 60 1345 
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C 

50 LE--LIN 

60 niNC-=XINC*(10.»*LE) 

IF (MINC-5 .LT. 0 ) GO TO 70 
XINC“5,/(10.**LE) 

GO TO 100 

70 IF (niNC-4 ,LT* 0 ) GO TO 80 
XINC-4./(10.**LE) 

GO TO 100 

80 IF cniNC-S .GE. 0 ) GO TO 90 
XINC=1./(10.4C3»:LE) 

GO TO 100 

90 XINC-a./(10.*^^LE) 

100 DX-XINC 

C COMPUTE START AND FINISH XMIN AND XMAX 
R-XniN/XiNC 
KAR-ABS(R) 

k:r=r 

XKAR*KAR 

C TEST FOR FIRST CELL STARTING UALUE 

IF (KR .NE. 0 ) GO TO 130 
IF (R .LT. 0 ) GO TO 120 
110 XniN=0. 

GO TO 160 
120 XMIN*-XINC 
GO TO 160 

C COMPUTE NON-ZERO MINIMUM STARTING CELL LOCATION 

130 IF (R) 140,110,150 
140 XMIN-=-(XKAR+l . )*XINC 
GO TO 160 

150 XMIN-XINC*(XKAR-1 . ) 

C END CELL LOCATION 

160 XMAX-X(N)+XINC 
RETURN 
END 
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*DECK RNDX 1283 

SUBROUTINE RNDX (JK) 1283 

SUBROUTINE FOR COMPUTING THE RANDOM NORMAL DEUIATE 1284 

RANF IS CDC BUILT-IN PSUEDO-RANDOM NUMBER GENERATOR 1285 

K IS THE CURRENT RANDOM SEQUENCE INTEGER UHICH CHANGES 1286 

UITH EACH NEU RANDOM NUMBER GENERATION 1287 

COMMON K,RXD(30) 1288 

DO 30 I-1,JK 1289 

Q-RANF(K) 1290 

C THE FOLLOUING CODING UILL GIUE A SUITABLE REPLACEMENT FOR RANF 1291 

C KA>KX23 1292 

C JJ-1.0E+7 1293 

C N=KA/JJ 1294 

C M-N*JJ 1295 

C K»KA-M-N 1296 

C AK»K 1397 

> C Q=1.0E-7*AK 1398 

i C IN THE ABOUE OPTION K SHOULD BE A 7 DIGIT NUMBER 1399 

^ C TRANSFORM 0-1 UNIFORM DISTRIBUTION TO -0.5 TO 0.5 UNIFORM 1300 

IF (Q-.5 .LE. 0 ) GO TO 10 1301 

Q-1.0-Q 1303 

AB--1.0 1303 

GO TO 30 1304 

10 AB>>1.0 1305 

C TRANSFORM UNIFORM TO GAUSSIAN N<0,1) 1306 

30 XPT“-3.*AL0G(Q) 1307 

XRT-ABS(XPT) 1308 

Y-SQRT(XRT) 1309 

D-3.515517+.803853:»:V+.010338*Y»Y 1310 

E-1 .0+1 .433788*Y+. 189369*Y*Y+.001308»Y*Y*Y 1311 

C COMPLETE TRANSFORMATION TO GET RANDOM NORMAL DEUIATE 1313 

30 RXD(I )-AB*(Y-(D/E) ) 1313 

RETURN 1314 

END 1315 
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»DECK RUNQE 1316 

SUBROUTINE RUNGE (L, I) 1317 

THIS SUBROUTINE PERFORMS THE RUNGE -KUTTA INTEGRATION IN 1318 

COMBINATION UITH THE CALCU SUBROUTINE 1319 

DIMENSION RAT(4 ),SAU(4) 1320 

COMMON/DDV/YDOT <4),V(4),T,DT 1321 

I-I+l 1322 

GO TO (10,20,40,60,80), I 1323 

10 L-1 1324 

RETURN 1325 

20 DO 30 J-1,3 1326 

SAU(J)-V<J) 1327 

RAT(J)-VDOT(J) 1328 

30 Y(J)-SAU(J)+0.5*DT*VDOT(J) 1329 

T-T+0.5*DT 1330 

• L»1 1331 

^ RETURN 1332 

I 40 DO 50 J-1,3 1333 

RAT(J)-RAT(J)+2.*YD0T(J) 1334 

50 Y<J)-SAU(J)+0.5XDT*YDOT(J) 1335 

L-l 1336 

RETURN 1337 

60 DO 70 J-1,3 1338 

RAT(J)-RAT(J)+2.*YD0T(J) 1339 

70 Y<J)-SAU(J)+DT»YDOT(J) 1340 

T-T+0.5*DT 1341 

L-l 1342 

RETURN 1343 

80 DO 90 J-1,3 1344 

90 Y(J)-SAU(J)+<DT/6. )*<RAT(J)+YDOT<J)) 1345 

L-2 1346 

1-0 1347 

RETURN 1348 

END 1349 
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»DECK TABGEN 

SUBROUTINE TABGEN (Q0,U0,GAH0) 

THIS SUBROUTINE PROPAGATES THE TRAJECTORY PARAMETERS 
VELOCITYiDVNAPlIC PRESSURE, GANNA, AND ALTITUDE DURING 
A GRAUITY TURN BOOST. INITIAL UELOCITY AND DYNAMIC 
PRESSURE ARE USED TO DETERMINE ALTITUDE. THRUST, UEIGHT 
DRAG COEFFICIENT AND INITIAL FLIGHT PATH ANGLE(GANMA) ARE 
USED TO SOLUE THE EQUATIONS OF NOTION USING THE RUNGE 
AND CALCU SUBROUTINES. RETURNED TABLES ARE UIA COHHON DOUT 
DIMENSION DUN(90) 

CONNON/DOUT/NQT,TBL1(90),NUT,TBL5(90),NGT,TBLG(90),NHT, 

1 TBL11(90),TF,NTC,TBL4(90) 

C0NN0N/DRH/RH0(18),NRH 

COMNON/DDI/TBLa(90),TBL3(90),NTHT,NUT,CDS,Ma,M3,UBO 
C0NM0N/DDY/YD0T(4),Y(4),T,DT 
DATA DT,NP/0.5,4/ 

RHO IS A TABLE OF Z,P,Z,P,Z,P UALUES UHERE 
Z » ALTITUDE IN FEET 

P « NATURAL LOG OF INUERSE ATMOSPHERIC DENSITY(SLUGS/FT3 ) 
ALOGd/DENSITY) 

DATA (RHO( J),J« 1,18 )/0., 6. 04a, 100000., 10. 3417, 130000., 

1 11. 871, 140000., 13.1335, 160000., 18.9839,800000., 14. 3163, 

8 840000., 15. 8883, 388000., 20. 665, 557700., 87. 856/ 

DO 10 J»l,4 
10 YDOT(J)»0. 

MRH-1 
H2-1 
N3-1 

C COMPUTE INITIAL ALTITUDE FROM UELOCITY AND DYNAMIC PRESSURE 
CALL ALTI (ALTO,UO,QO) 

T-0. 

KP-NP 

C SET INITIAL TRAJECTORY STATES 
Y(1 )-UO 
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1371 

1372 

1373 

1374 

1375 

1376 

1377 

1378 

1379 

1380 

1381 

1382 

1383 




V(2)-ALT0 

1384 


Y(3)-GAnO 

1385 


Y(4)-Q0 

1386 


BEGIN INTEGRATION LOOP AND FILLING OF TABLES 

1387 

80 

IF (KP-NP .LT. 0 ) GO TO 30 

1388 


NQT-NQT+a 

1389 


NTQl-NQT-1 

1390 


TBLICNTQI )-T 

1391 


TBLSCNTQl )-T 

1392 


TBL6(NTQ1 )=T 

1393 


TBLIKNTQI )-T 

1394 


TBL11(NQT)-Y(2) 

1395 


TBL6(NQT)«=Y(3) 

1396 


TBL5(NQT)»Y(1) 

1397 


TBL1(NQT)»Y(4) 

1398 


KP-0 

1399 

30 

KP-KP+1 

1400 


I *0 

1401 


CALL RUNGE-KUTTA INTEGRATION SUBROUTINE 

1402 

40 

CALL RUNGE (L,I) 

1403 


IF (L.EQ.2) GO TO 50 

1404 


CALL STATE EQUATION UPDATE SUBROUTINE 

1405 


CALL CALCU 

1406 


GO TO 40 

1407 

50 

CONTINUE 

1408 


IF (T-TF .LE. 0 ) GO TO 20 

1409 


END OF INTEGRATION LOOP 

1410 


NUT-NQT 

1411 


NGT-NQT 

1412 


NHT-NQT 

1413 


ND-1 

1414 


CHANGE INPUT TABLE OF THETA-GAflMA TO NEU THETA TABLE 

1415 


DO 60 J-2,NGT,2 

1416 


T-TBL6( J-1 ) 

1417 






CALL TBLU < NTC, TH> T, TBL4 . flD ) 

1418 


DUn(J-l )-T 

1419 

60 

DUM(J )»TH+TBL6(J ) 

1480 


DO 70 J-1,NGT 

1481 

70 

TBL4< J )-DUrU J ) 

1488 


NTC-NGT 

1483 


RETURN 

1484 


END 

1485 
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*DECK TBLN 1436 

SUBROUTINE TBLN ( X, T, A, NT, N) 1437 

THIS SUBROUTINE IS A TABLE LOOKUP FROM ABSCISSA TABLE T 1438 

AND ORDINATE TABLE A . Y IS ORDINATE AT GIUEN ABSCISSA X. 1439 

NT IS LENGTH OF TABLES T AND A. MIS LOCATION OF LAST UALUE S0UGHT1430 

DIMENSION T(1),A<1) 1431 

10 IF (T(M)-X) 50,30,30 1433 

30 Y“A(M) 1433 

RETURN 1434 

30 IF (T(l)-X .LT. 0 ) GO TO 40 1435 

1*1“ 1 1436 

GO TO 30 1437 

40 M-M-1 1438 

GO TO 10 1439 

50 MM=M+1 1440 

IF (MM-NT .LE. 0 ) GO TO 60 1441 

M=NT 1443 

GO TO 30 1443 

60 IF (T(MM)-X .GT. 0 ) GO TO 70 1444 

M“MM 1445 

GO TO 50 1446 

70 M»MM-1 1447 

DT*T(MM)-T(M) 1448 

IF (DT .NE. 0 ) GO TO 80 1449 

Y“A(M) 1450 

RETURN 1451 

80 DY=A(MM)-A(M) 1453 

DDT“X-T(M) 1453 

Y“A(M)+DY)tcDDT/DT 1454 

RETURN 1455 

END 1456 
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:<cDECK 

TBLU 

1457 



SUBROUTINE TBLU (NT,V,X,T,n) 

1458 

C 


SINGLE TABLE LOOKUP SUBROUTINE 

1459 

C 


NT - NUMBER OF UALUES IN ARRAY 

1460 

C 


Y - RETURNED ORDINATE 

1461 

C 


X » ABSCISSA UALUE CALLED 

1468 

C 


T - INPUT TABLE OF ALTERNATING ABSCISSAS AND ORDINATES 

1463 

C 


ORDINATES MUST BE MONOTONICALLY INCREASING 

1464 

C 


n - PREUIOUS INDEX USED IN THIS TABLE LOOKUP 

1465 

C 


THIS INDEX GETS CHANGED TO CURRENT UALUE 

1466 



DIMENSION T(1 ) 

1467 


10 

IF (T(M)-X) 50,80,30 

1468 


80 

Y»T(M+1 ) 

1469 



RETURN 

1470 


30 

IF (T(l)-X *LT. 0 ) GO TO 40 

1471 



M-1 

1478 



GO TO 80 

1473 


40 

M«M-8 

1474 



GO TO 10 

1475 


50 

MM=n+8 

1476 



IF (MM-NT-1 .LE. 0 ) GO TO 60 

1477 



n-NT-1 

1478 



GO TO 80 

1479 


60 

IF (T(MM)-X .GT. 0 ) GO TO 70 

1480 



n»MM 

1481 



GO TO 50 

1488 


70 

M=nM-8 

1483 



DT-T(MM)-T(M) 

1484 



IF (DT .NE. 0 ) GO TO 80 

1485 



Y-T(M+1 ) 

1486 



RETURN 

1487 


80 

DY-T(MM+1 )-T(f1+l ) 

1488 



DDT-X-T(M) 

1489 



Y-T(M+1 )+DY*DDT/DT 

1490 



RETURN 

1491 



END 

1498 
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SUBROUTINE THENIN ( THEC, GR> DR, H,T1 > TE^THIIIN ) 

C SUBROUTINE FOR COMPUTING THE MINIMUM DISTURBING ACCELERATION 

C RESULTING IN ONE SIDED LIMIT CYCLE MOTION 

C THEC - CONTROL ACCELERATION 

C GR - RATE TO DISPLACEMENT GAIN RATIO 

C DR - DEADBAND HALFUIDTH 

C H - HYSTERESIS 

C T1 - EFFECTIUE TURN-ON DELAY TIME 

C T3 » EFFECTIUE TURN-OFF DELAY TIME 

C THMIN= MINIMUM ACCELERATION (RETURNED) 


THEC - CONTROL ACCELERATION 

GR - RATE TO DISPLACEMENT 

DR - DEADBAND HALFUIDTH 

H - HYSTERESIS 

T1 - EFFECTIUE TURN-ON DEI 

T3 » EFFECTIUE TURN-OFF DE 

THMIN= MINIMUM ACCELERATION 
D-DR/57,3 
A»a.i:GR-Tl-Ta 
B»(GR-Tl/a. )*T1 
C=(GR-Ta/2. 'itTZ 
Z-GR-Tl 
ABC“(B-C)/A 
CA-C/A 

Z2»(D*H+C*THEC)/A 
Zaa“(D*H-C*THEC)/A 
APR* B-Z*ABC- . 5* ( GR*GR+ ABC^ABC ) 
BPR=a.*D-Z*Z2-ABC*CA*THEC 
CPR- ( .5*Zaa:*:Zaa+ABC*D*H/A )/APR 
XPR=.5*BPR/APR 
ZPR= ABS ( XPR*XPR+CPR ) 

ZPRa»SQRT(ZPR) 

THMINl«=ABS(XPR+ZPRa) 
THMINa«ABS(XPR-ZPRE) 

IF (THMINl-THMINa . GE . 0 ) GO TO 10 
THMIN-THMINl 
GO TO 20 
10 THMIN-THMIN2 
20 RETURN 
END 
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515 

516 

517 

518 

519 

520 
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:t:DECK UBL 1527 

SUBROUTINE UBL (N,Y,NG,F) 1528 

THIS SUBROUTINE COMPUTES THE STATISTICAL DISTRIBUTION PARAMETERS 1529 
MEAN, STANDARD DEUIATION, SKEUNESS, KURTOSIS, GOODNESS 1530 

OF FIT TO NORMAL, LOG-NORMAL, AND UEIBULL 1531 

DIMENSION Y( 1),Q(20),R( 18 ),S(24 ),X(1000),F(20) 1532 

DATA <0(1 ),I = l,20)/-0.95,3.46, 1533 

1-0.87,3.0,-0.64,2.3,-0,37, 1.8,-0.09,1.39,0.16,1.1,0.63,0.7,0.93, 1534 

20.5,1.52,0.2,2.0,0.001/ 1535 

DATA (R(I ), 1 = 1, 18)/0., -0.001, 0.4, -0.95, 0.8, 1536 

1-1. 7, 1.2, -2. 43, 1.8, -3. 48, 2. 4, -4. 55, 3. 0,-5. 65, 3. 66, -6. 92, 4. 5, -8. 5/ 1537 
DATA (SCI ), I -1,24 )/0. , 1., 0.2, 0.932, 0.4, 0.902, 0.6, 0.888, 0.8, 0.886 1538 

1 ,1.1,0.893,1.8,0.9275,2.2,0.947,2.6,0.963,3.0,0.974,3.5, 1539 

2 0.984,4.5,0.994/ 1540 

MQ»1 1541 

MR-1 1542 

MS-1 1543 

AUI-N-1 1544 

XK-1./SQRT(1.-1.96*SQRT(2./AUI ) ) 1545 

NQ-20 1546 

NR-18 1547 

NS-24 1548 

AN-N 1549 

S1=0. 1550 

52- 0. 1551 

53- 0. 1552 

54- 0. 1553 

DO 10 J»1,N 1554 

X(J)-Y(J) 1555 

10 S1«S1+X< J ) 1556 

C COMPUTE THE MEAN UALUE 1557 

XM-Sl/AN 1558 

DO 20 J-1,N 1559 

DX-(X(J)-XM) 1560 
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DS-DX*DX 1561 

Sa-S2+DS 1563 

S3-S3+DS*DX 1563 

30 S4-S4+DS*DS 1564 

C COMPUTE THE UARIANCE, SKEUNESS,KURTOSIS, AND STANDARD DEUIATION 1565 

UAR-S3/AN 1566 

SKEU»(S3/AN)/(UAR2|:*1.5) 15^7 

SKUR- (S4/AN)/(UAR*UAR) 1568 

SIG = SQRT(UAR ) 1569 

C COMPUTE THE A,B,C PARAMETERS FOR A UEIBULL DISTRIBUTION FUNCTION 1570 

CALL TBLU ( NO, CL,SKEU, Q, MQ ) 1571 

C-EXP(CL) 1573 

CALL TBLU ( NR, BL, CL, R, MR ) 1573 

B»SQRT(UAR/EXP(BL) ) 1574 

CALL TBLU ( NS, A1 , CL, S, MS ) 1575 

A-XM-B*A1 1576 

OC-l./C 1577 

F(1)-XM 1578 

F(3)=SIG 1579 

F(3)=SKEU 1580 

F(4)-SKUR 1581 

F(9)»XM+XK*3.576*SIG 1583 

F(ll)=A 1583 

F(13)=B 1584 

F(13)*C 1585 

F( 14 )»A+B*(0. 694**00 1586 

F( 15 )«A+B*(4. 6**00 1587 

F(16)“A+B*(5. 3**00 1588 

F(17)-A+B*(6,9**0O 1589 

USE CHI-SQUARED GOODNESS OF FIT TEST TO ESTIMATE LEUEL OF 1590 

SIGNIFICANCE (CR) 1591 

TEST UEIBULL DISTRIBUTION 1593 

CALL CHISQ (N,X,CR,NG,XM,SIG,A,B,C,1 ) 1593 

F(30)-CR 1594 
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C TEST NORMAL DISTRIBUTION 1595 

CALL CHISQ (N,X,CR,NG,XM,SIG,A,B,C,0) 1596 

F(18)-CR 1597 

Sl-0. 1598 

Sa-0. 1599 

S3-0. 1600 

S4>0. 1601 

C COMPUTE THE LOGARITHM OF UALUES 1603 

DO 30 J-1,N 1603 

X( J)=ALOG(V( J) ) 1604 

30 S1=S1+X(J) 1605 

XM-Sl/AN 1606 

DO 40 J-1,N 1607 

DX=(X(J)-XM) 1608 

DS = DX3t:DX 1609 

SB-S2+DS 1610 

S3-S3+DS*DX 1611 

40 S4=S4+DS*DS 1612 

UAR-S2/AN 1613 

SIG»SQRT(UAR) 1614 

SKEU=S3/(AN*UAR**1.5) 1615 

SKUR»S4/(AN*UAR^:UAR) 1616 

F(5)=XM 1617 

F(6)=SIG 1618 

F(7)“SKEU 1619 

F<8)“SKUR 1620 

F(10)-EXP(XM+XK*2.576*SIG) 1621 

C TEST LOG-NORMAL DISTRIBUTION FOR LEUEL OF SIGNIFICANCE 1622 

CALL CHISQ (N,X,CR,NG,XM,SIG,A,B,C,0) 1623 

F(19)-CR 1624 

RETURN 1625 

END 1626 
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